THE JOURNAL 


OF 


CHEMICAL PHYSICS 


8 


MARCH, 1940 


Numser 3 


The Infra-Red Absorption Spectrum of Methyl Amine 
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(Received December 27, 1939) 


The parallel type methyl amine band at 10u has been resolved. Its center is at 1045.27 cm™, 
and the very regular structure of its P and R branches indicates a close approximation to axial 
symmetry. The mean of the two largest moments of inertia is 37.43 X 10~* g-cm*. Assuming 
angles and distances for NH: as in ammonia, and for CH; as in methane, the C—N distance is 
computed to be 1.48 X 10-* cm. Certain enhanced lines suggest a superposed perpendicular band. 
The 13, band is partially resolved, and assigned to a motion corresponding to v3 in NH3;. The 
bands arising from vibrations characteristic of the amino and methyl groups have been ex- 


amined but could not be resolved. 


BSERVATIONS upon the infra-red spec- 

trum of methylamine vapor which have 
recently been published by Cleaves and Plyler! 
include the envelopes of a number of unresolved 
absorption bands, the frequencies of which have 
been assigned to various normal modes of vibra- 
tion of the CHsNHe molecule. The same spec- 
trum has been under investigation in this labo- 
ratory with a slightly higher dispersion, the 
results being identical in some respects, but 
sufficiently different in others to justify further 
discussion. 

The molecule approximates a rigid rotator 
which is almost, but not quite, axially sym- 
metrical, and its spectrum is similar in many 
respects to those of CH;OH and CH;F. It may 
also be thought of as an ammonia derivative in 
which the methyl group has replaced one of the 
hydrogen atoms. There is a single plane of 
symmetry through the CN axis and bisecting the 
HNH angle. The motions associated with the 
fifteen vibrational degrees of freedom may be 
subdivided for convenience into three classes. 


Cleaves and Plyler, J. Chem. Phys. 7, 563 (1939). 


Six of them represent the mutual oscillations of 
the methyl and amino groups; three others are 
characteristic of the amino group alone, while the 
remaining six belong primarily to the methyl 
radical. 


MUTUAL VIBRATIONS OF THE METHYL AND 
AMINO GROUPS 


The molecule may first be considered as con- 
sisting of two approximately rigid groups of 
atoms attached to one another through the C—N 
bond. Of the six characteristic motions of this 
structure (rotation and translation excluded) the 
simplest is the valence oscillation, which is 
accompanied by a change in electric moment 
nearly parallel to the CN axis. The corresponding 
absorption band lies close to 10u, and, as 
represented by Cleaves and Plyler, has a narrow 
intense zero branch. Under somewhat higher 
resolution the rotational structure of this band 
appears, showing the remarkably regular pattern 
indicated in Fig. 1. The spectral range included in 
the slit was about 0.5 cm~'. The zero branch 
remains quite narrow and unresolved, and the 
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Fig. 1. The band at 10x. 


individual lines composing the wings show no 
observable splitting. This indicates that axial 
symmetry is closely approximated. Table I 
presents the positions of the observed lines, in 
wave numbers, reduced to vacuum. 

An unusual feature of this band is the enhance- 

“ment in intensity of certain lines, particularly of 
+1, although +7, +15 and others stand out 
above the envelope, and a strong narrow line 
appears at 1085 cm~, at the high frequency end 
of the band. This suggests that a weak perpen- 
dicular type band with relatively wide spacing is 
superposed upon the principal absorption. Fur- 
ther observations with different amounts of gas 
may reveal more detail. 

A second mode of motion associated with the 
CN axis is the torsional oscillation. Cleaves and 
Plyler assign to it the frequency 620 cm™, and 
report a band envelope resembling that of a 
parallel vibration, with an intense though rather 
wide central maximum. Our observations seem 
hardly compatible with this. While an irregular 
pattern of rather weak absorption extends from 
500 to 700 cm-!, we have observed no region of 
particularly high intensity, and no characteristic 
band structure. It is of course possible that this 
absorption represents the fundamental band due 
to torsional oscillation, but we are inclined to 
believe that the frequency is too high* to be so 
interpreted, and that instead one or more differ- 
ence bands are involved. 

A pair of perpendicular vibrations are expected 
in which the C—N bond is deformed either in the 
plane of symmetry or normal to this plane. These 
are perhaps represented by rather weak bands at 

*This supposition is supported by observations of 
Randall and Lawson, as yet unpublished, but kindly com- 


municated to us. They find a region of rather strong ab- 
sorption in the neighborhood of 270 cm™. 


about 1130 cm and 1150 cm. Here the 
structure is quite complicated, with no very 
intense lines, and is not sufficiently resolved to 
permit an analysis. 

The other two vibrations characteristic of a 
molecule considered as a pair of approximately 
rigid groups are torsional oscillations in which 
the hydrogens belonging to NH2 move, either in 
the same sense or in opposite senses, along lines 
perpendicular to the N—H bond and parallel to 
the plane of symmetry. Thinking of the structure 
as a substituted ammonia, the first of these 
corresponds to the symmetrical NH; vibration 
v3, which gives rise to a pair of bands near 10x, 
the doubling being a consequence of the two 
minimum potential function. This interesting 
characteristic must persist in methyl amine, 
though somewhat modified by the fact that 
whenever the N atom passes through the 
H—H-—CHs; plane from one equilibrium position 
to the other, the methyl group must rotate 
through an angle of 60° about its own symmetry 
axis. The resulting absorption band extends from 
about 710 cm-! to about 850 cm-, its central 
portion having the structure indicated in Fig. 2. 


TABLE I. Observed wave numbers, reduced to vacuum, for lines 
in the 10u band. J for the lower rotation state. 


R BrANcH 


1061.52 
1062.77 
1064.08 
1065.39 
1066.70 
1067.90 
1069.10 
1070.24 
1071.38 
1072.59 
1073.73 


P BRANCH 


1026.23 
1024.56 
1022.90 
1021.19 
1019.54 
1017.84 
1016.10 
1014.35 
1012.56 
1010.82 
1009.06 


1046.44 
1047.86 
1049.30 
1050.71 
1052.11 
1053.45 
1054.85 
1056.21 
1057.56 
1058.86 
1060.22 
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ABSORPTION SPECTRUM OF METHYL AMINE 


Fic. 2. A part of the 


band. 


740 760 


There is a rather strong central region around 
780 cm—', obviously consisting of many lines, and 
flanked on either side by a fairly sharp maximum 
suggesting a Q branch. The positions of these 
maxima are 763 cm~ and 795 cm-'. The change 
in electric moment will include both parallel and 
perpendicular components. 

The corresponding unsymmetrical oscillation 
introduces an angular momentum which must be 
compensated by a tilting of the methyl group. 
Thus the change in electric moment is principally 
in a direction normal to the symmetry plane of 
the molecule. This nondegenerate perpendicular 
vibration corresponds to the 6u band, v4, of NH3. 
It has not been identified in the spectrum, but 
may possibly lie in the 10u region, or even as far 
out as 


VIBRATIONS CHARACTERISTIC OF NH, 


The symmetrical valence oscillation of the 
amino group should have both perpendicular and 
parallel components. The latter, though rela- 
tively less intense, is unmistakable because of its 
sharp Q branch, and the frequency may thus be 
determined. Our observations are so nearly 
identical with those presented by Cleaves and 
Plyler,* that it is not necessary to reproduce the 
curves. The zero branch is plainly indicated at 
3360 cm—. The absorption pattern is not sym- 
metrical with respect to this ordinate, however, 
but is much more intense on the high frequency 
side. Its appearance suggests a superposed per- 


100 200. +300 400 
Fic. 3. ApF/2(2J+1) vs. (J+4)2. 


* Reference 1, Fig. 6. 


800 cm 


pendicular band at about 3410 cm~ due to the 
unsymmetrical valence oscillation. 

The symmetrical deformation oscillation of 
the NH: group, which includes both perpen- 
dicular and parallel components, is represented 
by a band at 1625 cm~. Our observations show 
no more detail than those of Cleaves and Plyler, 
as represented in Fig. 4 of their report. 


VIBRATIONS OF THE METHYL RADICAL 


The symmetrical and unsymmetrical valence 
oscillations in CH; yield superposed bands in the 
neighborhood of 3.5u. A well-defined zero branch 
appears at 2820 cm, with a somewhat similar 
maximum at 2964 cm. These are also observed 
in Raman scattering as strong sharp lines.? 
Whether or not they are related through reso- 
nance with 2écu is not clear, but the question 
of resonance in this case is complicated by the 
lack of true axial symmetry or complete de- 
generacy. The center of the strong perpendicular 
type band appears to be around 2930 cm. 

The symmetrical and unsymmetrical deforma- 
tion vibrations produce bands in the neighbor- 
hood of 74, which we find very poorly resolved. 
There seems to be a weak zero branch at 1385 
cm~, and a broad complex absorption with its 
center near 1470 


MOLECULAR MOMENTS OF INERTIA 


When the molecule is considered as a sym- 
metrical rotator the term values may be written 
W/hce=G(v) + F(J,K,v), 

where 
F(J,K ,v) 
1)?+X 
The mean of the two larger and nearly equal 
moments of inertia for the ground state is equal 
to h/8x*cX,°, which may be evaluated by the 
usual combination relations. 
A2F°= R(J+1) —P(J—1) 
= 
? Kirby-Smith and Bonner, J. Chem. Phys. 7, 880 (1939). 
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Fig. 3 shows A:F/2(2J+1) plotted against 
(J+), the upper line corresponding to v=0, 
and the lower one to v=1. While the points 
scatter considerably, due in part to the superposi- 
tion of one or more weak bands, they approxi- 
mate straight lines with the constants 


X,°=0.7385 and X2°=12X10-’, 
X;'=0.7300 and X2!=36X10-7. 


The moment of inertia in the ground state is, 
therefore, 37.43X10-* g cm’, and the corre- 
sponding C—N distance is 1.48X10-* cm, as- 
suming angles and distances for NHe as in 
ammonia, and for CH; as in methane. 

The band center is at 1045.27 cm“. Its position 
is obtained by extrapolating to J=0 the mid- 
points between successive pairs of lines chosen 
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from the P and R branches. These midpoints are 
well represented by the equation 


b= vy + — 


The smallest moment of inertia would be 
obtained at once from the spacing between Q 
branches in any nondegenerate perpendicular 
band. If the enhanced lines in the 10u region 
constitute such a band, with a mean spacing of 
about 6.2 cm—, the value of C is about 7.2 X 10-*° 
g-cm?. This can hardly be considered as well 
established, however. 

The examination of methyl amine was sug- 
gested to us by Professor J. G. Aston, of Pennsyl- 
vania State College, and it is a pleasure to 
acknowledge our indebtedness to him for a supply 
of ammonia-free material. 
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Spectral Transmission of Glucose Glass in the Near Infra-Red 


VOLUME 8 


WiLL! M. 
El Cerrito—Berkeley, California 


(Received November 13, 1939) 


The transmission of samples of glucose glass varying in thermal history was previously 


investigated. This work covered the visible part of the spectrum. Transmission data in the near 
infra-red are presented now as far as could be obtained photographically (12,000A). No limit 
of absorption was found at long wave-lengths as far as tested. The possible use of organic 


I. THE METHOD 


HIS paper is a continuation of my paper.’ 
Two series of glucose glasses have been 
tested prepared by somewhat different methods, 
see paper I. The series include samples GG 1, 2; 
and 3, 5, 6, 8, 9, respectively. (Glasses GG 4 and 7 
were now omitted because of increased crystal- 
lization.) The samples range from clear to dark 
brown in color. The Hilger grating spectrograph 
of paper I was adapted for work in the near 
infra-red. The first order was used, the over- 
lapping ultraviolet of the second order being 
eliminated by means of a No. 88 Wratten filter. 
A neutral wedge equal in length to the length 
of the slit was inserted in front of the slit as 
discussed in paper I. The width of the slit 
varied from 0.1 to 0.3 mm. 


1 Willi M. Cohn, “A Spectral Investigation of Glucose 
Glass,” 4, Chem. Phys. 6, 65-67 (1938). (Cited as 


glass for investigations in the near infra-red is discussed. 


Wedge spectra were obtained on Eastman 
plates, emulsions K, P, M, and Z using sensi- 
tivities I or 144. Plates of series M and Z were 
hypersensitized by means of ammonia. A refer- 
ence spectrum of an incandescent lamp operated 
at a constant voltage, together with spectra of 
one or two glass samples were taken on a plate. 
The photometric reduction was carried out as 
discussed in paper I. 


II. THE RESULTS 


Transmission values D, as a function of the 
wave-length are presented in Table I. These 
data include results of both former and present 
investigations. We further computed the “‘rela- 
tive transmission” of the samples by arbitrarily 
setting the transmission at \6000 equal to 1.0. 
This permits an easier comparison of the two 
groups of samples. Fig. 1 shows the relative 
transmission as a function of the wave-length. 
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SPECTRAL TRANSMISSION OF GLUCOSE GLASS 
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Fic. 1. Relative spectral transmission of glucose glasses varying in thermal history. 


The results may be summarized as follows: 
Glucose glass shows, with increasing carameliza- 
tion, a general absorption over the entire 
spectrum; a rather flat maximum of trans- 
mission; no narrow bands of absorption; and a 
shift of the maximum of transmission with an 
increase in coloring matter. Glasses GG 5 and 8 
show secondary maxima of transmission at 8100 


TABLE I. Glucose glass; transmission D) in the visible and 
near infra-red parts of the spectrum; thickness 
of layer 25.0 mm. 
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and 8600A, respectively. None of the samples 
reaches the limit of transmission at long wave- 
lengths within the region investigated : Quantita- 
tive data could be determined from wedge 
spectra as far as 9800A, and the visual examina- 
tion of spectra on plates Z could be carried out 
as far as 12,000A. (Since wedge spectra on plates 
Z are of less contrast than the remaining spectra 
we consider the results of the visual determina- 
tions of less weight than the remaining data.) 

The comparatively high values of the trans- 
mission of glucose glass in the near infra-red 
seem of interest in connection with investigations 
in this region of the spectrum, particularly in 
regard to the ease of preparing containers of 
various shapes from this material. 


III. SUMMARY 


The presence of small amounts of colloidal 
carbon formed during caramelization causes char- 
acteristic changes in color and spectral trans- 
mission of glucose glass. All samples show 
comparatively high values of transmission in the 
infra-red as far as tested (12,000A). For a com- 
parison of the color and spectral transmission of 
glucose glass and inorganic glass prepared with 
the addition of carbon, see our paper.? My thanks 
to Dr. George S. Parks of Stanford University in 
connection with these investigations were ex- 
pressed in my former paper dealing with glucose 
glass. 

2 Willi M. Cohn, “Color and spectral transmission of 


some inorganic and organic glasses containing carbon.” 
To be published soon. 
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Intensities of Electronic Transitions in Molecular Spectra 


IX. Calculations on the Long Wave-Length Halogen Spectra* 


RoBeErRT S. MULLIKEN 
University of Chicago, Chicago, Illinois 
(Received December 8, 1939) 


The dipole strength Dyg for the perpendicular-type 
N—@ transitions in the long wave-length spectra of the 
halogen molecules is computed theoretically according to 
each of two approximations (LCAO and AO) for the wave 
functions. The work illustrates the behavior of antisym- 
metrical wave functions in transition moment calculations, 
and justifies the ignoring (as an approximation) of closed 
shells. The overlapping and dipole moment integrals 
S [o(npo) o(npo’) and S[¢(npo) are 
evaluated using Slater’s approximate atomic orbitals for 
the ¢’s. The calculated Dyg values for both approximations 
show the same rapid variation with atomic number as do 
the observed values, except for iodine, where the anomalous 
smallness of the observed value may reasonably be 
attributed to strong case c effects. The LCAO calculated 
values are consistently much larger than the AO values; 
the latter agree surprisingly well with the experimental 


values (cf. Table I). The calculations show Dyg to be very 
sensitive to variations of the internuclear distance 7, and to 
variation of the effective charge Z* in the Slater atomic 
orbitals used. It is pointed out that the former variation 
may require reconsideration of work in which repulsive 
potential energy curves are deduced from intensity data on 
absorption continua. A qualitative discussion is given of the 
effects on Dyg of various possible improvements of the 
wave functions used. In connection with this, calculations 
have been made which show that the energy of the oyls 
antibonding state of H.* in LCAO approximation is 
minimized when Z* for 1s is (roughly) 0.84, for normal 
internuclear distance. This contrasts with Z*=1.228 de- 
termined by Finkelstein and Horowitz for the og1s bonding 
state, and suggests that to improve molecular wave 
functions Z* should be increased for bonding but decreased 
for antibonding molecular orbitals. 


I. INTRODUCTION 


earlier papers,'? the possibility of 
making significant quantum-mechanical cal- 

culations of absolute intensities in molecular 
spectra has been explored, and its importance for 
the interpretation of experimental data has been 
emphasized. The calculations previously reported 
have been for transitions of the special very 
intense ‘‘N-—V” type. In diatomic molecules, 
N-V transitions are polarized parallel to the 
axis of the molecule. 

Calculations of the intensities of the N-V 
spectra of the halogens and comparisons with 
observation have already been given in one of 
the previous papers,' except that the calculation 
of the overlapping integral S has not been de- 
tailed. This is included in Section II of the 
present paper. 

The present paper is devoted to calculations 
on the weak perpendicular-type transitions 
N-—"'Q, of the long wave-length absorption 
spectra of the halogen molecules. The results of 
these calculations are compared with observed 

* Reported at the 1938 Chicago meeting of the American 
Physical Society; Phys. Rev. 55, 239(A) (1939). 


1R.S. Mulli mf Chem. Phys. 7, 14, 20 (1939). 
R.S, Mulliken, J. Chem. Phys. 7, 121, and later (1939). 


absolute dipole strengths.* These were evaluated 
in a separate paper.‘ 

The calculations are based on certain assumed 
electron configurations for states N and Q; these 
have been discussed in previous papers® and now 
seem to be well assured. The electron configura- 
tion of any molecular state can be expressed 
either in terms of AO’s (atomic orbitals) or of 
MO’s (molecular orbitals). These correspond to 
two different approximations to the electronic 
structure. Calculations are made below corre- 
sponding to both. 

Before proceeding further, it should be men- 
tioned that there really must be two Q complexes, 
i.e., two +3Q sets of states in the halogen mole- 
cules.> One comprises **II,, the other '*II,, 
states; it is predicted that the latter should be 
somewhat the higher in energy. According to the 

The comparison of observed with calculated dipole 
strengths D (cf. reference 1 for definitions) has advantages 
over the comparison of observed and calculated f values, 
which was used in previous papers of this series (references 
1, 2). Namely, a complete theoretical calculation can 
readily be made for each D, whereas for a purely theoretica 
f, the frequency of the transition would have to be calcu- 
lated in addition. In previous papers (cf. reference 2, 
footnote 2), this was avoided by using ‘‘semi-theoretical 

’s” combining calculated D’s with observed »’s. 


4R.S. Mulliken, Phys. Rev. (in press). 
5 R.S. Mulliken, Phys. Rev. 46, 549 (1934). 
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INTENSITIES IN MOLECULAR SPECTRA 


selection rules for homopolar molecules, transi- 
tions from the normal state N 'Z+, can occur 
only to the »*II,, set, which is then to be identified 
with the observed Q complex. Hence we have 
430, (observed), and '*Q, (predicted). Although 
only the former is of immediate interest, it is 
desirable for a systematic understanding to con- 


235 


sider both. This is true especially because in the 
mixed halogens (see X of this series), transitions 
to both Q complexes are allowed, and are ap- 
parently observed. 

The electron configurations and state-species 
of states NV, Q., and Q, are as follows (inner-shell 
electrons are omitted): 


N: (npo-npo')(npr)*(npr’)4, 1+, 


AO approximation: 


npo* npx® npo' npr’! 


10.510, 
\npo npr npo”? npr’® 


N: np) zt, 
MO approximation: {*'0,: rnp  ounp, * (1b) 


In (1a) and (1b), n=2, 3, 4, 5, for Fe, Cle, Bro, 
I,. The primed and unprimed symbols refer to 
the two atoms of the molecule. 

For the Q, and Q, states jointly, two equivalent 
electron configurations are given for the AO 
approximation in (1a). The actual Q, and Q, 
wave functions correspond to “resonance” be- 
tween these two configurations, that is they are 
two 50,50 linear combinations of wave functions 
corresponding to the two configurations. These 
wave functions are discussed at the beginning of 
Section ITI. 

If we use AO chemical bond terminology,‘ 
two three-electron bonds (one o, one 7), equiv- 
alent in strength to one (weak) ordinary bond, 
may be said to exist in the Q, states (both the 
*Q,, and the 'Q,). The Q, states may be said to 
have two three-electron antibonds, equivalent 
to one weak ordinary antibond. [The term 
“antibond,” not hitherto used, seems suitable. ] 

In MO chemical bond terminology, the Q, 
states (both *Q and !Q) as described by (1b) have 
two o (strong) and four 7 (weak) bonding elec- 
trons, and one o (strong) and three 7 (weak) 
antibonding electrons, giving a net bonding effect 
equivalent to one strong o and one weak x 
bonding electron, or two half-bonds, one o 
(strong) and one w (weak). Similarly in the Q, 
States, the net effect is equivalent to two half- 
antibonds, one a, one 7. 


‘L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, 1939), Chapter VIII. 


II. CALCULATION OF N->1°Q DIPOLE STRENGTHS, 
MO anv LCAO TueEory 


Referring to Eq. (18) of a related paper,' it is 
seen that the sum of the dipole strengths of 
N—'Q,; and N—'0Q in the halogen spectra should 
be equal to the quantity Dyg.’ This represents 
the dipole strength which N-—'Q, or N-'Q, 
would have for A,S coupling. Dye is equal to the 
square of Qyg, which we wish here to calculate 
quantum-mechanically. By Eq. (17) of the paper 
just mentioned, ‘ 


Ove= f Ho 


where the summation 2%; is over all the electrons 
in the molecule. To evaluate this integral, we 
proceed as in Sections 4 and 5 of a previous 
paper.! We shall begin with the MO approxima- 
tion. The discussion will be given in considerable 
detail, in order to make clear certain points that 
will arise repeatedly in the calculation of Q 
integrals. 

The antisymmetrical wave function y(N) 
corresponding to the MO electron configuration 
in (1b) is as follows 

7The antisymmetrical wave functions of molecules in 
terms of MO’s are constructed in the same way as are 
antisymmetrical wave functions of atoms in terms of 
AO’s. See J. C. Slater, Phys. Rev. 34, 1293 (1929); 38, 
1112-1115 (1931). 

§ The normalization factor (!)~4 is strictly correct onl 


if all the MO’s used are mutually orthogonal (cf. J. C. 
Slater, reference 7). 
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¥(N) = (n 1)? (3) 


where . . . stands for 1*2¢(3)2*2(4)rya(5)ya(6) 
Xtra(7)29(8)oa(9)o8(10) and so on down 
through all the closed shells and the x electrons. 
The single-electron wave function symbol 
m*ya(1) is an abbreviation for r,np,(1)a(1) ; and 
so on; my, mz, stand here for r,np,, 
TuNpy, Tunpz, and o,np, respectively, and 
a and 8 for the two possible states of the electron 
spin. P refers to any permutation of electrons 
among the single-electron wave functions. The 
quantity (—1)? is +1 for even, —1 for odd per- 
mutations. For the degenerate t MO’s we are 
here using real forms 7, and 7,: 7:=f(p,2) cos ¢, 
y= f(p,2) sin ¢, with axes x, y perpendicular to 
the molecule axis, p= (x?+y’)!, and ¢=tan7! y/x. 

The antisymmetrical MO configuration wave 
function for Yo('Q,) is® ° 


with 


where the dotted parts are the same as in Eq. 
(3). The symbol o* is an abbreviation for o,np. 
¥o(Qz) is just like Yo(4Q,) except that, in 
and y, x and y are exchanged in the three z* 
factors in the basic permutation. In passing, it 
may be noted that the ¥(°Q:) wave functions 
differ from the ¥o('Q) only in having a negative 
sign before y; in the first of Eqs. (4). 

Putting the results of Eq. (3) and (4) into Eq. 
(2), and reducing the double 2p summations to 
single ones in the manner given by Slater 
(reference 7, p. 1115), one finds 


(5) 


* For the very best MO gee y, analogous to 
the Hartree y of an atom (cf. Slater, reference 7), some of 
the MO’s would no longer be quite orthogonal. (Many are 
always orthogonal for reasons of symmetry.) However, the 
corrections then required in the normalization factor and in 
the value of Qvg from Eq. (5) would be relatively unim- 
portant, and may properly be neglected here. 
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If the MO’s in Eggs. (5) are all mutually orthog- 
onal, as may be assumed with negligible error,° 
it is found that all but two of the very numerous 
integrals contained in Eq. (5) vanish. These 
two are 


X [a*ya(1)o* (2) f 
and 


J 


Hence 


Qve=2! f (x onpy)y(ounp)dv 


f (6) 


This result is exactly the same as if from the 
beginning we had ignored all shells of electrons 
which remain closed during the transition. Eq. 
(6) is exactly the same as if we were dealing with 
a two-electron problem corresponding to a 
transition 


(x (oun) 
(x 


For if we treat this two-electron problem just 
as we treated the transition (¢,1s)?—>(o,1s)(ou1s) 
of He in a previous paper,! writing equations 
analogous to Eqs. (8) and (9) of that paper, we 
arrive at the present Eq. (6), which is analogous 
to Eq. (10) of that paper. 

To make as accurate as possible a calculation 
of Qve by Eq. (6), we should use the most 
accurate possible (Hartree-type) forms for the 
MO’s z,np and o,np in Eq. (6).° If we knew 
these, Qveg might turn out to be fairly accurate, 
but unfortunately we do not know them. For a 
rough calculation, however, we may use the 
well-known rough LCAO forms to approximate 
our MO’s. 

The normalized LCAO MO’s needed, in real 

10For the present purpose, 7,mpz and x,np, belong in 


effect to two different shells, of which one or the other 
remains closed during the transition. 
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form, are :51 
= (7) 
oyunp = (npo —npo’) 


where the overlapping integrals S,,, and S,,: are 


f 
(8) 
f (npo)(npo)dv. 
The AO’s npo and npo’ are here to be taken 
both to have the same sign in the region where they 
overlap in the middle of the molecule. 


Introducing the LCAO forms for z,”p and 
ounp (cf. Eqs. (7)), Eq. (6) reduces to 


—Qve= Sex) Soo’) Ore), (9) 


where 
= f (nprz)x(npo’)dv 

= f (10) 
f 


f (npx,)y(npo)dv=0. (11) 


Since Q,,=0, it remains to evaluate Q,.. For 
this purpose, we must introduce analytical forms 
for the AO’s np, and npc. Suitable are the radially 
nodeless forms given by Zener and Slater."' The 
lack of radial nodes can be seen to be of negligible 
importance for our calculations. Nevertheless, the 
use of these simple forms for the AO’s introduces 
an additional element of approximation into our 
calculations. 

If A and B refer to the two atoms, the Slater 
functions are of the forms 


npz= Nywxara™— exp [—hra], (12) 


npo’ = N exp [—hrz ], (13) 
where 
h=Z*/aon*, (14) 


4 J. C. Slater, Phys. Rev. 36, 57 (1930). 


with ad» equal to the Bohr radius of the normal 
hydrogen atom, and n*=2, 3, 3.7, 4, Z*=5.12, 
6.1, 7.6, and 7.6 for the respective atoms F, Cl, 
Br, and I.” The normalizing factors Ny» are 
readily determined to have the values (h5/r)!, 
(2h7/157)*, and (h°/1057)!, for F, Cl, and I, 
respectively ; the calculation for Br has not been 
made here. In Eqs. (12), (13), x4 and r, refer to 
the x or r coordinate measured from the center 
of the A atom, with similar meanings for zg and 

Taking the origin of coordinates at the center 
of the molecule, with the z axis directed from B 
to A, and letting r be the distance between the 
nuclei, we have zg=2+7r/2, x4=x, and so on. 
Introducing these relations into Eqs. (12), (13), 
and the latter into Eq. (11), we have 


f 49) (rare) 
Xexp (15) 
The integrations can be effected without 
obstacle after transformation to confocal elliptic 
coordinates 7, where £=(ra4+rz)/r, 
n=(ra—rg)/r. The results for Fe, Cle, and I, 
(Br2 has not been done) are as follows: 
Fe: (16) 


where 
y=rh; (17) 


Cle: = +y+34y2/75+3y3/25 
+31y'/1575+7°/525); (18) 
+297°/49 X 1350+7/49 X675). (19) 
For the respective molecules F2, Cle, Iz, one finds 
y¥=2.56r/ao, 2.03r/ao, and 1.9r/ap. 


In order to complete the evaluation of Qwe 
using Eq. (12), we need the quantities S,,. and 
Sx’. One could of course neglect these without 
gross error, but they can be evaluated fairly 
easily, and will be needed later. Substituting 
Eqs. (12), (13) into Eqs. (8), transforming to 
The values of the parameters are those given by 


Slater, except for fluorine, where Zener’s Z*=5.12 (Phys. 
ees” 51 (1930)) has been used instead of Slater’s 


(6) 
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elliptic coordinates, and integrating, the results 


are as follows 


It is of interest that the expression for S,,” in 
Eqs. (20) is exactly the same as that for Q,e in 
Eq. (16) except for the omission of a factor 3r 
which in the latter betokens its dipole-moment 
character. The same relation between S,,' and 
Q,-' holds in general,!* so we need not give the 
Sy’ equations explicitly for Cle and Ie. 


Clo: Seo 
+34y4/1575 —2y3/75 
—9y?/25—y—-—1). (21) 
To: Seo +23y7/66150 
+67y°/66150+ 75/294 
— 1/735 —17y*/210 
—29y?/70—y—1). (22) 
III. CaLcuLaTion oF DIPOLE 
STRENGTHS: AO THEORY 


Before calculating numerical values of Qy@ and 
Dye, let us evaluate Qve in terms of the AO 
approximation. The antisymmetrical AO wave 
functions for states NV and pure !0 are as follows: 


¥(N) =A (23) 
Yo(20,) (24) 
where A and B are normalization factors, and 
vat = 
X myp(4) ya’ 


(25) 


13 For 2p AO’s, Soo and are given (under the re-. 


spective designations Soo? and S1:*) by J. H. Bartlett, Jr., 
Phys. Rev. 37, 512 and 522 (Table IX), 1931. 

4 This can be seen as follows: in Eq. (12), the part of the 
integral which comes from the z in the factor (z+ }r) of the 
integrand vanishes. The remainder of the integral, except 
for the factor 4r, is exactly S;x-. 
pit De method is that of J. C. Slater, Phys. Rev. 38, 1109 
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Yir=N 

X 

Viv = + 
X 

= Nyt 

X 

X 


The ... stands for (9) (10) 
and so on with suitable additional single-electron 
functions occupied by inner-shell electrons 
11---n. The symbols o, z,, zz stand here for 
npo, npr,, and npz,z of atom A, o’, m2’ for 
the same AO’s of atom B. The AO’s nupo and 
npo’ are to be taken with the same sign in the 
region where they overlap strongly. o(!Q.) is 
like Yo(!Q,) except that x’s and y’s replace each 
other throughout. 

The wave function given by Eq. (24) is one 
of four possible independent linear combinations 
of the four y’s of Eqs. (26). The four independent 
combinations, differing only in choice of signs 
(and slightly in the value of B) in Eq. (24), may 
be described and identified as follows. The sign 
combinations ++++, ++--, +-+-, 
+-—-—-+ correspond to the y component of the 
respective states *Qiu, (cf. Section I 
and Eqs. (1) in regard to the Q, states). 

Introducing ¥(N) and yYo(!Q) from Eqs. (23)- 
(26) into Eq. (2), and reducing the double p's 
to single Yp’s (thereby at the same time removing 
the (n!)-? factors) in the manner given by 
Slater,!® we obtain 


(27) 


(26) 


where, for example, 


f 


[= (1) (2)a(3) yp’ (4) 


X ++ (28) 
16 Reference 7, p. 1115. 
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To obtain Qyve, we have to evaluate the eight J 
integrals, also A and B. Roughly, A=2-! and 
B=i, and N,;=---Nyi=1. These relations 
would be exact if yx and yn in Eq. (23) were 
orthogonal, likewise in Eq. (24), and 
the various AO’s in Eqs. (25) and (26). The lack 
of orthogonality arises mainly’ from the fact 
that AO’s of two different atoms occur in each 
y of Eqs. (23)-(26), and any two such AO’s, if 
alike in \ and (if \>0) in x or y character, depart 
more or less from orthogonality. 

Now it can be shown that, if we neglect unim- 
portant factors, each of the J’s can be reduced 
to a four-electron integral,’* for example 


Jim = f 


x Typ (4) jr. (29) 


If the best Hartree AO’s are used, AO’s even of the 
same atom’ may not be orthogonal if alike in /, \, and x or v 
type; but for our purposes we can safely neglect this 
additional possibility. 

**The proof is as follows. Consider for example 
Ji, mr of Eq. (28). In the first place, all terms coming 
from yr°+*¥n are easily seen to vanish. Next, let us 
write instead of tya(3)ryp(4) 
Xtya’(5)ayg’(6) in the first part of Jz, m1. This does not 

alter the latter’s value, as can be seen by going back to yy 
in Eqs. (25) and first making the change there. We see that 
the terms from ys+. now vanish, leaving only those from 
yi to y4. Now consider the contributions to J of all permu- 
tations P of wave functions among electrons 5---n. It is 
found that the total effect of these is to introduce a factor 
_ which is a product Ix, :(1—Sx*) of all possible quantities 
(1—S;2), where Si:= that can be formed from 
pairs of wave functions k,/ selected from among those 
occupied by electrons 5- --m in the basic or identity permu- 
tation. Note also that electrons 5---n are assigned to the 
same set of closed shells in the part of the integrand pre- 
ceding as in that following the symbol P in Eq. (28). 
Now if we similarly consider the quantities Nj and Ny, 
each of them is found to contain a factor [1/Ix, 1(1—Sy22) }}. 

In the calculation of J, these factors just cancel the above 

factor I1(1—S,2). This suggests that we would obtain 

correct results by simply omitting all but four electrons and 

a set of four one-electron wave functions from the y’s of 

Eqs. (25)—(26). The wave functions to be dropped are those 

occupied by electrons 7---n in Eqs. (25)-(26), together 

with two out of those occupied by electrons 1 to 6; the 
choice of these two varies, however, from one J to another. 

The foregoing proof that J can be reduced to a four- 
electron integral chosen as above described is exact and 
rigorous except for one thing. Among the permutations in 

Eq. (28), we have not yet considered those which exchange 

the arguments of electrons 1 to 4 with those of 5---n. 

These have the effect of introducing a further factor 

Th, i(1— Sh, ?), where h corresponds to arguments of 

electrons in the o 1---4, / to arguments of electrons in 

the group 5---m. They also introduce some additive 

Contributions to J which are not expressible simply by a 

multiplying factor. It appears that these effects are 

ely but not exactly canceled out by similar effects in 
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From each n-electron integral Jr, s all those n—4 
one-electron wave functions whose electrons 
belong to the same closed shells'® both in Yr and in 
Ys are eliminated in forming the four-electron 
integral.'* In obtaining Eq. (29), for example, 
the basic permutation in yy is first altered by 
exchanging one-electron functions between elec- 
trons 3, 4 and 5, 6 (this does not alter the value 
of ¥1), whereupon it is seen that all the electrons 
5-++m are now in identical closed shells in yy; and 
Yim, while the electrons 1 to 4 are not in like 
closed shells in yy and Ym. To a sufficient ap- 
proximation, in Eq. (29); the »v’s 
have the same relation to 4-electron w’s as do 
the N’s to the n-electron y’s of Eqs. (25), (26). 

The four-electron J’s are now easily reduced 
to expressions in terms of one-electron integrals. 
Each J yields two nonvanishing contributions. 
When all these are put into Eq. (27), the result is 


where Q,, and Q,, are as in Eqs. (10), (11). One 


finds 


(31) 
1/B?=4(1 — + 


from which, to a close approximation, 


As was noted above (cf. Eq. (10)), Q..=0 if 
ordinary np AO’s are used. 


IV. CALCULATION OF N->!°Q DIPOLE 
STRENGTHS; COMPARISON WITH 
EXPERIMENT 


Using Eqs. (9), (32) together with Eqs. (10), 
(14), (16)—(22), we can now compute the quan- 
tities Qvg and Dye theoretically for the transi- 
tions N-1'Q in Fe, Cle, and Ix, for both the 
LCAO MO and the AO approximations. The 
values so computed for F2, Cle, and I, for their 
respective equilibrium internuclear separations, 
also interpolated values for Bre, are given in 


the N product. The conclusion then is that the four- 


electron J expression is a good approximation to the 
n-electron expression, but is not quite exact. (It appears 
that in special cases the error might be serious.) 

19 For this pur , a closed shell consists of two electrons 
with opposite spins occupying the same AO, e.g., npry. 
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Table I, together with corresponding experi- 
mental values taken from the analysis in another 
paper.* In the case of F2, Qve and Dye have 
been computed for three different values of r. 

Everything considered, the correspondence 
between observed and calculated quantities is 
extremely gratifying, if for the moment we omit 
iodine from consideration. In the other halogens, 
the AO calculated values are rather close to the 
observed. Notably, the AO calculated D values 
parallel very closely the twenty-fold increase 
in the observed D values from F: to Bre. The 
LCAO calculated values agree less well than the 
AO both in being several times too large, and in 
showing too great a ratio of intensity between 
F, and Br2; however, the progressive increase 
from F, to Bre is qualitatively correct. 

In calculations on N-V transitions, the 
scanty experimental intensity values seemed to 
fall about midway between LCAO and AO cal- 
culated values.? The present data and evidence 
on N-(@Q transitions, and further similar evidence 
to be presented in X of this series, tend to in- 
dicate that the AO approximation is the better.” 
This indication cannot be taken as proof, how- 
ever, because of the additional approximations 
’ involved in the use of Slater AO’s for N-@Q 
transitions, and not involved for N—V transi- 
tions, where the forms of the AO’s drop out.” 

The observed agreements for the N—-Q 
halogen spectra give confirmation to the existing 


interpretation of the observed bands and con- 


tinua in terms of electron configurations and 
electronic states (cf. Eqs. (1)). 

It remains to explain why the observed Dye 
for iodine is so.out of line with those for the 
other halogens, and with the calculated values. 
In this connection, it is probably significant that 
the Dyg values for HI and CH;I (to be discussed 
in X of this series) are likewise abnormally small 
compared with those of the other hydrogen and 
methyl halides. 

Empirically, there is a striking correlation in 
iodine of abnormally high intensity of N—-Qp, 
attributable only to strong case c effects, with 
the abnormally low intensity of N—*1Q,. There 
are also indications of a similar correlation in 


20 In this connection, the writer would like to emphasize 
again that the LCAO approximation, used here, represents 
only a crude approach to the best MO approximation. 
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CH,lI (see X of this series). It seems fairly safe, 
then, to attribute the low Qwe of iodine to some 
kind of case c disturbance. It might for instance 
be that one or both of the N and Q wave func- 
tions become so much admixed with others that 
considerable added contributions appear in the 
expressions for Qyve in Eqs. (9) and (32), and 
happen to be opposite in sign to Q,., thus partly 
canceling the latter. 


V. ConcLusIons. DISCUSSION OF 
APPROXIMATIONS 


The calculations for F2 in Table I illustrate 
that Dyg may vary very rapidly with the inter- 
nuclear distance r. It would be of interest to 
make similar calculations for the other halogens.* 
If the result for F is at all typical, existing deter- 
minations of upper potential energy curves from 
intensity distributions in continuous N—-@Q 
spectra (halogens, hydrogen and alkyl halides) 
should undergo reconsideration. Rapid variation 
of Dye with r would considerably alter ymax and 
the intensity distribution for a continuum cor- 
responding to a given upper and lower potential 
energy curve.” 

Referring to Eqs. (9) or (32) and (16), it will 
be seen that the variation of D with r comes in 
mainly through the parameter y=rZ*/n*ap. 
Hence Dyg should be about equally as sensitive 
to change in the assumed value of Z*/n* as to 
variation of ry. Thus for Fs, an increase of 20 
percent in Z* would cause the computed Dye in 
Table I to be decreased to about one-tenth its 
value according to the AO, or to about one- 
fourth its value according to the LCAO method. 
It is therefore of interest to consider whether 
our Z* values were well chosen. 

For the computations reported in Table I, the 
Z* values suggested by Slater for free atoms were 
used. It is, however, known that when the vari- 
ation method is used to minimize the energy of 
molecular wave functions, treating Z* as a 
parameter, the optimum value of the latter may 
be considerably different than for free atoms. 

* Note added in proof—Calculations being made by Dr. 


A. Sayvetz, show similar rapid variation for all the halo- 
gens. 

21H. M. James and A. S. Coolidge, Phys. Rev. 55, 184 
(1939) have also found a considerable variation of D with r 
in the case of the transition 2sa,, souls, of 
the hydrogen molecule. Their calculations were made with 
accurate wave functions. 
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The method has been applied principally to the 
normal states of H2*+ and of He, where (both in 
AO and in LCAO approximations), values near 
1.2 are found for Z*, as compared with 1.00 for 
the free atoms.” 

In the normal states of H,* and He, orbitals 
or interactions are purely bonding. In the halogen 
N and Q states, antibonding as well as bonding 
interactions are present. It therefore seemed de- 
sirable to investigate minimization of the energy 
by varying Z* for an antibonding state. 

The simplest example of the latter is the o,1s 
antibonding MO of H2*, given in LCAO ap- 
proximation by 


o,1s= (dis— dis) /(1 — Sis, 1s’). 


On minimizing the energy for r=2do (near r, of 
normal H¢*), the best Z* was found to be roughly 
0.84. This may be contrasted with the optimum 
Z* value 1.228 found by Finkelstein and 
Horowitz for the o,1s bonding MO of Het. The 
comparison indicates that the optimum Z* is 
increased for bonding but decreased for anti- 
bonding states.”*¢ 

2 Cf., e.g., L. Pauling and E. Bright Wilson, Jr., Intro- 
duction to Quantum Mechanics (McGraw-Hill, New York, 
1935), pp. 331, 349. 


%e Analogy to Hz and He, assumed here, may of course 
possibly not be a reliable guide for estimating in the case 


In the N—@Q halogen spectra in MO approxi- 
mation, an electron jumps from one antibonding 
MO (a,np) to another (co,p). Using, in the 
LCAO forms, Z* values less than for the free 
atoms would increase the computed Dyg values, 
thus making their agreement with experiment 
worse than before ;?* quite possibly not much 
worse, however, since it appears probable that 
the Z* changes may be relatively much smaller 
than in the case of He*. 

Turning now to the AO approximation (cf. 
Eqs. (1), (23)-(32)), a qualitative examination 
shows that a nearly nonbonding z AO of state NV 
and a bonded oa AO of state Q are involved in the 
integrals S,,, and Q,., which appear in the ex- 
pression for Qyg. Presumably, then, the op- 
timum Z* is about equal for z to the Slater value, 
but larger for ¢. The calculated Dye is then 
somewhat decreased. Thus the effects on Dwg of 
optimizing Z* are apparently opposite for the 
LCAO and AO approximations, giving improved 
agreement with experiment only in the latter 
case.?24 

A fault of the Slater AO’s, not removed by 


of halogen molecules whether the optimum Z* is greater or 
less than for free atoms. Also, the absence of radial nodes in 
the Slater functions may conceivably affect our arguments 
(except for fluorine). 


TaBLeE I, Calculated and observed dipole strengths for N—>Q in the halogens. 


MOLECULE PARAMETERS? INTEGRALS? DIPpoLe STRENGTH* 
LCAO Catc.5 AO Catc.é Oss.” 
re(A)! §r(A) See’ Sug’ |} Ong(A) | Dne 
F; 


0.64 6.21 0.235 0.077 0.0494 0.083 | 0.0069 0.033 0.0011 
1.45+.05| 0.70 6.8 0187 0.053 0.037 0.060 | 0.0036 0.020 0.00040; 0.018 0.00034 
0.75 7.28 0.151 0.038 0.0285 0.045 | 0.0020 0.012 0.00014 
Cl, 
1.983 0.99 7.64 0.302 | 0.100 0.0985 0.178 | 0.0317 0.084 0.0071 0.053 0.0028 
Bre 
0.077 0.0060 
or or 
2.284 1.14 8.88 [0.204 [0.0416]| [0.110] | [0.0121]} 0.084 0.0070 


2.667 1.33 9.57 0.392 | 0.091 | 0.121 


0.230 | 0.0529 | 0.136 | 0.0185 


(0.025) | (0.0006) 


I 1 The equilibrium internuclear distances re are from band spectrum data except for F2, where the value given was obtained by L. O. Brockway, 
‘a 


-Am. Chem. Soc. 60, 1348 (1938), and is rather uncertain. 
6.1, 7.6, and 7.6, n* =2, 3, 3.7, 4 for F, Cl, Br, and I, respectively. 


2 The parameters are: 7 is the internuclear distance, and y =rZ*/n 
3 Cf. Eqs. (16)-—(22). 


o, with Z* =5.12, 
‘The dipole strength Dyg is the square of Ove, ; 


the latter, given by Eq. (2) of the text, is conveniently measured in angstrom units (A). Qnq has not been calculated for Br2; the ‘‘calculate 
aes given in brackets for vg have been interpolated between those calculated for Cl2 and for Iz, and the bracketed Dwg values are the squares 


these. 5 Cf. Eqs. (9), (10) of the text. 
e's and N-—!Q given in the last column of Table I of reference 4 (cf 


6 Cf. Eqs. (32), (10) of text. 


7 The observed Dwg values are obtained by adding the D values for 


. text above, preceding Eq. (2)). For Br2, the two values given here corre- 


nd to the alternative interpretations I and II of reference 4. For Iz, the value is an average of those corresponding to interpretations I and II 


N 
pyre: ane 4, but is in any case based 


on very rough data, and may be in error by a factor as large as 3 or 5. 
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optimizing Z*, is that they correspond to an 
inverse square field. The best AO’s would cor- 
respond to a Z* decreasing from inside out, to an 
asymptotic value 1 at large r. Qualitative con- 
sideration indicates that such best AO’s of two 
atoms would overlap more, and thus increase 
Dye, as compared with those based on constant 
Z*. Hence it appears that this kind of improve- 
ment in the wave functions would make the 
calculated Dyg’s agree less well than before with 
experiment, for both AO and LCAO approxima- 
tions.?* 

On the basis of the foregoing discussion, it 
appears not very useful to try to improve the 
wave functions by improving the forms of the 
AO’s used in the AO or LCAO approximation. 
Fortuitously, perhaps, the simple unimproved 
Slater AO’s give relatively good results. It may 
seem anomalous that improved wave functions 
can give worse calculated Dy@g values,”* but this 
is entirely possible so long as the wave functions 
_ used are still at best only rough approximations. 
Apparently it happens that those parts of the 
wave functions which are of predominant im- 
portance for Dyg are more nearly correct using 
the original Slater AO’s than when each wave 
function as a whole is improved? by varying Z*. 
It should be borne in mind, however, that the 
observed Dyg values, since they are rather small, 
and thus relatively sensitive to perturbing 
effects, should be used with some caution as a 
test of the goodness of various approximations. 

For the sake of completeness, we may consider 
the improvement of the Slater AO’s in molecular 
wave functions by ‘“‘polarization.”” In the MO 
approximation, we may reasonably expect this 
to give better calculated Dyg values, since in the 
case of MO’s “polarization” means essentially 
getting away from the limitations of the very 
simple LCAO approximation, in the direction of 
using “true” or ‘best’? MO’s.?* Taking as an 
example the and MO’s of a halogen 
molecule, we may write expansions in pairs of 
AO’s as follows, writing the presumably most 

** Compare the work of H. M. James and A. S. Coolidge, 
J. Chem. Phys. 1, 825 (1933); 3, 129 (1935). Here, to 
sure, an expansion in terms of elliptical coordinates is used 
(also containing extra terms in 712 and 2 carrying the 


wave functions even beyond the best MO approximation), 
- > principle is the same as for an expansion in pairs 
s. 


MULLIKEN 


important terms first : 
o.np=a(npo —npo’)+B(ndo—ndo’) 
+y(nso—nso’)+---, (33) 
np, 


The first pair of terms in each of Eqs. (33) 
gives the LCAO approximation, further pairs 
may be thought of as corresponding to the use 
of polarized AO’s, thus: 


por), 


(34) 


where 
(npo poi) (35) 


is an example of a polarized AO. 

In Eq. (33), a, 8, y and so on can be treated as 
parameters in a variation problem. To minimize 
the energy of an antibonding o,np MO expressed 
in the form (33) or (34)—(35), overlapping of the 
two AO’s and po: in Eqs. (34) should 
be diminished. From the forms of npo and ndo 
one readily sees that for this purpose B/a in (33), 
or 1/k in (35), should have negative sign. Numer- 
ical values, and further coefficients, might be 
obtained by calculations. If we substitute 
polarized AO’s obtained in this manner for the 
ordinary AO’s in Eqs. (9)—(11), the latter ex- 
pressions for Qye@ are converted from rough 
LCAO-approximation-MO equations to best-MO 
equations. 

In Eq. (11), the integral Q,,, hitherto vanish- 
ing, now assumes importance. Considering just 
the effect of ndc, which except in F2 is probably 
the principal polarizing term in Eq. (35), we now 
have 


f 


with / negative (see above). Since 
JS (npry)y(ndo)dv 


is large and negative, it is seen that Q,, assumes 
a considerable positive value. At the same time, 
Qre of Eq. (10) is somewhat diminished because 
of the diminished overlapping of the polarized 
AO’s now involved. Thus without much doubt 
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polarization causes Q,, partly to cancel** 
and so to reduce the calculated MO approxima- 
tion Dye, making it agree better with the ob- 
served value.” This improvement most probably 
much outweighs the contrary effects” obtained 
by optimizing Z*. 

The effect on Dye of using polarized AO’s is 
quite different in the AO than in the MO ap- 
proximation. Referring to an earlier paragraph 
on the variation of Z* in the AO approximation, 


% Presumably Qyo does not become so large as to more 
than cancel Qze:, but this cannot be shown by the present 
qualitative arguments. 

.*8 This applies to the heavier halogens. For F;: the result 
is not qualitatively obvious, but presumably still holds, 
since the best-MO approximation pretty surely should give 
better Dwg values than the LCAO approximation. 


similar remarks now apply for improvement of 
the AO’s by polarization. The only AO’s for 
which improvement is important are the npo 
AO’s in state Q. These are occupied by electrons 
which are somewhat bonding, so to minimize the 
energy the polarization should be such as to give 
increased overlapping. This would be obtained by 
using Eq. (35) with 1/k positive, instead of 
negative as for best polarization of the LCAO 
MO’s. 

The effects of such polarization on Qy@g can 
be seen from Eq. (32). S,. and Q,. are in- 
creased, while Q,, takes on a negative value, so 
that Dyg is increased. Agreement with experi- 
ment is made worse, opposing*** the effect ob- 
tained by optimizing Z*. 
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The Hydrogen Isotope Exchange Reaction 


J. O. HaLForp 
Department of Chemistry, University of Michigan, Ann Arbor, Michigan 


(Received August 25, 1939) 


The possibility of deriving the equilibrium constants for the higher deuterium substitutions 
from the isotope exchange coefficient is examined with the aid of calculations on the ammonia 
and methane reactions with water. If the experiment is carried to solutions at equilibrium with 
water containing as much as 25 percent D.O and the error in a single measurement does not 
exceed 1.0 percent the second constant is determined to within 10 percent of the correct value. 
As a standard of comparison against which the experimental results may be examined to ad- 
vantage, an equation is developed showing the relation between the equilibrium constant for 
the first substitution and the exchange factor at higher deuterium fractions, subject to the 
assumption that each successive substitution of D for H produces the same frequency shift. The 
evaluation of frequency shifts from exchange data is considered. 


HE study of hydrogen isotope exchange 
equilibria has been restricted almost exclu- 
sively to systems of low deuterium fraction in 
which only a negligible proportion of the mole- 
cules can acquire more than one deuterium atom. 
For exchange reactions involving water this 
restriction limits the interpretation to the first 
of the reactions 


AH,+pDOH=AH,_,D,+pHOH. (1) 


In principle the measurement and interpretation 
need not be so limited. It is possible to perform 
exchange experiments over a wide range of 
isotope fractions involving systems in which the 
higher deuterium substitution products are pres- 


ent at considerable concentrations. From such 
experiments, if they were sufficiently accurate, 
all the equilibrium constants of the reactions 
expressed by Eq. (1) could be derived. The 
success of an actual interpretation, however, is 
limited by the accuracy of the individual dis- 
tribution measurement and by the fact that the 
relative scarcity of deuterium makes it imprac- 
tical to obtain data over the entire range of 
isotope compositions. 

In the following paragraphs it is shown by 
calculations with the ammonia and methane 
reactions that the constant for the second sub- 
stitution can be obtained with a fair degree of 
accuracy even with a liberal allowance for experi- 
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mental error and that the satisfactory determina- 
tion of the higher constants would involve a 
greater precision than is normally achieved. 
The equilibrium constants are derived from the 
distribution by trial. The extrapolation to zero 
deuterium fraction for K, and the selection of 
first approximations to Ke and K; are facilitated 
by the use of some limiting algebraic relations. 
For their development it is convenient to define 
Ba by the equation 
p=n p=n 
pA,/ LA), (2) 
p=1 p=0 
in which A, represents the number of molecules 
of AH,_»D, in the equilibrium mixture. Eq. (2) 
then describes the exchange experiment in terms 
of the average number of D atoms per molecule 
of AH, and its deuterium substitution products. 
Each substance AH,_,D, will come to equi- 
librium with HzO, DOH and AH, in accordance 
with reaction 1 for which the mass law ex- 
pression is 
A,=AoK>R?. (3) 
The concentration ratio DOH/HOH has been 
represented by R. Eqs. (2) and (3) are combined 
to yield the relation 


p=n , p=n 
Ba= bKpR?/1+ K,R’. (4) 
p= p= 
The corresponding expression for BH20 takes the 


form 
Bu.o=R(1+2KR)/(1+R+KR?), (5) 


with K denoting the equilibrium constant of the 
reaction 


2DOH = H.0+D.,0. (6) 


In order to develop Eq. (4) any further it is 
necessary to assume relations between the n 
constants K,. The numerical or symmetry effect 
can be separated from the other statistical factors 


by means of the equation 
kp n! 


K,=————_. (7) 
2” (n—p)!p! 
If k, is set equal to k,” it follows that 
Ky P n! 
K,=(— (8) 
(n—p)!p! 


Eq. (8) describes the limiting case for which the 
same factors of mass, moment of inertia and 


frequency could be used to calculate the equi- 
librium constant of each successive reaction in 


the series 
AH,-»,D,+ DOH =AH,_»-1D p41 +H20. 


This condition is fulfilled closely enough by the 
mass and inertia factors but may be far from 
true for the vibrational effect. Consequently the 
derivation leads to a relation in which Ky,’ varies 
with B64 in a manner determined primarily by the 
change in frequency shift with successive sub- 
stitution of D for H. 

Equation (8) is equivalent to a conclusion 
reached by Brodskii.' It is combined with Eq. (4) 
and simplified as follows: 


=(—) n! 
n (n—p)'(p—1)! nK,'R 


Ba= = 
p=n nn! 
1+ (—) 
n (n—p)!p! 


It follows that 
Ky’ =nB4/R(n—Ba) (10) 


for the limiting case under consideration. 

The numerical factor introduced into Eq. (7) 
requires no justification, since it merely sets the 
probability of occurrence of AH,_»D, propor- 
tional to the number of permutations of (n—p) 
atoms of the one isotope with p atoms of the 
other. It is of interest, however, to consider the 
relation of this factor to the symmetry number 
which occurs in the simplified form of the 
partition function or sum over states. 

In the general case, AH,_,D, may represent 
several molecules, the number depending on the 
symmetry number oy of AH,. To the ith molecule 
is assigned the symmetry number ¢a,;. The con- 
stant K,; of the reaction 


AH,+pDOH (11) 
then becomes 
K (Rpi/2”) (00/0 i). (12) 


The relations equivalent to Eqs. (7) and (8) are 
obtained by summing over the 7 molecules: 


1 kp 
i 2” Opi 2? Opi 


1 Brodskii, Trans. Faraday Soc. 33, 1180 (1937). 
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methane-water exchange reactions at 25°C. 


(“)' 
Mo —Av Kp (n—p)!p! 


H.0 1.000 1.000 0 _ — 
DOH | 1.085 1.659 11334 | — ee 
D.0 1.172 2.632 23774 | — ome 


NH; 1.000 1.000 0 

NH:D | 1.090 1.457 1097* | 1.214] 1.214 
NHDz | 1.183 2.043 2417* |0.804} 0.491 
ND; 1.278 2.752 3850* | 0.223] 0.066 


CH, 1.000 1.000 0 

CH;D | 1.096 1.342 1161¢ | 1.747] 1.747 
CH2D, | 1.194 1.776 2441f 1.145 
CHD; | 1.296 2.269 3725¢ |0.553| 0.333 
CD, 1.400 2.827 5004f | 0.073} 0.036 


4 Barker and Sleator, J. Chem. a. 5 660 (1935). 
* Howard, J. Chem. Phys. 3, 208 (193. 
ad t _— Morikawa, Barnes and Taylor, J. Chem. Phys. 5, 1 


TABLE II. Ammonia exchange reaction at 25°C. Distribution 
from molecular constants. 


R %D20 B Ky 
0.6 24.09 0.7027 1.529, 
4 17.22 4825 1.437; 
3 13.40 3653 1.3865 
9.27 1.3323 
4.81 .1223 1.274, 
05 2.45 .0610 1.244, 
0 0 0 1.214 


The factor kp is a weighted average of the 
values of k,:. If both Eqs. (7) and (13) are 
correct it is necessary to set }°j;00/o»; equal to 
n!/(n—p)!p!. 

In the limiting case for which oo=1 each op; 
must be equal to unity because each substitution 
of D for H will result in a recognizably different 
molecule. The ratio o0/¢; is equal to unity and 
the sum over the 7 forms will be equal to the 
number of forms n!/(n—p)!p!, that is, to the 
number of distinguishable permutations of pD 
atoms with (n—p)H atoms. If, however, oo 
equals abc, the product of a set of numbers repre- 
senting equivalent arrangements produced by 
internal rotation or rotation of the entire 
molecule, the substitution of pD atoms will 
eliminate one or more of these factors, so that 
7p: might be equal to a and ao/e»: equal to be. 
The rotational operations corresponding to 6 and 
c, when performed upon the substitution product, 


TaBLE I. Statistical constants for ammonia-water and 


are no longer symmetry operations, but instead 
show the equivalence of bc of the total possible 
permutations. In this way they give to the form 
of symmetry number @ a statistical weight equal 
to the number of permutations which can be 
shown to be equivalent by performing just those 
rotations which were eliminated as symmetry 
operations by the substitution. None of the 
n'/(n—p)!p! permutations are lost and each one 
must belong to one of the 7 forms. The quantity 
must therefore be independent of the 
actual symmetry numbers and equal to the 
number of permutations. It can be shown, in 
addition, that this relation is a necessary conse- 
quence of the accepted principle that the 
difference between isotopes of slightly different 
mass, such as Cl® and Cl*’, may be taken into 
account or neglected at will without introducing 
significant statistical errors. 

For the purpose of illustrating the sensitivity 
of Ke and K; to variations in the distribution 
values from which they have been calculated, 
it has seemed preferable to use the ammonia 
and methane reactions, rather than to set up a 
purely hypothetical case. All molecular data 
necessary to a complete description of these 
systems are available. Table I summarizes the 
values of K, at 25°C together with the mass 
ratios, moment ratios and frequencies used in 
their calculation. K, has been obtained from the 
relation K,=L,/Lpon” where 


A,B,C, n! 
Mo AoBoCo (n—p) \p! 


Factors of the form (1—e~**”/*7), associated with 
individual frequencies, have been omitted from 
Eq. (14) because their contribution to K, is 
negligible. The moment ratios are based upon 
the assumption that interatomic distances and 
bond angles do not change with the isotope sub- 
stitution. The tetrahedral angle has been used 
for methane, and for water and ammonia, re- 
spectively, the angles taken are 105° and 113°. 
It is possible that some of the frequency shifts 
may be in error, since they are based in part 
upon calculated quantities. For example, the data 
of Ginsburg and Barker? set the first methane 
shift at 1088, which would give a constant about 


? Ginsburg and Barker, J. Chem. Phys. 3, 668 (1935). 
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TABLE III. Methane exchange reaction at 25°C. Distribution 
from molecular constants. 


R %D:0 B Ky 
0.6 24.09 0.928, 2.0145 
A 17.22 .6535 1.952. 
3 13.40 5:02 1.913; 
9.27 1.8666 
Al 4.81 1.8117 
05 2.45 .0871 1.780; 
0 0 0 1.747 


15 percent lower than the value of 1.747 which 
appears in the table. This difference is not 
important for the present purpose. 

The constants of Table I are used to yield 
values of 8, in accordance with Eq. (5), for 
selected values of R. These quantities, together 
with K, from Eq. (10) and, for reference, the 
deuterium oxide concentration of the equilibrium 
water in the conventional units, are presented in 
Tables II and III, respectively, for the ammonia 
and methane reactions. 

The practical experimental range has been 
taken to involve equilibration with water con- 
taining not more than 25 percent D,O, although 
for highest accuracy in determining the constants 
K, the distribution should cover the entire range 
of isotope fractions, and each constant should be 
derived mainly from the distribution in the 
region in which AH,_,D, contains its maximum 
fraction of the deuterium. Tables II and III 
show the character to be expected for typical 
isotope exchange data involving small molecules 
in which the replaceable hydrogens are linked to 
a single atom. Such data could be conveniently 
extrapolated to zero deuterium fraction to yield 
the limiting value of K, by plotting both K,’ and 
8/R against R and taking the best intercept 
common to the two functions. For selecting a 
preliminary value of K2 from experimental data, 
a comparison of K,’ with the quantities appear- 
ing in Tables I, II and III should prove to be 
useful. In this connection, it is observed that the 
deviation of K,’ from K; at a given value of R is 
about the same for the ammonia and methane 
reactions, in spite of the fact that Ke is much 
closer to K,?/n in the latter case. This indicates 
that the deviation of K, from constancy is larger 
for larger values of m, other factors being 
constant. 

In testing the sensitivity of Ke to uncertainties 
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in the data from which it is derived, it is assumed 
that each experimental value of 8 will be accurate 
to 1.0 percent, and that 6 values used in the 
interpretation will be taken from a smooth 
curve of B against R, so that the sum of the 
errors in two 8 values separated by 0.3 unit on 
the R scale will not exceed 1.0 percent. In Tables 
IV and V are shown Kz and K; as derived from 
K, and values of 8 at R=0.3 and R=0.6 with 
these quantities differing from those of Tables II 
and III in the manner indicated in the first 
three columns. The fourth and sixth columns 
show Ke and K3, the fifth and seventh the 
corresponding frequency shifts when the entire 
deviation from Eq. (8) is assumed to be due to 
the frequencies. This assumption introduces an 
error of 18 cm in Ave for ammonia and would 
have a much smaller effect in the interpretation 
of data for larger molecules. 

The largest tabulated deviation of Kz from the 
correct value is 8.2 percent. The corresponding 
error in Ave is 53 cm—, of which 18 cm could be 
accounted for in a more detailed calculation. 
It is interesting to observe that this is consider- 
ably smaller than the difference of 73 cm 
between the calculated and observed values of 
Av; for methane to which reference was made in 
connection with Table I. It is evident that the 
determination of Ky and Ave from exchange 
equilibrium data of reasonable accuracy is 
feasible, but that satisfactory evaluation of the 
corresponding quantities for the third substitu- 
tion would require measurements of higher 
precision extending over a wider range of 
isotope fractions. Exchange data can thus pro- 
vide an additional condition which must be met 
by an assignment of frequencies to a series of 
molecules >. 

There are no recorded exchange experiments 


TABLE IV. Ammonia exchange reaction at 25°C; Variation 
of and K; with K, and 


1.0% VARIATION 


Ki} (0.3) B(0.6) Ke Av2(2417) Ks Avs(3850) 


0.223| 3796 
-738 2364 -285| 3899 
.763 2377 3815 
855 2424 150} 3632 
805} 2400 3853 


| 


re 


Id 
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with any substance for which >1 which extend 
over a sufficient range of isotope fractions to 
permit an accurate interpretation along the lines 
illustrated above. Some results’ obtained in this 
laboratory with acetone, however, may be ex- 
amined with profit in spite of their inadequacy. 
They are shown in Table VI. 

From the last four listed measurements it 
appears that the limiting K, at 25°C should be 
about 2.45, giving 0.82 for the exchange factor. 
This is in very close agreement with the reciprocal 
of Brodskii’s exchange factor of 1.22 and agrees 
well with the calculated result of K. Wirtz.‘ 
A factor for internal rotation or torsional vibra- 
tion of the methyl group, equal in the first 
approximation to (4/3)!, is introduced into Eq. 
(14), and Ay; is calculated to be —1217 cm. 
This method of taking internal rotation into 
account, which has been used by Hobden, 
Johnston, Weldon and Wilson,’ by Halford and 


TABLE V. Methane exchange reaction at 25°C; Variation 
of Kz and K; with K, and B. 


1.0% VARIATION 


Ki} B(0.3) B(0.6) K2 Av2(2441) Ks Av3(3725) 


1.495 2433 0.553 3694 
1.605 2461 458} 3615 
1.555 2450 558 3698 
1,412 2409 680; 3780 
1.514 2438 479 3634 


oo++o 


* Halford, Anderson, Bates and Swisher, J. Am. Chem. 
Soc. 57, 1663 (1935). 

‘ Wirtz, Zeits. f. physik. Chemie, B31, 309 (1936). 
61 (IS) Johnston, Weldon and Wilson, J. Chem. Soc. 
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’ Halford and Pecherer, J. Chem. Phys. 6, 571 (1938) 


TABLE VI. Acetone exchange reaction. 


Temp. (°C) %D:20 R B Ky 
80 5.17 0.1080 0.2789 2.71 

80 4.14 0857 2153 2.61 

80 1.89 0384 .0860 2.27 
80* 0.919 01852 0400 2.17 
60 1.136 02293 0524 2.31 
40, 35 1.167 02356 2.40 


* Reversed reaction, subject to larger experimental error. 


Pecherer,® and by Wirtz, is equivalent to an 
assumption that the contribution of this type of 
motion to the exchange equilibrium constant is 
not sensitive to the height of the potential barrier. 
The frequency shift determines the temperature 
dependence of the equilibrium constant and per- 
mits the calculation of K, at 80°C to be 2.37, a 
value considerably higher than the experimental 
one at this temperature. The calculated value is 
taken to be the more reliable and is combined 
with the first two measurements of Table VI to 
yield K2=4.1 and Ave= —2650 cm-. These are 
minimum magnitudes consistent with available 
information. 

With this interpretation the acetone reaction 
falls at least qualitatively into line with other 
exchange reactions. The frequency shift pro- 
duced by single substitution in the methyl group 
is numerically greater than Av; for methane by an 
amount comparable to the parallel difference 
observed between the alcohols and water. In 
addition, the relation between the indicated 
values of Av; and Av: for acetone is similar to that 
observed for water, ammonia and methane. 
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(Received December 13, 1939) 


The rate of conversion of parahydrogen and of ortho- 
deuterium is measured on paramagnetic crystals, viz., 
CuSQ,, NiCl, and neodymium oxalate. In the last case, 
crystals of approximately uniform size were used so that 
the total converting area was known. The conversion rate 
has a maximum at 0°C and it is shown that below this 
temperature the desorption rate of the Hz and the D» 


is rate-determining. Hence, the ratio of the rates for Hz and 
D: is not far from unity. The rate of desorption increases 
with increasing temperature and at higher temperatures, 
the rate of conversion on the surface becomes the slow 
process. The ratio of the rates for Hz and D2 becomes nearly 
11 in this region. The heat of adsorption of the hydrogen 
on neodymium oxalate is about 2.4 kcal./mole. 


HE collision of hydrogen or deuterium mole- 
cules with paramagnetic ions or molecules 
causes with a certain small probability’-> an 
ortho- paraconversion in the former. This con- 
version is induced by the inhomogeneous mag- 
. netic field of the paramagnetic particle.® 
The ortho- paraconversion on the surface of 
certain catalysts is also of paramagnetic nature. 
This is indicated by the experiments of 
Bonhoeffer, Farkas and Rummel’ which have 
made it very probable that in the case of char- 
coal, unsaturated, paramagnetic carbon atoms 
are responsible for the conversion. The conversion 
takes place in this case “inside the adsorption 
layer,” with a monomolecular mechanism. Taylor 
and Diamond® have demonstrated that the ortho- 
paraconversion is much more rapid on the surface 
of paramagnetic crystals than it is on diamagnetic 
materials. The present experiments deal with the 
ortho- paraconversion on paramagnetic crystals 
in a temperature region in which the amount of 
hydrogen adsorbed on the crystal is negligible 
and the conversion takes place during the 
collision process of the hydrogen molecule with 
the surface. 


1L. Farkas and H. Sachsse, Zeits. f. physik. Chemie 
B23, 1 (1933). 

?L. Farkas and H. Sachsse, Zeits. f. physik. Chemie 
B23, 19 (1933). - 

3 A. Farkas, L. Farkas and P. Harteck, Proc. Roy. Soc. 
A144, 481 (1934). 

a . 35) Farkas and L. Farkas, Proc. Roy. Soc. A152, 152 
a 938) Farkas and U, Garbatski, J. Chem. Phys. 6, 260 
*E. Wigner, Zeits. f. physik. Chemie B23, 28 (1933). 

7 Bonhoeffer, A. Farkas and Rummel, Zeits. f. physik. 
Chemie B21, 225 (1933). 
ass and Diamond, J. Am. Chem. Soc. 57, 1251 


There are two different types of collision 
processes during which an ortho- paraconversion 
may take place. (1) Simple collisions for which 
the hydrogen (or deuterium) molecule remains 
only for a short time in interaction with the 
crystal. The order of magnitude of this time is 
a/3c~10-* sec. denotes the diameter of 
the hydrogen molecule, ¢ its mean velocity). 
(2) Adsorptive collisions, when the mean time of 
sojourn of the hydrogen on the surface is still 
short compared with the time needed for estab- 
lishing the ortho- para-equilibrium but long com- 
pared with 10~" sec. 

In the simple collision processes we should 
expect an ortho- paratransition with approxi- 
mately the same probability as it occurs in the 
gaseous state. This is about 5-10" per collision. 
If, however, the time of interaction of the hydro- 
gen with the crystals is longer than 10~" sec. the 
transition probability during each of such a 
process may become larger. This will be discussed 
further below. The approximate amount of 
crystal surface required to obtain a measurable 
ortho- paraconversion by simple collision proc- 
esses is easily estimated. If N is the number of 
hydrogen molecules per cc, they strike Nc/4 
times a surface of 1 cm?. The half-lifetime of the 
conversion will be given by the expression 


41n 2 Nv 
(1) 
yNc s 


if y is the transition probability, v the total volume 
of the vessel, and s the total surface. Inserting the 
numerical values we obtain 7r=3.3-107v/s sec. 
Hence in a vessel of 20 cc and with a surface of 
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110 min. The surface of 105 cm? could be realized, 
e.g., by 1.7 cc substance consisting of little cubes 
of 10-” cm? (= 1y’) volume. 


EXPERIMENTS WITH PARAHYDROGEN ON CuSO, 
AND NiCl, (ANHYDROUS) 


In a vessel of 27 cc 65 percent parahydrogen 
was brought into contact with 3 g of CuSO, 
powder at 25°C. The copper sulphate was 
dehydrated before the experiment by heating it 
in vacuum. It disintegrated into a fine powder 
during this process. Parahydrogen was then 
introduced into the vessel and the rate of con- 
version followed by the microconductivity 
method.® The conversion proceeded according to 
the equation 

**, (2) 


where u; and wo denote the concentration of 
parahydrogen above the equilibrium concen- 
tration at the times ¢ and t=0. The half-time of 
conversion was 45 minutes, so that k=In 2/45 
min.—'. Using the same crystal powder, Eq. (2) 
was well reproducible on successive days. No 
special precautions were needed to prevent 
poisoning of the surface of the catalyst and it 
could be left in the vacuum without making 
efforts to exclude tap grease, etc. 

The conversion of p-Hz on the paramagnetic 
CuSO, is exactly of first order—the lifetime 
being independent of the pressure in the range 
(0.2 mm-—60 mm) investigated. On the other hand, 
the rate of conversion is increased proportionally 
with the amount of CuSO, used. The rate was the 
same no matter whether the CuSO, was in one 
heap at the bottom of a thick tube, or whether it 
was spread over a larger surface. This shows that 
the hydrogen penetrates rapidly into the space 
between the particles of the powder and the rate 
of this diffusion is fast compared with the rate of 
conversion. 

Using NiCl, as catalyst at 25°C similar results 
were obtained. However, the diamagnetic ZnSO, 
powder, dehydrated in the same way, failed to 
induce any measurable conversion within 6 days 
even when used in ten times greater amounts 
than the CuSO, was used. 


* A. Farkas, Zeits. f. physik. Chemie B22, 344 (1933); 
A. Farkas and L. Farkas, Proc. Roy. Soc. Al44, 467 (1934). 
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105 cm? the half-lifetime of the conversion will be 
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CONVERSION OF ORTHODEUTERIUM ON CuSO, 
AND NiCl. 


The ortho- paraconversion of deuterium, in- 
vestigated in a similar way, proved to be about 
ten times slower than that of hydrogen. The same 
ratio for the rate constants has been observed in 
presence of paramagnetic gases. The slower 
rate of conversion of deuterium is due to the 
smaller magnetic moment of the deuteron, the 
rate of conversion being, ceteris paribus, pro- 
portional to the square of the nuclear magnetic 
moment of the material to be converted.® 

These experiments prove that in the present 
case the kinetics of the conversion is analogous 
to that in the gaseous state, the absolute rate 
being proportional to the number of molecules 
striking the paramagnetic surface, each collision 
inducing with a small probability an ortho- 
paratransition. 

In the heterogeneous catalyses hitherto in- 
vestigated the conversion took place at low 
temperatures in the van der Waals adsorption 
layer and the rate of conversion was determined 
by the desorption velocity. Consequently, as soon 
as the adsorption layer became saturated with 
hydrogen, the absolute rate of conversion became 
independent of pressure. The reason for this is 
that the time of sojourn of the Hz or Dz on the 
surface is so large that the ortho- para-equilibrium 
is practically established for the molecules leaving 
the surface. This is the situation already at 
moderately low pressures and since the rates of 
desorption for light and heavy hydrogen are not 
very different, a low ratio of conversion rates was 
observed." 


EsTIMATE OF THE COLLISION EFFICIENCY 


In order to gain some information on the transi- 
tion probability of the ortho- paraconversion 
during the collision process with the paramagnetic 
solid, two series of experiments were carried out. 
(1) Estimation of the collision efficiency at ordi- 
nary temperature. (2) Investigation of the tem- 
perature dependence of the conversion rate. 

In order to determine the collision efficiency of 

10 Kalckar and Teller, Proc. Roy. Soc. A150, 520 (1935). 
11 The difference of the desorption rates is caused by the 
difference in zero point energies of Hz and Ds, that of the 


latter being somewhat smaller and the desorption rate 
correspondingly lower. 
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the conversion it is essential to estimate the total 
surface of the paramagnetic solid. 

Neodymium oxalate-hydrate crystals were 
chosen for this purpose as they form very small 
crystals of fairly well defined shape, when pre- 
cipitated from nitrate solution. As they are not 
hygroscopic, they can conveniently be handled 
and also examined under the microscope. They 
were obtained by precipitating neodymium 
nitrate with oxalic acid at 50°C in dilute (0.05 
molar) solution. The half-time of the parahydro- 
gen conversion with 1 g of substance in a reaction 
vessel of 27 cc volume was 7.3 minutes. Exami- 
nation of the crystals under the microscope 
showed that they formed prisms of a length 
between 0.5 and 50u. As there was a possibility 
that still smaller crystals were present, which 
could not be seen under microscope, the substance 
was sedimented several times in water, rejecting 
alternately the upper and lower fractions. The 
remaining substance—approximately ¢ of the 
original amount—gave a half-time of conversion 
of 5.4 min./g in the vessel of 27 cc volume. The 
crystal size was found to be much more uniform, 
the larger ones missing entirely. The average 
length of the prisms was 1y. Since the rate of 
conversion was not much changed by the process 
of sedimentation, it appears that submicroscopic 
crystals were not present in appreciable amounts 
in the original sample. 

According to Eq. (1), 1 g of neodymium 
oxalate (0.57 cc), consisting of particles of 1° 
size, should give in a vessel of 27 cc volume a 
half-time of conversion of 260’. The actually 
observed conversion was about 45 times faster. 
This discrepancy can be established, independ- 
ently of Eq. (1), by comparing the conversion 
caused by Nd*** ions in water solution, with the 
conversion on neodymium oxalate crystals. As- 
suming that the concentration of the Nd ion on 
the surface of the crystal powder corresponds to 
its concentration in the bulk of the solid, 1 g of 
our powder has 2.5-10!* ions on its surface. In a 
vessel of 27 cc this corresponds to a concentration 
of 1.8-10-* mole/liter which in aqueous solution 
would induce a paraconversion with the half- 
time? of 160’. Thus the rate is about 30 times 
slower than found on the surface of the powder. 
In fact the half-lifetime of the conversion on the 
solid should be even slower than 160’ if we take 


L: FAR'KAS AND L. SANDLER 


into account that the paramagnetic ions on the 
crystal can be approached from one side only. 
From these comparisons it is evident that the 
rate of conversion on the paramagnetic solid is 
about 30-50 times quicker than we would expect 
if it proceeded by simple collisions, as described 
in the first section. The simplest explanation of 
this result is that the time of sojourn of a 
hydrogen molecule on the surface is longer than 
the time of interaction in simple collision processes 
which we estimated to be 10~" sec. 

In order to understand the implications of this 
finding, we must consider the probability of 
conversion as a function of the time of sojourn 3 
on the surface. According to a formula previously 
derived,® the probability of conversion is pro- 
portional to the square of #, if # is very small 
compared with the rotational period of the 
molecule, ?<2J/h, where J is the moment of 
inertia of the molecule. If 3 approaches this 
quantity, the probability of conversion increases 
less rapidly and should, according to these 
formulae, cease to increase after # has reached 
the value of about 2//h. Actually the approxi- 
mations used in reference 6 break down if 
3>2J/h and a closer investigation shows that 
the conversion probability continues to increase, 
for 3>2J/h, proportionally with the first power 
of 3. Since even in the liquid or in the gas phase 
the time of collision is not far from 2J/h, we can 
assume that the roughly 50 times greater con- 
version probability of the collision with the 
crystal surface is caused by the 50 times longer 
duration of this collision, i.e., a 50 times longer 
time of sojourn of the Hz or Dz atom on the 
crystal surface than would correspond to a 
“simple collision.”” This indicates that the con- 
centration of Hz molecules on the crystal surface 
is about 50 times higher than in the gas phase, 
corresponding to a heat of adsorption Q of 


e@/kT — 50, 


i.e., Q=2.4 kcal. 


(3) 


TEMPERATURE DEPENDENCE OF 
THE CONVERSION 


In Fig. 1 the values of the half-lifetime of the 
conversion at 10 mm pressure are plotted as a 
function of the temperature. They represent the 
average values obtained in various runs with 1 g 


min. 
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Fic. 1. Conversion rate of parahydrogen as function 
of temperature. 


L 
-80° 


neodymium oxalate hydrate of 1p grain size in 
our vessel of 27 cc volume. The curve in Fig. 2 
represents” the ratio of the half-lifetime obtained 
with orthodeuterium to that obtained with 
parahydrogen. The curve in Fig. 1 has a mini- 
mum at 0°C (fastest conversion). At this point 
the ratio of conversion rates of light and heavy 
hydrogen, which is at higher temperatures con- 
stant (about 11) begins to decrease with de- 
creasing temperature. It is obvious that on the 
two sides of the minimum point two different 
processes are rate determining for the con- 
version. At the low temperature side the rate of 
conversion is determined by the rate of desorption, 
the number of hydrogen or deuterium molecules 
which are given off by the surface in the para and 
ortho states corresponding to thermal equilib- 
rium. In this temperature region, the conversion 
rates for light and heavy hydrogen differ only 
because of the somewhat lower desorption rate of 
deuterium. On the other hand, at high tempera- 
tures the mean time of adsorption is becoming 
shorter and the simple collisions begin to pre- 
dominate. This causes the decrease of the rate of 
conversion. The ratio for the conversion rates of 
light and heavy hydrogen evidently approaches 
the value 11 in this region, corresponding to the 
ratio of conversion probabilities in simple 
collisions. 

In the high temperature region, not the rate of 
desorption, but the rate of conversion on the 
surface is rate determining. The probability that 
a molecule is converted while on the surface is 


”% More recent experiments indicate that the ratio of 
the two conversion rates is somewhat higher at high 
temperatures than given in Fig. 2. It is 10.5 at 0°C; 
12.5 at 25°C and 13.5 at 50°C. 
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-80° 
Fic. 2. Temperature dependence of the ratio of conversion 

rates of parahydrogen to orthodeuterium. 

proportional to its time of sojourn &. The total 
amount of conversion per cm? surface is, there- 
fore, proportional to #Nc/4, where Nc/4 is the 
number of molecules striking 1 cm? of the 
surface per second. On the other hand, the total 
amount of hydrogen, adsorbed on 1 cm? of the 
surface is also 3Nc/4 so that the rate of conver- 
sion in the high temperature region is proportional 
to the adsorbed amount. This latter quantity 
again is proportional to exp (—Q/RT) where Q 
is the heat of adsorption. 

Figure 1 shows that the rate of conversion in 
the high temperature region decreases by a factor 
of about 2 if the temperature is raised by 50°. 
This gives the equation 


exp [Q(R/T —R/(T+50))]=2 


or 
Q=R-In 2(0.00057) ~2.4 kcal. (4) 


in surprisingly good agreement with (3). 

It is not certain up to what temperature the 
above described decrease of the conversion rate 
will continue, since at high temperatures other 
types of conversion, due to the activated 
adsorption, may become predominant. 

The present experiments are promising for 
measuring very small amounts of adsorption, the 
measurement of which has caused great diffi- 
culties thus far. At 10 mm pressure, the adsorbed 
layer is 0.001 molecular. The work here outlined 
allows also a measurement of both the amount of 
surface which is blocked for Hz by other adsorbed 
gases, and also the amount of surface left open. 
Most previous methods permit only the determi- 
nation of the area covered by foreign gases, not of 
the area left open. 
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IX. Further Studies of the Vapors of the Fatty Acid Series 
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The monomer spectra of the vapors of C;H;COOH, 
C;H;COOD and C,H;COOH have been obtained in the 
near infra-red with a rocksalt spectrometer. The association 
of each of the following molecules: HCOOH, HCOOD and 
CsH;COOD has been studied in the vapor phase by 
measuring the variation with temperature and pressure of 
the intensity of some association band. The bond strengths 
obtained for these fatty acids are 6.2, 6.4 and 6.9 kcal. per 
mole per bond, respectively. The fact that the boad 
strengths for light and heavy formic acids agree within the 
experimental error although the determinations were 


carried out with different bands is further evidence that the 
spectroscopic method is satisfactory. The lower value 
obtained for the bond strength in formic acid by this 
method compared to the vapor density value is discussed. 
In addition, it is suggested from the experimental data 
available at the present time, that it does not seem possible 
to claim the existence of a simple linear relationship 
between the bond strength and Av/v. Some results con- 
cerning »v(CH) absorption in the first five members of the 
fatty acid series are given. 


OME time ago Badger and Bauer! suggested 
that the frequency shift which a band 
v(X—H) experiences when hydrogen bonding 
takes place might be a semi-quantitative measure 
of the energy of bonding. They plotted Av/y 
corresponding to the frequency shift of a band 
due to hydrogen bonding against the bond 
strength in the hope that some approximate 
proportionality existed between these two quanti- 
ties and found that a straight line roughly fitted 
the data known at that time. Any such simple 
relationship would be of great importance be- 
cause it would lead to a knowledge of the bond 
strength from the relatively simple determination 
of the frequency shift which occurs when hydro- 
gen bonding takes place. However, the existence 
of such a relationship must be checked by the 
addition of reliable points to their plot. It is 
desirable to know the bond strength more accu- 
rately than is possible by estimations from latent 
heats. It would also be preferable to obtain the 
frequency shift by studying the vapor where 
possible so that the question of shifts due to the 
interactions present in the liquid phase and at 
times in solution would be avoided. From the 
results reported in this paper, as well as from 
earlier work?‘ it is possible to attempt the check 
an M. Badger and S. H. Bauer, J. Chem. Phys. 5, 839 
é' 3 ws. Bonner and R. Hofstadter, J. Chem. Phys. 6, 531 
(1938). 


R. C. Herman and R. Hofstadter, J. Chem. Phys. 6, 
534 (1938). 


mentioned above. In the present work, values of 
Av/v have been obtained for the shift of »(OH) to 
v(OHO) for the vapors of butyric and valeric 
acids. Comparing these values with those ob- 
tained for the first three members of the fatty 
acid series one finds that they are only roughly 
the same, their average value being about 0.136. 
This is in fair agreement with Av/yv calculated 
from the results of other investigations on acetic 
acid.® ® There is a frequency shift of the funda- 
mental »(OH) band in the first four members of 
the series of about 70, 63, 68 and 66 cm™ per kcal. 
of bond energy, respectively. Fig. 1 is a plot of 
Av/v vs. bond strength which includes the original 
points given by Badger and Bauer. In addition 
are included points which represent the first four 
members of the fatty acid series as well as other 
points for which we have employed bond 
strengths estimated by Pauling from latent heats 
of sublimation.’ The difficulty of making such 
estimates is shown by the different values ob- 
tained from latent heats of sublimation and 
latent heats of vaporization.’ Since we have no 


4R. C. Herman and R. Hofstadter, J. Chem. Phys. 7, 
460 (1939). 

5 A. M. Buswell, W. H. Rodebush and M. F. Roy, J. Am. 
Chem. Soc. 60, 2239 (1938). 

6M. M. Davies and G. B. B. M. Sutherland, J. Chem. 
Phys. 6, 755, 767 (1938). 

7L. C. Pauling, Nature of the Chemical Bond (Cornell 
Univ. Press, 1939). 

8 Badger and Bauer (reference 1) obtained the bond 
strength in the first members of the alcohol series by 
plotting the latent heat of vaporization against the number 
of carbon atoms in the molecule and extrapolating back to 
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we must tentatively employ the bond strengths 
that are given. An inspection of the plot shows 
that the points representing the few substances 
for which both Av/v and bond strength are known 
do not lead to any clear relationship. The absence 
of any definite regularity is indeed not at all in 
disaccord with general considerations since there 
seems to be no @ priori reason to expect a 
correlation between the bond strength and Av/v». 
If it is possible that the bond strengths for water 
and the alcohols are in error, a precise knowledge 
of the bond strengths might lead to some 
empirical relationship. It would be very desirable 
to devise some direct method of obtaining ac- 
curately the bond strengths of substances like 
water and the alcohols which are not sufficiently 
associated in the vapor phase to employ known 
methods with any degree of success. 

The present work shows that the spectroscopic 
method is a satisfactory method of determining 
bond strengths. The procedure of determining 
the amount of association at various tempera- 
tures and pressures and therefore the bond 
strength by measuring the intensity of some 
association band in the vapor phase could of 
course not be regarded as satisfactory if different 
bands led one to different values for the heat of 
association or if the absorption constant of the 


zero carbon atoms. The 7.5 kcal. of energy remaining was 
attributed to hydrogen bonding and since it was assumed 
that each molecule has two neighbors a bond strength of 
7.5 kcal. per mole per bond was obtained. However, 
according to the x-ray diffraction experiments of G. G. 
Harvey, J. Chem. Phys. 6, 111 (1938); 7, 878 (1939) both 
methyl and ethyl alcohol are not completely coordinated 
at room temperature so that the estimate of the bond 
strength would be expected to be too low. From latent 
heats of sublimation, on the other hand, Pauling (reference 
7) estimates that the bond strength in methyl and ethyl 
alcohol is 6.2 kcal. per mole per bond. Thus, from present 
evidence one can only conclude that the bond strength 
for alcohol probably lies between 6.2 and something over 
7.5 kcal. per mole per bond. Likewise, the estimate of the 
bond strength of water (about 5 kcal. per mole per bond) 
made by assuming that the latent heat of vaporization is 
almost entirely due to hydrogen bonding and that the 
number of hydrogen bonds per molecule is two, would be 
expected to be too low since Cross. Burnham and Leighton, 
J. Am. Chem. Soc. 59, 1134 (1937) have found that water 
is not completely coordinated near and above the melting 
point. In fact it has been estimated that at 40°C the average 
number of hydrogen bonds per water molecule is already 
only somewhat greater than half the maximum number 
er yore If Pauling’s (reference 7) estimate of a hydrogen 

nd strength of 4.5 kcal. per mole per bond for water is 
correct, this would imply that an appreciable part of the 
latent heat of vaporization of water is due to van der Waals 
and other forces. 
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way of knowing how reliable these estimates are, 
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band under investigation varied markedly with 
the temperature. Recently, Davies and Suther- 
land® studied the association of several carboxylic 
acids in CCl, solution by measuring the variation 
with temperature of the integral absorption of 
the monomeric v(OH) band. They obtained heats 
of association which varied from 6 to 10 kcal. per 
mole compared to the value of about 15 kcal. 
per mole known from vapor density measure- 
ments. They attribute this variation to a vari- 
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Fic. 1. Plot of Av/v for a fundamental band vs. bond 
strength in kcal. per mole. The closed circles are some of 
the original points given by Badger and Bauer. The squares 
represent points which have been plotted by using Pauling’s 
value for the bond strength, obtained from heats of sub- 
limation. The bond strength’s employed for formic and 
acetic acids are those given by the vapor density method 
while those for propionic and butyric acids have been 
obtained spectroscopically. 


ation of the absorption constant with the temper- 
ature when working with solutions. We have 
found no evidence for any variation with temper- 
ature of the absorption constant of the association 
bands that we have studied in the vapor, and 
from a study of HCOOH and HCOOD find that 
to within the experimental error different bands 
lead to the same value of the heat of association. 

During the past few years several spectroscopic 
investigations of association in the carboxylic 
acids have been made. Badger and Bauer® in- 
vestigated the variation with temperature of the 
second harmonic of the OH band in acetic acid 
and obtained rough agreement with the data of 
MacDougall’® found by the vapor density method. 
Studies made in this laboratory on the vapor of 
~ 9R. M. Badger and S. H. Bauer, J. Chem. Phys. 5, 605 


(1937). 
10 F, H. MacDougall, J. Am. Chem. Soc. 58, 2585 (1936). 
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Fic. 2. Equilibrium constant for C;sH7;COOD. 


acetic acid-d (CH;COOD) by means of the 
v(ODO) association band led to a heat of associ- 
ation of 15,900 calories per mole as compared to 
the value 16,400 calories per mole obtained by 
MacDougall.’ Apart from the work of Davies and 
Sutherland mentioned above the heat of associ- 
ation of propionic acid-d (C2:H;COOD) vapor has 
been obtained recently and found to be 14,100 
calories per mole.‘ It is the purpose of this paper 
to discuss further results that have been obtained 
for light formic acid, formic acid-d (HCOOD) 
and butyric acid-d (C;H;COOD). It seemed de- 
sirable to investigate the formic acids first 
because it would be possible to compare the 
results with Coolidge’s vapor density determi- 
nation" and second in order to see whether the 
method employed would give consistent results 
when two different bands are studied in the same 
compound.” The monomer spectra of the vapors 
of butyric and valeric acids and butyric acid-d as 
well as some results on v(CH) absorption in the 
first five members of the fatty acid series are 
reported here also. 


EXPERIMENTAL 


All the light acids employed were Kahlbaum 
products and were purified by vacuum distil- 
lation. Butyric acid-d (Cs;sH;COOD) was pre- 
pared in vacuum by mixing a slight excess of 
butyryl chloride (Eastman) with 99.6 percent 
heavy water and was purified by pumping and 


1 A. S. Coolidge, J. Am. Chem. Soc. 50, 2166 (1928). 


1 These are the same compounds for our purposes use 
the experimental error is too large to show any difference 
due to zero point energy. 


several distillations. Only middle fractions were 
retained for this as well as the other compounds 
employed. The vapor pressure of this compound 
as well as its infra-red spectrum indicated that 
less than several percent of water was present. 
The hydrogen-deuterium gas mixture obtained 
from the reaction of the heavy butyric acid with 
magnesium was analyzed in a mass spectrograph 
and the resulting D/H ratio was 40.1* Formic 
acid-d was prepared by exchanging light formic 
acid (Kahlbaum) with 99.6 percent heavy water 
several times. Between successive exchanges the 
water was removed by means of anhydrous 
copper sulphate. This method had been employed 
previously by Hofstadter in this laboratory.'4 In 
the attempt to prepare HCOOD with a greater 
D/H ratio than had been obtained previously, 
the anhydrous copper sulphate employed was 
initially exchanged with heavy water. Although 
the four exchanges carried out should have 
yielded a product with a D/H ratio of about 140, 
the final sample had a ratio of 16. Because a final 
dehydration had to be performed the sample 
eventually employed in the infra-red experiments 
had a D/H ratio of 12. The first of these D/H 
ratios was determined in the same way as that 
mentioned above for butyric acid-d. The other 
ratio was obtained from a comparison of the 
depths of the OH and OD bands found in a high 
temperature and pressure spectrum. This method 
has been employed by Hofstadter who found 
good agreement between the D/H ratio found 
in this way and by means of the mass spectro- 
graph. That both the light and heavy formic 
acids were free from water was known from a 
comparison of their vapor pressures with data 
obtained by Coolidge." 

The experimental arrangement employed to 
determine the association properties of the formic 
acids was the same as that used for acetic acid-d* 
and the determination of the heat of association 
of butyric acid-d was made with the same 
apparatus employed for propionic acid-d in this 
laboratory.‘ Different arrangements had to be 
employed because of the very different vapor 
pressure curves of these compounds. Intensities 


13 We wish to thank Mr. W. Guyer for analyzing our 
hydrogen-deuterium gas mixtures. 

4 R. Hofstadter, J. Chem. Phys. 6, 540 (1938). 

16 A. S. Coolidge, J. Am. Chem. Soc. 52, 1874 (1930). 
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were measured at various pressures and tempera- 
tures at the bottom of the O—D—O association 
band for HCOOD and C;H;COOD whereas for 
HCOOH the association band at 7.384 was used 
because the O—H—O association band falls in 
the region of C—H valence vibrations. Several 
times during a determination, the transmission 
was obtained with the sample frozen out of the 
absorption cell in order to check background 
absorption which remained constant. It was 
found necessary, as usual, to wait a considerable 
length of time between points in order for 
equilibrium to set in. Otherwise, points on the 
log K, vs. 1/T curve were found to scatter badly. 
In all the association measurements, pressures 
were measured directly with a mercury manome- 
ter and temperatures by means of a calibrated 
thermocouple. However, the pressures noted in 
the infra-red spectra shown in Figs. 5 and 6 were 
calculated from the temperature of the liquid in 
a side tube, except in the case of butyric acid-d 
(Fig. 5) where it was measured directly. The 
general technique employed will not be discussed 
here since it has been reported in previous 
communications.* 4 


RESULTS AND DISCUSSION 


As in the determination of the heat of associ- 
ation of propionic acid-d,* the absorption con- 
stant for the O—D—O association band in 
butyric acid-d was assumed to be the same as 
that in acetic acid-d. In a previous communication 
it was remarked that the addition of a CHe group 
in going from one fatty acid to the next would 
probably not affect appreciably the absorption 
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Fig. 3. Equilibrium constant for HCOOH. 
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Fic. 4. Equilibrium constant for HCOOD. 


constant of the O— D—O band.*‘ Fig. 2 shows the 
results obtained for butyric acid-d on the basis of 
this assumption. From the slope of the best 
straight line the heat of association obtained is 
13,800+700 calories per mole. This value is 
practically the same as that obtained for propionic 
acid-d. Calculations were made for this acid 
choosing an absorption constant differing by 20 


' percent from that above. The fact that there were 


no noticeable deviations from linearity of the 
curve on the log K, vs. 1/T plot shows that both 
the heat of association and the absorption con- 
stant are not in error by more than about 5 and 
20 percent, respectively.'® Moreover, over the 
higher temperature range pp/p (where pp is the 
partial pressure of double molecules and p the 
total pressure) was kept practically constant so 
that the slope of the curve in that region is 
practically independent of the choice of absorp- 
tion constant.!® 

Light formic acid was investigated in order to 
see whether this method would give results in 
agreement with the vapor density measurements 
obtained by Coolidge. The 7.384 band which was 
employed is a dimer band and is not mixed with 
any monomer band as can be seen from an 
inspection of the monomer and dimer spectra 
of this compound obtained by Bonner and 
Hofstadter.? Since Coolidge’s data were avail- 
able, an absorption constant was determined for 
this band in the same way that the absorption 
constant for the O—D—O band in acetic acid-d 
was determined from MacDougall’s data. The 


16R. C. Herman and R. Hofstadter, J. Chem. Phys. 7, 
630 (1939). 
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Fic. 5. spectra of C;H7COOH and C;H7;COOD. 
Absorption cell length 22 cm. 


value of the absorption constant for the 7.38 
band in HCOOH obtained in this way is 5104 
cm? per mole.” It should be remarked that by a 


17 All the absorption constants quoted in this paper have 
been obtained from the formula 
where pp is the density in moles per cm’ and / the length 
of the absorption path in cm. 
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method of trial and error one could find the 
value of the absorption constant for which the 
data on the log K, vs. 1/T plot would give a best 
straight line. This is effectively what was done in 
the case of butyric acid-d. The association data 
for HCOOH are shown in Fig. 3 and the heat of 
association obtained from the slope is 12,400+700 
calories per mole as compared with Coolidge’s 
value of 14,135 calories per mole. These two 
values differ by an amount larger than the 
experimental error and it was this discrepancy 
which in part prompted the further experiments 
on HCOOD for which a different band (the 
O—D-—O association band) could be employed. 
An absorption constant for the O—D—O band 
was obtained in the same way as for the 7.38u 
band in HCOOH and is 7.4104 cm? per mole.” 
Since there was some light formic acid present, 
the vapor was assumed to contain an equilibrium 
mixture of two types of single molecules and 
three types of double molecules. If one takes this 
into consideration it is simple to show that 


[1—(1+1/f)x}?p 
(1+1/f)x 
RT logio I/Io 


8 is the absorption constant, p the total pressure 
and f the D/H ratio in the acid position. We 
should mention that in obtaining B for HCOOD 
we assumed that double molecules with one D 


where 


x=pp/p=— 


TABLE I. 
HEAT OF ASSOCIATION 
AcIp INVESTIGATOR METHOD IN CAL./MOLE 
HCOOH®* Coolidge" V.D.* 14,135 
Ramsperger and Porter V.D. 14,100 (av. value) 
this paper an 12,400 
HCOOD this paper S 12 800 
CH;COOH* MacDougall? V.D. 16,400 
Nernst and von Wartenburg? V.D. 15,000 
Fenton and Garner! V.D. 13, 790 
Badger and Bauer® S(Harmonic) rough agreement 
with MacDougall 
Davies and Sutherland*® S(CCl, sol.) 9300** 
CH;COOD Herman and Hofstadter*® S 15,900 
C:H;COOD Herman and Hofstadter S 14,100 
C;H;COOD this paper S 13,800 
C.H:;COOH Fenton and Garner!® V.D. 7050*** 


* 
density method. 


nreliable. 
: ome other work on formic acid see H, C 


4 Nernst and von Wartenburg, Zeits. f. Electrochem. ba 37 (191 


troscopic method employing the fundamental band of a vapor unless stated otherwise. 


. Ramsperger and C. W. Porter, J. Am. Chem. Soc. 50, 3036 (1928). 
. Ramsperger and C. W. Porter, J. ‘Am. Chem. Soc. 48, 1267 (1926). 

° ° Hue earlier work on the vapor density of acetic acid see: W. Ramsay and S. 
Soc. 12, 374 (1909-1910); Drucker and Ullman, Zeits. f. physik. me 74, 604 ‘apie: Holland, Zeits. f. 


ng, J. Chem. Soc. 49, 790 (1886): S. Young, Proc. Roy. Dublin 
Electrochem. 18, 234 (1913). 
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atom and one H atom in the dimer ring contribute 
half the amount that double molecules with two 
D atoms contribute to the intensity of the 
O-—D-—O band. From the results for HCOOD 
shown in Fig. 4 we obtain a heat of association 
of 12,800+700 calories per mole. 

Although the formic acids in the absorption 
cells were examined over a period of days at 
temperatures as high as 150°C, no decomposition 
was observed, so that the pressures recorded 
were those of the acid alone. The fact that the 
heavy formic acid was not completely deuterated 
could be taken into account as mentioned above. 
From a plot of log J/Io vs. pp (the density of 
dimer molecules calculated by means of Coolidge’s 
data) it was evident that Beer’s law was obeyed 
for all but some of the points employed on the 
low temperature part of the curves shown in 
Figs. 3 and 4. Furthermore, the absorption con- 
stants were not found to depend upon the 
temperature in the range studied. The O—D—O 
band in HCOOD was mixed with a weak 
monomer band due to HCOOH. Under the con- 
ditions employed this monomer band was negli- 
gibly small. The presence of a monomer band, 
however, would cause the intensity of the dimer 
band to decrease more slowly than for the dimer 
band alone and would lead to a higher rather than 
a lower heat of association. Thus, the effect of the 
8.304 band on the long wave-length side of the 
7.384 band studied in light formic acid and the 
small monomer band superposed on the O—D—O 
band in heavy formic acid does not explain the 
lower heats of association obtained by this 
method. We do not think that the other factors 


TABLE II. Positions of monomer absorption bands. 


VALERIC ACID Butyric Acip Butyric Acip-d 
(cm~) A(u) (cm™) 
2.70 3703 2.80 3571 3.33 3003 
3.33 3003 3.31 3021 3.70 2703 
5.63 1770 5.65 1770 5.68 1760 
6.80 1471 6.91 1447 6.90 1449 
7.30 1370 7.30 1370 7.32 1366 
8.00 1250 7.90 1266 7.90 1266 
8.80 1136 8.12 1232 8.40* 1190 
9.17 1091 8.77 1140 9.25 1081 

10.5 952 9.23 1083 10.1 990 
12.2 820 11. 909 
12. 833 


wll My 


== 


Fic. 6. Absorption spectrum of C,H,COOH. Absorption 
cell length 22 cm. 


considered; namely, slit width (we employed 
constant slit width in each case) or adsorption 
account for this discrepancy either. Because of 
the care with which Coolidge’s determination 
was made it seems that the discrepancy resulting 
from these two different methods may be real. 
Perhaps it is‘caused by the fact that the vapor 
density method includes all types of association 
whereas this method does not. Or it may be due 
to the presence of higher polymers. Because of the 
scattering of points in Figs. 3 and 4 it is difficult 
to say exactly what this difference is but from our 
results on acetic acid it does not seem that our 
error would be as large as 1500 calories per mole. 

As far as we are aware the heats of association 
of only very few of the fatty acids have been 
obtained by the vapor density method, namely, 
formic, acetic and heptoic acids. In Table I we 
have gathered some of the results obtained by 
this method as well as others obtained spectro- 
scopically. Fenton and Garner" state that 
although the values they obtained for the heat of 
association of heptoic acid do not possess much 
value by themselves, the results indicate that the 
heat of association of this acid is decidedly lower 
than that of acetic acid and suggest the investi- 
gation of acids with intermediate chain length. 
From our results there is no indication that the 
bond strength is decreasing with chain length as 
far up the series as butyric acid. It should be 
interesting to investigate carefully a member 
higher in the series to see whether the bond 


* Asymmetry on long wave-length side. 


#8 Fenton and Garner, J. Chem. Soc. 694 (1930). 
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NUMBER OF C ATOMS 
Fic. 7. Plot of the absorption constant for the band 


v(CH) vs. the number of carbon atoms in each fatty acid 
molecule (reference 22). 


strength remains roughly constant as _ this 
investigation shows or decreases as Fenton and 
Garner’s work on heptoic acid indicates. 

The monomer spectra of the vapors of valeric 
acid and butyric acid and butyric acid-d are 
shown in Figs. 5 and 6 and band centers are 
listed in Table II. Since the degree of dissociation 
is rather large for these compounds under con- 
ditions for which the vapor pressure is high 
enough to obtain sufficient absorption, it was not 
possible to obtain good dimer spectra. As in the 
other fatty acids studied the high frequency 
region is interpretable. An inspection of the 
spectra in this and earlier papers?-* 4 shows that 
the position of the »(C=0) frequency in the 
carboxyl group is practically independent of the 
length of the attached hydrocarbon chain. This 
result has been obtained previously by means of 
the Raman effect.’ The question of assigning 
bands beyond »(C=0) must be left until further 
work is done with higher resolution. However, it 
is likely that the band appearing near 6.9y in 
these higher acids is due to some hydrocarbon 
chain frequency and perhaps the band in the 
vicinity of 10.54 in the spectra of the higher 
members of the series which have D atoms in 
the acid position, contains the OD bending 


19K. W. F. Kohlrausch, F. Koppl and A. Pongratz, 
Zeits. f. physik. Chemie B21, 242 (1933); B22, 359 (1933). 


frequency which appears at about 114 in heavy 
phenol (CsH;OD).?° 

If one studies the monomer spectra of the fatty 
acids with progressively longer hydrocarbon 
chains it is evident that the band associated with 
CH valence vibrations becomes deeper. We have 
obtained absorption constants for this band in 
the first five members of the series by measuring 
the depth of the band below the background with 
a wide slit. The following relation holds ap- 
proximately for all but the first member : 


By = Bacetic+0.77(N — 2) 


where JN is the number of carbon atoms in the 
molecule and Bacetic is about 3.2104 cm? per 
mole. We also find that the »(CH) absorption 
constants for formic and acetic acids are about 
the same. These results are shown in Fig. 7. It 
seems then that each CHe group added to the 
hydrocarbon chain increases CH absorption by 
about the same amount and also that the 
contribution to CH absorption of the CH group 
in formic acid is about the same as that of the 
CH; group in acetic acid. Similar results have 
been found in the case of the alcohols by Henri. 
These results are probably only approximately 
correct since Bonino in the case of the alcohols 
has shown that as one goes to higher members of 
the series, the deviations from linearity become 
greater and greater.”! 
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20V. Williams, R. Hofstadter and R. C. Herman, J. 
Chem. Phys. 7, 802 (1939). : 

21 A complete account of this point may be found in 


—— and Matossi, Das Ultrarote Spektrum (1930), p. 
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Reactions of Free Radicals with Organic Compounds Containing Atoms with 
Unshared Electron Pairs* 
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An experimental test for the reaction CH; +CH;OCH.zCH;~CH;0CH;+CH;CH; was made 


VOLUME 8 


by decomposing methyl ethyl ether both alone and in presence of azomethane, and examining 
the products for dimethyl ether; none was detected, and if this reaction occurs at all it is only 
to the extent of a very few percent. A similar test for the reaction CH;+CH;CH2CH2NH:2 
—CH;NH:+CH;CH:CH: was made by decomposing pure -propylamine in a flowing system, 
but this also gave a negative result. The products formed in the promoted decomposition of 
dimethyl ether were found to be essentially the same as those formed in the normal decom- 


position. 


1. INTRODUCTION 


N a recent paper,! the whole question of the 

nature of the elementary reactions between 
radicals and molecules was considered, and the 
conclusion was reached that, in general, the most 
favored reaction would be removal of a hydrogen 
atom from the organic molecule by the free 
radical. On the other hand, the reaction that 
might be expected to have the smallest proba- 
bility of occurring is the reaction between a free 
radical and a carbon atom having four valences 
joined to four univalent atoms or groups. Such a 
“shielded” carbon atom would be difficult to 
attack, and the reaction would be by far the 
least likely sort to occur. 

However, atoms like chlorine, oxygen, or 
nitrogen are not shielded to the same extent as 
carbon, all having at least one pair of unbonded 
electrons. Since we know that free radicals 
attack doubly and triply bonded carbon atoms, 
it appeared desirable to investigate experi- 
mentally whether any similar sort of reaction 
occurs between radicals and atoms having one or 
more unshared electron pairs. . 

If free radicals react with an atom having an 
exposed pair of electrons, we would expect that 
introduction of methyl radicals into an ether, 
amine, or alkyl chloride would lead to the 
production of the corresponding methy! deriva- 


*In part, from a dissertation submitted by W. D. 
Walters to the Johns Hopkins University in partial fulfill- 
= < oe requirements for the degree of Doctor of 

ilosophy. 

Wm. R. Warner Fellow in Chemistry, 1934-1938. 

sent address, University of Rochester. 

t Hormal Foundation Fellow, 1934-1938. 
Rice and Teller, J. Chem. Phys. 6, 489 (1938). 
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tives, according to the following equations: 


CH;+RCI-CH;CI+R. 


In view of the importance of this question, we 
attempted to obtain experimental evidence for 
reactions of this type by investigating three 
decompositions, namely the ordinary decompo- 
sition of methyl ethyl ether, the promoted 
decomposition of this compound, using azo- 
methane as promoter, and the ordinary decompo- 
sition of normal propylamine. 


2. EXPERIMENTAL 


The method used in this investigation was to 
introduce a known amount of substance into a 
heated reaction vessel until the desired amount of 
decomposition, as measured by the pressure 


‘change on the manometer, had occurred. The 


sample was then quickly removed from the 
furnace, condensed at liquid nitrogen tempera- 
ture, and separated into various fractions by 
distillation in a Podbielniak still. In order to 
carry out such a distillation, a larger sample is 
needed than can be contained in the small 
reaction vessel ordinarily used for rate measure- 
ments. Therefore a large furnace was constructed 
containing a 12-liter flask surrounded by two 
graphite hemispheres approximately 5 cm thick, 
contained in an oven made of 3” transite board, 
and heated by inside windings of 16-gage 
Nichrome wire. The oven was insulated by 
Sil-O-Cel brick. We found that under these 
conditions, with a working temperature range of 
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THE RMOCOUPLE 
A = 
MANOMETER: 
M. 
FLASK 
TO DIFFUSION AND 
VACUUM PUMPS 
POOBIELNIAK 
COLUMN 


TO COMPRESSION SIDE 
OF OIL PUMP 


450-550°C, there was very little oxidation of the 
graphite hemispheres. The temperature was 
measured by means of a Chrome-Alumel thermo- 
couple introduced into a well in the reaction 
vessel; the thermocouple checked to +1°C, at 
the boiling point of sulfur. The reaction flask 
was connected directly to the distillation appa- 
ratus so that a sample could be readily removed 
and separated into fractions. 

In Fig. 1 we show a sketch of the complete 
apparatus with the exception of the furnace and 
the reaction vessel. At the beginning of a run, the 
reaction flask was pumped out and brought up 
to temperature, and the substrate was introduced 
from the small flask A. When we were studying 
the promoted reaction with methyl ethyl ether, 
it was necessary to make up mixtures of known 
composition in flask A and evaporate the whole 
of the material into the reaction flask to prevent 
error due to selective evaporation. We could 
follow the course of the reaction by means of the 
manometer M,, and at any selected point during 
the reaction, could quickly pump out the con- 
tents of the reaction vessel by connecting it to 
flask A, cooled in liquid nitrogen, and connecting 
flask A in turn to the flasks V;, V2, V3, which 
could be used as a sort of Toepler pump. We did 
not try to remove more than about 95 percent of 
the contents of the reaction flask, and simply 
applied a correction for the amount not removed 
which we could determine from the reading of the 
manometer M,. The contents of the small flask 
A, and of the bulbs Vi, Ve, Vs, were then 
transferred to the flask B and distilled through 
the Podbielniak column, the products being 


collected again in the bulbs Vi, V2, Vs; these 
bulbs were in a glass box and were calibrated so 
that a measurement of the pressure and the 
temperature gave the molar contents. The results 
of the Podbielniak distillation were supplemented 
when necessary by gas analyses. 

The decomposition of m-propylamine was 
studied by a flowing method in an apparatus 
described in a previous paper.” The apparatus 
used permitted us to make a complete over-all 
analysis, but we did not do this since we 
were interested primarily in finding out if any 
methylamine was present. Accordingly, we 
simply distilled the contents of the liquid-air 
trap, using a Podbielniak column, and looked for 
indications of a break in the distillation curve at 
the boiling point of methylamine. 


3. DECOMPOSITION OF METHYL ETHYL ETHER 


The methyl ethyl ether used in these experi- 
ments was a commercial sample which was 
distilled carefully in the Podbielniak column and 
stored at —80°C until needed. The results of 
experiments 1 and 2, which were carried out with 
pure methyl ethyl ether, are summarized in 
Table I.* Considerable amounts of aldehydes are 
formed as intermediates, and some formaldehyde 
polymerizes in the first trap A;. We placed some 
aniline in the distillation flask B in order to keep 
the aldehydes in this flask and to prevent 


TABLE I. Decomposition of pure methyl ethyl ether. 


Exp. TEMP. In. PREss. TIME. % PREss. 
No. Ss MM MIN. NCREASE 
1 448 596 40 20.0 
2 473 401 63 50.5 


Products found by Podbielniak distillation. The three fractions were 
not analyzed; fraction 1 should contain CH«, CO, and Hz; fraction 2, 
C2He and CsH,; fraction 3, CsHs, CaHio, and CHsOCHs. 


pa OF FRACTION 1 FRACTION 2 FRACTION 3 

CH;0C2Hs |B. p.<—157°C|B. p.—90+45°C| B. P.—73 TO 

Exp.| TAKEN CC cc AT cc AT +7°C cc AT 
No. at NTP NTP NTP NTP 
1 3236 609.5 401.1 113.2 
2 2160 1078 358.0 90.7 


2 Rice, Ruoff, and Rodowskas, J. Am. Chem. Soc. 60, 
955 (1938). 

3 These results may be compared with those of previous 
workers who analyzed the gaseous products. See Glass and 
Hinshelwood, J. Chem. Soc. 1804 (1929); Ure and Young, 
J. Phys. Chem. 37, 1169 (1933). 
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TABLE II. Induced decomposition of methyl ethyl ether using azomethane. Experiment 4. 


VoL. VoL. 

TEMP. ETHER CHsNNCHs In. PREss. PREss. TIME 
a cc NTP cc NTP CH NCH: MM NCREASE MIN. 
297 3351 19.1 0.57 473.5 16.4 92 

Results of distillation and analysis of products cc at NTP 
CO CH, CoH, CH;0C;:H; 
124.7 44.3 347 10.2 94.4 2878 


* Ether decomposed =3351—2878 =473 cc at NTP. 


>From the above data we estimate that there were about 400 cc (at NTP) of aldehydes and higher boiling compounds that were not de- 


composed. 


TABLE III. Induced decomposition of methyl ethyl ether using azomethane. Experiment 5. 


VoL. 

Temp. oumen CHINN Hs % In. PREss. % PRESS. TIME 
~~ cc NTP cc NTP CH:NNCH; MM NCREASE MIN. 
300 2884 102.3 3.43 426.0 55.8 280 

Results of distillation and analysis of products cc at NTP 
co CH, CoH, CoH. CH;0C;3H,; 
630 132 1347 5.9 264 1553 


* Ether decomposed =2884 — 1553 =1331 cc at NTP. 


>From the above data we estimate that there were about 600 cc (at NTP) of aldehydes and higher boiling compounds that were not de- 


composed. 


them from polymerizing and clogging up the 
Podbielniak column. 

There was no break whatsoever in the distil- 
lation curve at the boiling point of dimethyl 
ether, and the whole fraction collected between 
—73 and 7°C amounted to barely 10 percent of 
the total products. Since there is the possibility 
that in these experiments we are not dealing with 
radical reactions, or at least, that the radical 
chain‘ occurs only to a small extent, we decided 
to study the promoted decomposition of methyl 
ethyl ether under conditions where the pure 
substance is stable. 


4. PROMOTED DECOMPOSITION OF METHYL 
EtHuyL ETHER 


We first made a run with pure methyl ethyl 
ether at 299°C and an initial pressure of 410 mm. 
After 277 minutes, the pressure was 413.4 mm, 
so that there was only 0.8 percent pressure 
increase. We then performed two experiments at 
approximately the same temperature by using a 
mixture of methyl ethyl ether and azomethane. 
As in runs 1 and 2, we distilled the products, but 
in these experiments we also made a complete 


‘See in this connection Staveley and Hinshelwood, Proc. 
Roy. Soc. A159, 198 (1937), who have shown that at 560°C 
the reaction is inhibited by nitric oxide and that the 
minimum chain length is 7.5. 


over-all analysis of all the fractions; the results 
are summarized in Tables II and III. In none of 
the experiments was there any break in the 
distillation curve at the boiling point (— 24°C) of 
dimethyl ether. While we did not look so carefully 
for a break in the case of propane and butane, the 
distillation curve showed no trace of inflection at 
these boiling points either. The absence of these 
compounds indicates that radicals attack neither 
the carbon nor oxygen atoms of the ether. Both 
acetaldehyde and formaldehyde were formed in 
the promoted decomposition ; in a fraction in the 
correct boiling range, the presence of acetaldehyde 
was proved by its characteristic odor and a test 
with Schiff’s reagent ; formaldehyde appeared as 
its polymer, which was identified by its odor and 
by the red color produced by its solution in the 
presence of resorcinol and sodium hydroxide. 

As in experiments 1 and 2, there was a small 
fraction of gas (90-100 cc) collected over the 
temperature range —73 to +7°C. This fraction 
appeared to be more or less independent of the 
size of the original sample, as well as of the 
amount decomposing. In view of the appearance 
of these samples which constitute a small, but 
appreciable amount of the products, we decided 
to make up known mixtures of ethane, methyl 
ethyl ether, and dimethyl ether, with the object 
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TABLE IV. Thermal decomposition of n-propylamine. 


Moles passed through furnace (39.3 g) 
Moles decomposed 
Percent decomposed 
Temperature, °C 
Pressure, mm 
Time of run (minutes) 
Moles of gas not condensed by liquid air 
3 
CH;NH:; 


of determining the smallest amount of dimethyl 
ether that we could detect under conditions 
similar to those in the distillation of the products 
of decomposition. We found that we could just 
detect 47.5 cc of dimethyl ether, which gave a 
slight break in the distillation curve, whereas 
with 97 cc we could obtain a definite and true 
boiling point. 

In experiment 5, 1332 cc of methyl ethyl ether 
were decomposed, and 5 percent should give a 
detectable break in the distillation curve. In 
experiment 4, where we decomposed only 472 cc 
of methyl ethyl ether, we can conclude that any 
production of dimethyl ether was less than 10 
percent of the amount of methyl ethyl ether 
decomposed. 

Assuming that each molecule of azomethane 
gives two methyl radicals, the number of ether 
molecules decomposed by the addition of each 
radical in experiment 4 is 17 and in experiment 
5 is 7. Although no special weight can be given 
to our two results, it is interesting to note that 
they indicate that the chain length at a given tem- 
perature is inversely proportional to the square 
root of the concentration of the azomethane. 


5. DECOMPOSITION OF n-PROPYLAMINE 


The decomposition of propylamine has been 
investigated by Sickman and O. K. Rice.’ They 
concluded from a study of the decomposition at 


&’Sickman and O. K. Rice, J. Am. Chem. Soc. 57, 22 
(1935). 
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low pressures that the reaction is probably a 
chain. The decomposition of ethylamine pro- 
ceeds® through the intermediate formation 
of diethylhydrazine according to the equation 
However the de- 
composition of diethylhydrazine was shown to be 
a radical chain, so that we may infer that even if 
propylamine follows a similar course, free propyl 
radicals will be formed at some step in the reac- 
tion and will decompose into methyl radicals and 
ethylene. 

In our experiments, we used a commercial 
sample of propylamine which was carefully dis- 
tilled, and the middle fraction, boiling over a 
range of 0.5°C was used. The -propylamine was 
thermally decomposed to the extent of 2-3 
percent at 650°C and 10 mm pressure and the 
results obtained are shown in Table IV. Since we 
were primarily interested in looking for the 
presence of lower amines, we did not make a 
complete over-all analysis of the products formed. 
There was no indication whatsoever in the 
distillation curve of the products of a methylamine 
plateau at —7°C, or of an ethylamine plateau 
at 17°C. 

While our results indicate that reactions of the 


type 
R CH +CH 3CH.CH, 


either do not occur or are followed by re- 
dissociation into the original components, it 
should be pointed out that there is evidence for 
attack by methyl radicals on the nitrogen atoms 
in simple azocompounds’ according to the 
equation 


CHs+ CHs— N= = N—-CHs; 


followed by further reaction with radicals to 
give higher molecular weight compounds. 


a 936) A. Taylor and Ditman, J. Chem. Phys. 4, 212 

7 See in this connection Davis, Jahn and Burton, J. Am. 
Chem. Soc. 60, 10 (1938); Riblett and Rubin, ibid. 59, 
1537 (1937). 
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On the Predissociation of Several Diatomic Molecules 
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From the observed anomalous intensity distribution in 
the S» bands it is concluded that the dissociation limit for 
the first predissociation of S. lies appreciably below the pre- 
dissociation limit (4.4 ev). It is found that the anomalous 
intensity as well as the diffuseness persists even at very 
low pressures (0.01 mm) and that therefore this is a case 
of genuine predissociation and not pressure broadening 
as assumed by Lochte-Holtgreven. The second predis- 
sociation of S, starting with the 18—0 band leads to an 
upper limit of 3.6 ev for the heat of dissociation of So. 
An anomalous intensity distribution similar to that of S» 
is found for NO assuming that the e absorption bands are 
a continuation of the y bands. This indicates a predissocia- 
tion of NO in the upper 22 state of the y bands into normal 
atoms with a limit 1.3 ev above the dissociation limit. 
The remarkable intensity distribution found by Schiiler 


and co-workers in the CuH emission band \4280 at low 
pressure is explained as due to a predissociation into the 
82+ state derived from normal Cu(2S) and normal H(2S). 
Such a state cannot cause the predissociation of the level 
J’=0 in agreement with observations. The heat of dis- 
sociation of CuH comes out to be <23,325 cm™ (i.e., 
<2.892 ev). The method previously suggested by one of 
the authors for the detection of maxima of potential 
curves is applied to AID (AIH) and BH. From the break- 
ing off points observed in the visible bands of these 
molecules it is established that their 'Il states have po- 
tential curves with maxima. Therefore the dissociation 
energies are appreciably lower than the energies of the 
observed breaking off points, i.e., Do(AIH) <24,775 cm™, 
Do(AID) <25,010 Do(BH) <28,350 


A. DETECTION OF PREDISSOCIATION IN So AND 
NO sy ANOMALOUS INTENSITY DiIsTRI- 
BUTION IN ABSORPTION 


T is well known that a diffuse absorption line 
(or band) under small dispersion appears 
with greater apparent intensity than a sharp line 
(or band) of the same true intensity. Therefore, 
if the true intensity is known in some way, a 
diffuseness may be detected indirectly by an 
anomalously large apparent intensity even if it is 
too slight for direct detection by the spectrograph 
used. This method has been applied, for example, 
in a number of cases to the detection of induced 
predissociation, where it is found that certain 
bands for the same number of absorbing mole- 
cules (that is, the same true intensity) appear 
more intense at higher pressure due to a predis- 
sociation (broadening) induced by collisions. In 
only one case, it seems, has this method of 
detection as yet been applied to ordinary pre- 
dissociation, namely in the case of the D system 
of the Hz molecule by Beutler, Deubner and 
Jiinger.! This system exhibits in absorption only 
the bands with v’=3 whereas those with v’ <3 are 
so weak that they are not observed. Since the 


* Present address: Department of Physics, Purdue 
University, Lafayette, Indiana. 

‘H. Beutler, A. Deubner and H. O. Jiinger, Zeits. f. 
Physik 98, 181 (1935). 


true intensity of the 2—0 and 3—0 bands must 
be of the same order of magnitude the sudden 
jump in apparent intensity can only be explained 
by the assumption that the bands with v’=3 are 
broadened in consequence of predissociation even 
though no visible broadening is found. This 
conclusion is confirmed by the absence of emis- 
sion bands with v’=3. 

Since, as in the above example, this method 
of detecting ordinary predissociation by anoma- 
lously large intensity in absorption is considerably 
more sensitive than the method employing an 
actual diffuseness of the bands it should prove 
useful also in other cases particularly when 
extremely large dispersion is not available. 

In the following we shall describe two further 
cases in which the method has been applied, 
namely the absorption spectra of S, and NO. 


(a) Se 


We have reinvestigated the absorption spec- 
trum of diatomic sulphur vapor with a very long 
absorbing path (6 meters) hoping to find new 
weak systems of absorption bands (forbidden 
transitions) at longer wave-lengths than the well- 
known main system in the ultraviolet. Although 
we did find such absorption bands, extending as 
far as 8700A their analysis, due to small dis- 
persion, has not yet been possible. However we 
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Fic. 1. Absorption spectrum of the S, molecule at low pressure (about 0.05 mm). 
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Fic. 2. Estimated intensities in the absorption band progression v’’=0 of Sy. 


also photographed the main absorption system 
with the long path at correspondingly low pres- 
sures and obtained some results that, we believe, 
are of interest particularly in view of recent 
controversies concerning the predissociation of 
the molecule. 

As is well known, there are two points in the 
main absorption system of Se where the bands 
become diffuse: the first beginning with the 
11—0 band and the second with the 18—0 band. 
This is shown by the spectrogram reproduced in 
Fig. 1. The fine structure lines of the 11—0 band 
appear distinctly broadened compared to those 
of the 10—0 band. Such a broadening is also 
present in the 12—0, 13—0, --- 17—0 bands 
although to a lesser and lesser extent. The 18—0 
band and the subsequent bands however are 
quite difftise, no fine structure at all being 
noticeable. 

While these facts concerning the predissoci- 
ation of Ss have been known for a long time, 
another point, which until now seems to have 
escaped consideration, is brought out by our 
spectrogram Fig. 1: The intensity distribution in 
the series of bands is quite different from what 
one would expect from the Franck-Condon 
principle. Since all the bands have v’’=0 the 
intensity should first rise gradually and smoothly 
and then drop again in the same manner. 
Actually there is a sudden jump in the intensity 
at the 10—O0 band: the bands from 10—0 on are 
very strong whereas the preceding bands are 
quite weak. Within each group of bands, 
those with v’<10 and those with v’=10, the 
intensity does not vary very greatly. Fig. 2 
represents graphically the intensity distribution 
from estimated intensities. 

It seems to us that the sudden intensity jump 
in the progression of S. absorption bands can be 
explained only by the assumption that this is a 
case, similar to the one observed for He (see 
above), in which bands appear with an anoma- 
lously large intensity because of diffuseness of 
the fine structure lines. With the dispersion used 
the 10—0 band appears to have a sharp fine 
structure whereas the following bands are visibly 
diffuse. However from the large intensity of the 
10—0 band we conclude that a diffuseness is 
present in this band as well, although to a lesser 
extent, and that therefore predissociation sets in 
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at the v'=10 level and not at v’=11. This con- 
clusion is in agreement with results of Olsson? who 
investigated the spectrum with very large dis- 
persion and found the 10—0 band to have 
slightly diffuse lines. 

It is significant that as shown by the spectro- 
gram Fig. 1 the fine structure lines of the 9—0 
band appear much more intense at some distance 
from the head than at the head. This distribution 
is different from that in any other band and 
indicates that the higher rotational lines of this 
band are also broadened, i.e., that predissociation 
occurs for the higher rotational levels of the 
v’=9 state. 

In this connection it should be mentioned that, 
in emission, only bands with v’=9 occur (see 
Olsson?) and that the bands with v’=9 break off 
abruptly.’ This confirms the conclusion that the 
limit for this particular predissociation lies be- 
tween some high rotational level of the v’=9 
state and the lowest rotational level of the 
v’=10 state.* 

There are essentially two cases of predis- 
sociation that consist of radiationless electronic 
transitions from a discrete to a continuous state. 
In the first case (Ib*) the point of intersection of 
the potential curves of the two electronic states 
in question lies below the asymptote of the 
potential curve of the state causing the predis- 
sociation (see Fig. 3); in the second case (Ic) the 
point of intersection lies above the asymptote 
(see Fig. 4). In the first case predissociation sets 
in with maximum probability at a point corre- 
sponding to the asymptote (dissociation limit). 
In the second case the maximum probability of 
the radiationless transition corresponds to the 
point of intersection of the potential curves which 


*E. Olsson, Zeits. f. Physik 108, 40 (1937); also Thesis, 
Stockholm, 1938. We had reached the above conclusion 
before this paper appeared but unfortunately publication 
of our results was greatly delayed. 

*G, Herzberg, Ann. d. Physik 15, 701 (1932). 

‘Rosen, Désirant and Duchesne (Phys. Rev. 48, 916 
(1935)) and Rosen and Bouffieux (Bull. Acad. Belg. 1936, 
p. 885) have observed a termination of the fine structure 
of emission bands even with v’=8 at a higher value of the 
rotational quantum number (see however Olsson, reference 
2, and further below). 

*See G. Herzberg, Ergebn. exact. Naturwiss. 10, 207 
(1931) ; for a slight revision of the nomenclature concerning 
the subcases of case I of predissociation (radiationless elec- 
tronic transition) see G. Herzberg, Molecular Spectra and 
Molecular Structure I. Diatomic Molecules (Prentice-Hall, 
— York, 1939). Case Ja is the limiting case between 6 
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is above the asymptote of one of them; but the 
predissociation will set in, though with smaller 
probability, somewhat below this point in conse- 
quence of the quantum-mechanical penetration of 
potential barriers. To be sure, this effect is very 
slight indeed for molecules as heavy as Se, except 
when the potential hill formed by the two inter- 
secting potential curves is very narrow. 

In the first predissociation of Se, however, as 
shown above, the level v’ = 10 is not predissociating 
as strongly as the level v’ = 11, and therefore this 
predissociation does represent an example of case 
Ic with a rather narrow potential hill. Thus we 
conclude that the dissociation limit corresponding 
to this predissociation (i.e., the asymptote of the 
potential curve of the state causing the predis- 
sociation) lies appreciably below the predissociation 
limit. 

This conclusion is the opposite to the one 
obtained by Rosen, Désirant and Duchesne‘ who, 
from the observed breaking off in the emission 
bands with v’ =8 and 9, according to a method of 
one of us,® concluded that this predissociation 
belongs to case Jb, i.e., that the predissociation 
limit coincides with the corresponding dissoci- 
ation limit. However these authors have not 
analyzed in detail the fine structure of the bands 
in question. Also there are violent perturbations 
in the rotational levels of the states v’=8 and 9 
(see Olsson?) so that the above-mentioned 
method cannot be applied safely, at least not 
without a detailed analysis of the rotational 
structure. We consider the arguments for case Ic 
given above as sufficiently strong to rule out 
definitely case Ib. Therefore the energy of the 


Fic. 3 (left). Potential curves for case Ib of predissociation 
predissociation limit = dissociation limit). 
Fic. 4 (right). Potential curves for case Ic of predissocia- ° 
tion (predissociation limit >dissociation limit). 


6 G, Herzberg, Ann. d. Physik 15, 677 (1932); G. Biitten- 
bender and G, Herzberg, Ann. d. Physik 21, 577 (1935). 
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predissociation limit E(v’ = 10) = 4.4 ev represents 
only an upper limiting value for the heat of 
dissociation of the Se molecule. 

Lochte-Holtgreven’ has obtained the absorp- 
tion spectrum of S2 in CS, vapor at a temperature 
sufficient to decompose the latter. He claims 
that in this way it is possible to obtain lower 
vapor pressures of S. with less admixture of 
polymers than by the usual method. From the 
spectra obtained he comes to the conclusion that 
at low vapor pressure the first point of diffuseness 
in the Se spectrum (starting with the 11—0 band) 
disappears altogether, and therefore is not due to 
a genuine predissociation.2 However, in our 
experiments even at the lowest pressures used the 
lines of the 11—0 band were still broadened 
compared to those of the preceding and some of 
the following bands and the anomalous intensity 
distribution was still present. Since our ab- 
sorbing layer was 600 cm compared to Lochte- 
Holtgreven’s 65 cm his Se pressure was at least 9 
times larger than that used by us for the same 
intensity of absorption. From the temperature of 
the side tube of our absorption tube (100—180°C) 
the vapor pressure is found to be between 0.01 
and 1 mm. Most of this vapor at a temperature of 
950°C of the main absorption tube is in the 
diatomic form. Furthermore we have carried out 
experiments also with a short tube (20 cm long) 
at.a correspondingly higher pressure and have 
not found any noticeable difference in the 
broadening of the lines or, which is a more 
critical test, in the intensity distribution of the 
bands. 

Therefore, on the basis of our experiments, 
covering a much wider range of pressures than 
Lochte-Holtgreven’s, we conclude that the first 
point of diffuseness in the S2 absorption spectrum 
represents a genuine predissociation which supplies 
a definite upper limiting value for the dissociation 
energy of S2, namely 4.4 ev. Lochte-Holtgreven’s 
value D(S2)=4.7 volts therefore cannot be 
maintained.®* 


The fact that Lochte-Holtgreven arrived at conclusions 
contrary to ours may be explained in the following way: 
7 W. Lochte-Holtgreven, Zeits. f. Physik 103, 395 (1936). 


_ See also the controversy between Lochte-Holtgreven 
and Olsson in Zeits. f. Physik 109, 147 (1938). 


8* Compare also M. Burton and G. K. Rollefson [J. 
Chem. Phys. 6, 416 (1938)] who accepted Lochte-Holt- 
greven’s conclusion as correct in their discussion of 
photochemical primary processes. 


HERZBERG AND L. G. MUNDIE 


Apparently he has not himself taken spectra at consider- 
ably higher pressures with correspondingly thinner layers 
and seems to have been under the impression that at high 
pressures the bands 11—0 etc. are completely diffuse, 
which is not the case, the broadening of the fine structure 
lines even at high pressure being slight when observed with 
moderate dispersion. Therefore when he found only 
very slight or no broadening at the relatively low pressures 
used by him he considered the diffuseness found by others 
at higher pressures to be a case of pressure broadening. 
In support of our explanation of Lochte-Holtgreven’s 
results is the fact that his published spectrogram looks 
exactly like many of those obtained by us under various 
conditions when the intensity of absorption was the same 
as his. Unfortunately he does not reproduce the 11—0 and 
12—0 bands which show the maximum broadening whereas 
the bands 15—0, 16—0, 17—0 which he does reproduce are 
rather sharp anyway.® Also he uses a rather high intensity 
of absorption such that only a few lines of high rotational 
quantum number really come out on his spectrogram. 
The little that is left of the continuous background in such 
a spectrogram is more likely to appear sharp than when a 
moderate intensity of absorption is used. 


In contrast to the first predissociation of Se the 
second one seems to set in with maximum 
probability, occurring at the 18—0 band which is 
entirely diffuse. However on closer examination 
of the spectrogram it is seen that the lines of the 
17—0 band are not as sharp as those of the 
preceding bands.° This shows that the second 
predissociation, too, belongs to case Ic.'!° How- 
ever, since for this second predissociation the 
maximum probability is at any rate not farther 
away from the limit than in the case of the first 
predissociation the two states causing the two 
predissociations cannot lead to the same combi- 
nation of atoms, because, if they did, the potential 
hill for the second predissociation would be 
narrower than for the first and the second 
predissociation should set in more slowly. The 
simplest assumption then is that the first predis- 
sociation leads to normal atoms (#P+*P) whereas 
the second leads to one normal S atom (?P) and 
one Satom in the first excited state ('D, excitation 
energy A=1.1 ev). From the latter assumption 
it would follow (since, as shown, even the second 


8b It should be noted that Lochte-Holtgreven’s number- 
ing of the bands is not the same as the one adopted here 
and by most investigators of the S, spectrum. 

Cf. also the reproduced by Lochte-Holt- 
greven, reference 7. 

10The only alternative would be to assume that the 
broadening in the 17—0 band is due to a predissociation 
— than that in the 18—0 band and the subsequent 
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predissociation limit is somewhat above the 
corresponding dissociation limit) that the dis- 
sociation energy of Se into normal atoms is some- 
what smaller than E(v'=17)—A(!D)=3.6 ev. In 
other words the potential hill for the first 
predissociation would be at least 4.4—3.6=0.8 ev 
high. In view of this height, it is plausible that the 
hill is sufficiently narrow to allow a noticeable 
quantum-mechanical leaking through even for 
such a heavy molecule as So. 

Since the upper limiting value 3.6 ev just 
derived for Do(Sz) is based on an assumption 
about the products of dissociation at the second 
predissociation of S» it is not quite as certain as 
the upper limiting value 4.4 ev derived from the 
first predissociation. However, we consider it 
very likely that the smaller upper limiting value, 
3.6 ev, is also correct. Both upper limiting values 
are in agreement with the assumption that the 
pressure dependent broadening of the absorption 
bands 5—0 to 9—0 found by Kondratjew and 
Olsson" is a genuine induced predissociation and 
not just a pressure broadening since the energy 
of the v’=5 level is 4.2 ev above the ground state. 


(b) NO 

A case similar to that of S: seems to occur for 
NO. In recent compilations of band spectra four 
absorption systems of NO are given, the 8, y, 6 
and ¢ systems. An examination of the spectro- 
grams published by Leifson™ shows that the first 


TABLE I, Wave numbers of the y and ¢ absorption bands of 
NO according to Leifson.? 


v(cm~) 
0-0 44,156 
1-0 46,492 
2-0 48,809 
3-0 51,088 
53,271 
5-0" 55,568 
6—0* 57,813 
60,000 


Av =AG’ 


2336 
2317 
2279 
2183 
2297 
2245 
2187 


: ornnia sumbering is on the assumption that the e bands belong to the 


- ¥ Kondratjew and E. Olsson, Zeits. f. Physik 99, 671 


®S. W. Leifson, Astrophys. J. 63, 73 (1926). 
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four absorption bands of the y system have 
roughly equal intensities. Hitherto the succeeding 
bands were thought to be hidden by the e bands. 
These ¢ bands are very much stronger than the 
preceding y bands but they lie exactly in the 
positions expected for the succeeding y bands. 
This leads us to the conclusion, on the basis of 
the above discussion, that they simply are the 
4—0, 5—0---y bands which have anomalously 
large intensity in consequence of predissociation. 
Thus a predissociation limit is fixed between 51,088 
and 53,271 

In Table I the wave numbers of the bands of 
the y and e systems in absorption are given 
according to the measurements of Leifson. It is 
seen that the differences between successive 
bands are quite compatible with the assumption 
that they form one system. A slight unevenness 
in the differences may either be due to a pertur- 
bation which is likely to occur near the point of 
intersection of two potential curves or to the fact 
that the two groups of bands have so widely 
different intensities and were measured on 
spectrograms of very small dispersion. 

The dissociation energy of NO is known to be 
5.29 volts. The next dissociation limit above this 
is 7.26 volts which is higher than the above 
determined predissociation limit 6.6 volts. There- 
fore we have a predissociation into normal atoms 
and the point of intersection of the potential 
curves of the 2 upper state of the y bands and 
the state causing the predissociation lies high 
above the asymptote of the latter, that is we have 
case Ic of predissociation. This is not surprising if 
we consider that the 72 upper state of the y bands 
has an r, value appreciably smaller and an w, 
value appreciably larger than that of the ground 
state, so that the point of intersection of its 
potential curve with any repulsive curve is 
bound to be high above the asymptote. The state 
causing the predissociation is most probably the 
*> state arising from normal atoms N(#S)+O(°P). 
Apart from this 22 state and the “II ground state 
only quartet states, which are not likely to cause 
such a strong predissociation, result from normal 
atoms. According to the noncrossing rule the 
potential curves of the ?2 state just mentioned 
and of the *2 upper state of the y bands must 
avoid each other. But the fact that only a 
predissociation of the stable = state and no 
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violent perturbations take place indicates that 
the interaction of the two states is slight and that 
the two potential curves almost cross. 

Indications of other predissociations of the NO 
molecule have previously been found. Kaplan™ 
attributed the absence in emission of 8B bands 
with v’>4 to a predissociation. Flory and 
Johnston found a photochemical decomposition 
of NO by light of the wave-length of the 1—0 and 
2—0 6 bands, whereas they noted that 6 bands 
with v’=0 are the only ones observed in emission. 
This predissociation limit in the 6 bands is at 
about the same height as that in the y bands 
found above. It is quite likely that all three 
observed predissociations of NO are caused by 
the same repulsive *2 state resulting from normal 
atoms. It should be pointed out that Flory and 
Johnston’s experiments do not exclude the possi- 
bility of a photochemical decomposition by ab- 
sorption of light of the 4—0, 5—0 y bands which 
should occur according to the above interpreta- 
tion of the anomalous intensity distribution. The 
observation of y bands in emission with v’ >3 (see 
Guillery"), which would contradict our explana- 
tion, does not appear to be definite. 


B. DETECTION OF PREDISSOCIATION IN CuH By 
ANOMALOUS INTENSITY DISTRIBUTION 
IN EMISSION 


Schiiler, Gollnow and Haber" in a recent paper 
have found a remarkable intensity distribution in 
the CuH band at 4280 when excited in a hollow 
cathode at low pressure (0.02 mm). Instead of the 
usual maximum of intensity in each of the 
branches they find that the first line of the P 
branch P(1) is the most intense while all sub- 
sequent lines of the P branch and all lines of the 
R branch are much weaker and increasingly so 
with increasing J. The line P(1) is the only line 
having J’=0 in the upper state. It is thus this 
level J’ =0 which under the conditions mentioned 
is much more strongly populated than the others. 
Schiiler and his co-workers have suggested that 
this anomalous distribution is due to the process 
of formation of the CuH molecules. However it 

18 J. Kaplan, Phys. Rev. 37, 1406 (1931). 

“4 P. J. Flory and H. L. Johnston, J. Am. Chem. Soc. 
57, 2641 (1935). 

18M. Guillery, Zeits. f. Physik 42, 121 (1927). 


16H. Schiiler, H. Gollnow and H. Haber, Zeits. f. Physik 
111, 484 (1939). 
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Fic. 5. Predissociation of a stable 'Z+ state into an 
unstable *>*+ and an unstable 4II state. The level J=0 of 
the 12+ state cannot predissociate, according to Kronig’s 
selection rules, even if such a predissociation were ener- 
getically possible. 


seems impossible to picture a formation of CuH 
in a three-body collision at room temperature 
leading to CuH molecules practically only in the 
J’=0 level which would correspond to an 
extremely low temperature. Also the molecules 
would have to be formed directly in the upper 
state of the bands because, if first formed in the 
lower ground state, in a subsequent excitation 
about as many molecules with J’=1 should be 
formed as with J’=0. 

We would like to propose an entirely different 
explanation for this phenomenon: Numerous 
cases are known of a breaking off of the fine 
structure of electronic bands at a certain J’ value 
produced by predissociation. If above a certain 
energy (J’ value) in the upper state of a band 
system predissociation is possible and if the 
radiationless transition probability is larger than 
the probability of transitions to the lower state 
with radiation a breaking off or at least a sudden 
decrease in intensity of the lines in a band will 
occur. But this is just the phenomenon occurring 
for CuH except that here the sudden decrease of 
intensity occurs for J’=1 so that only one line in 
the band is strong. Therefore we propose to 
explain the anomalous intensity observed by 
Schiiler and co-workers as a predissociation. This 
explanation would of course imply that the 
dissociation energy of CuH is smaller than the 
excitation energy of the J’=1 level of the upper 
electronic state, that is <23,325 cm= (or <2.892 
ev). This is entirely possible since a linear 
extrapolation of the vibrational quanta of the 
ground state yields Do(CuH) ~ 24,300 and 
since such linear extrapolations are known to 
give Do values which are usually somewhat too 
large. 

It is however not necessary to make the 
assumption that the predissociation limit just 
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happens to fall between the levels J’=0 and 
J’ =1 of the upper 'Z+ state of the CuH band as 
long as it is below 23,325 cm~, because if the 
unstable state causing the predissociation is « Il or 
3>+ state the J=0 level of the '=* upper state of the 
CuH band cannot predissociate no matter where the 
predissociation limit lies. This is illustrated by 
Fig. 5. According to Kronig’s selection rules!’ 
only states with the same value of J and the 
same symmetry can predissociate into each 
other. But a 'II state has no level with J=0 and 
in a *2* state the only level with J=0 has the 
opposite symmetry (—) to that of the J=0 level 
(+) of 12+ (see Fig. 5). 

A Cu atom in the 2S ground state and an H 
atom in the *S ground state give rise to two 
molecular states and of which the former 
is the ground state. Therefore if a predissociation 
into normal atoms is assumed the only state that 
can cause the predissociation is the *2* state and 
such a state according to the above cannot cause 
a predissociation of the J’=0 level, in agreement 
with observation. Thus under this assumption 
the predissociation limit and, a fortiori, the 
corresponding dissociation limit may be ap- 
preciably below the energy of J’=1 (23,325 cm~?). 

It is not impossible although much less likely 
that the predissociation leads to a normal H 
atom and a Cu atom in the excited 7D state. In 
that case six states arise from the separated 
atoms. For two of them, *2+ and 'II the same 
arguments as above would apply for the non- 
predissociation of the J’=0 level (see Fig. 5). 
But the predissociation may also be caused by 
two of the other four states arising from 2S+2D 
(namely ‘2+ or *II) and then the nonpredissoci- 
ation of the level J’=0 would not follow auto- 
matically. Moreover in this case the heat of 
dissociation of CuH into normal atoms would 
have to be <12,125 cm™ (i.e., <1.503 ev) which 
appears to be too small although not impossible. 
Thus we consider it as most likely that the 
bredissociation is caused by the state re- 
sulting from normal atoms and that therefore 
Do(CuH) < 23,325 

Our assumption that the upper ‘2+ state of the 
CuH band predissociates into *E+ contradicts 


“See G. Herzberg, Molecular Spectra and Molecular 


= I. Diatomic Molecules (Prentice-Hall, New York 
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Kronig’s selection rule that the spin S should not 
change in a radiationless transition. But inter- 
combinations between singlet and triplets have 
been found to occur for predissociations of the 
heavier molecules just as for transitions with 
radiation.” 

In a predissociation '2+—*=*+ as well as in a 
predissociation the radiationless transi- 
tion probability increases roughly proportion- 
ately to J. While for '2*+—'Il this has been known 
for a long time (Kronig) it follows for '2+—*r* 
from a recent paper of Budé and Kovacs'* who 
calculated the perturbations between a '* and 
a *=*+ state. The observations of Schiiler and 
co-workers are in qualitative agreement with this 
theoretical conclusion since they found that the 
intensity continues to drop even after the jump 
from J’=0 to J’=1 whereas if the radiationless 
transition probability were independent of J one 
would expect to find after the breaking off point 
an increase of intensity up to the usual thermal 
maximum which would occur at P(5). 

There remains to explain the fact that at 
higher pressures the complete CuH band appears 
with the usual intensity distribution, the line 
P(i) having no longer any abnormally large 
intensity. It seems to us that this is the well- 
known pressure effect of predissociation also ob- 
served for AIH, CaH and SrH.'’:'* At higher 
pressures in consequence of the collisions the 
molecules in the J’=0 state, whose natural life is 
longer than in the states with J’=1, will go back 
to the ground state without radiation or will go 
over to the other rotational levels of the excited 
state, or will be dissociated by collisions (in- 
duced predissociation), and eventually a distri- 
bution corresponding to thermal equilibrium is 
set up. While this is the same explanation as for 
the other observed cases of pressure effect, the 
remarkable thing in the case of CuH is that such 
very low pressures are necessary in order to 
bring out the anomalous intensity distribution 
produced by predissociation. Even at a pressure 
as low as 2 mm the P(1) line is not the most 
intense line of the band. This occurs only at a 
considerably lower pressure. However, this can 
easily be understood because here (unlike the 


ass “y Budé and I. Kovacs, Zeits. f. Physik 111, 633 
19 See also Olsson, Dissertation, Stockholm, 1938. 
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Fic. 6, Observed rotational levels in the 4II state of AID 
according to Holst and Hulthén, reference 21. The broken 
a refer to predissociating levels (see, however, reference 


other cases) we have an intercombination for 
which the radiationless transition probability is 
not very much larger than the probability for 
transitions with radiation, as is also evidenced by 
the fact that lines with J’>0 do appear at the 
lowest pressure used even though with much 
smaller intensity than the line with J’ =0[P(1) ]. 
Therefore it requires only comparatively few 
collisions to produce a nearly thermal distribu- 
tion. To be sure the distribution is not yet 
entirely thermal at a pressure of 2 mm. 

Very probably it will be possible to detect 
weak predissociations of other molecules by 
observing their band spectra in emission at such 
low pressures as were used by Schiiler and co- 
workers for CuH. 

It should be remarked that according to the 
above interpretation one would expect to find a 
similar intensity distribution in CuH bands of 
the same system with v’>0, but such bands 
have not been investigated by Schiiler and his 
co-workers. 


C. ON PoTENTIAL HILLS IN THE PREDISSOCIATION 
or AIH anp BH 


If one wants to determine the dissociation 
energy of a molecule from predissociation data it 
is necessary to know whether the predissociation 
occurs across a potential hill or not. The criterion 
used in Section A for deciding this question in the 


case of the predissociation of the Se molecule 
(noncoincidence of predissociation limit with 
maximum probability of radiationless transition) 
is in general not applicable to predissociations 
observed by a breaking off of emission bands. 
But as pointed out by one of us” a decision is 
possible if a breaking off of the rotational fine 
structure can be observed in at least two but prefer- 
ably three successive vibrational levels. 

In Fig. 6 the observed rotational energy levels 
in the 'II state of AID are represented according 
to the data of Holst and Hulthén.”" The pre- 
dissociating levels are indicated by broken lines. 
The fact that predissociation starts at a higher 
energy for the lower vibrational levels is known 
to be due to the fact that the heights of the 
maxima of the effective potential curves increase 
with increasing rotational quantum number J. 

The effective potential energy Uys(r) for the 
rotating molecule (in cm~ units) is given by 


Us(r) = Uolr) + J(J+1), (1) 


where Up(r) is the potential energy without 
rotation (J=0). If Uo(r) has no maximum the 
maxima of U,(r) will shift to smaller r values 
with increasing rotation. If the energies of the 
maxima Umaz(J) are plotted against J(J+1) a 


1 
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Fic. 7. Breaking off points of AID and AIH for the 
various vibrational levels of the II state plotted as a func- 
tion of J(J+1). The circles represent the last nonpredis- 
sociating levels, the squares the first predissociating levels 
for the respective vibrational levels. The ordinate gives 
the energy above the minimum of the ground state 12. 


20 G. Herzberg, Ann. d. Physik 15, 677 (1932); J. Phys. 
Chem. 41, 299 (1937); Astrophys. J. 89, 290 (1939); see 
also reference 17. 

1930)" Holst and E. Hulthén, Zeits. f. Physik 90, 712 
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curve that is convex to the abscissa axis is 
obtained (limiting curve of dissociation). It 
intersects the ordinate axis with horizontal 
tangent at an energy value equal to U)(), ie., 
equal to the dissociation limit of the electronic 
state in question.” Jf, however, Uv(r) has a 
maximum the maxima of the effective potential 
curves U(r) will all be at approximately the same 
r value and therefore according to (1) approxi- 
mately a straight line will result for Umax(J) if 
plotted as a function of J(J+1). This will intersect 
the ordinate axis at Uo(rmax), the potential 
energy at the point of the maximum (without 
rotation). This conclusion is independent of 
whether the potential maximum for zero rotation 
is due to an intersection of two potential curves 
or not. 

From the experimental breaking off points a 
few points on the curve Umax(J) as a function of 
J(J+1) may be obtained since the breaking off 
will occur for those levels for which the sum of 
vibrational and rotational energy equals Umax(J) 
for the particular J value. In Fig. 7 the breaking 
off points in the 'II state of AlD, represented in 
Fig. 6, are plotted as a function of J(J+1) and 
thus the Umax curve is obtained. It is seen that 
the points lie very closely on a straight line which 
proves that the potential energy without rotation 
has a maximum.”’ Fig. 8 shows the corresponding 
plot for the 'II state of BH (calculated from the 
data of Almy and Horsfall*). Again a straight 
line is obtained proving that the potential curve 
without rotation has a maximum. This is not 
surprising since the two 'II states of AID and BH 
are analogous also in other respects. 

It is only in rare cases that a breaking off is 
observed (as for AID and BH) in three successive 
vibrational levels of an electronic state of a 

* See G. Biittenbender and G. Herzberg, Ann. d. Physik 
21, 577 (1935); G. Herzberg, J. Phys. Chem. 41, 299 
R. Schmid and L. Gerd, Zeits. f. Physik 104, 724 

* Olsson, reference 19, gives slightly lower breaking off 
points obtained from unpublished observations at low 
pressure. If his data are plotted they do not fall as exactly 
on a straight line as the points in Fig. 7 (taken from 
Holst and Hulthén, reference 21). However, the curve 
obtained does not have a horizontal tangent for J(J+1)=0 
and therefore the above conclusion as to the existence of a 
Fetentinl hill would not be changed by adopting his data. 

e curvature if real would be due to the fact that for 
large J values the position of Umax(J) shifts slightly even 
when there is a potential hill in the rotationless curve. 


(1935) M. Almy and R. B. Horsfall, Phys. Rev. 51, 491 
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Fic. 8. Breaking off points of BH in the 'II state. Almy 
and Horsfall whose data have been used here give only the 
‘J value for the sudden decrease of intensity but do not 
state where the decrease starts. Therefore only one point 
is given for each vibrational level. The ordinate gives the 
energy above the level v=0, J=0 of the ground state. 


diatomic molecule. If such a breaking off is only 
observed in two successive vibrational levels a 
decision as to whether the rotationless curve has 
a maximum or not may yet be obtained. The 
slope of the straight line through the two breaking 
off points in a diagram such as Figs. 7 and 8 
according to (1) gives the mean 7 value of the 
two corresponding maxima of the effective po-— 
tential curves. This 7 value is equal to the r value 
of the maximum of the rotationless curve if 
such is present. However, the maximum of a 
rotationless curve (or, in other words, the 
intersection of two potential curves) always lies 
at comparatively small r values, usually smaller 
than twice the equilibrium internuclear distance 
r. of the state considered, whereas the maxima of 
the effective potential curves for an ordinary 
rotationless curve (without maximum) lie at 
appreciably larger r values, at least for not too 
large J values. 

Thus a large slope of the line connecting two 
breaking off points® in a diagram like Fig. 7 
indicates a maximum of the rotationless potential 
curve or in case J of predissociation (radiationless 
electronic transition) it indicates the subcase (c) 
in which the point of intersection of the two 
potential curves lies above the asymptote of one 
of them (see Fig. 4). On the other hand, a small 
slope of the line in question indicates a normal 


* Or a large initial slope of the curve obtained from more 
than two breaking off points (see reference 23). 
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potential curve or in case I the subcase (b) in 
which the point of intersection is below the 
asymptote (see Fig. 3). 

As an example let us consider the ‘II state of 
AIH which (unlike AlD) shows only two breaking 
off points. These also are indicated in Fig. 7. The 
slope of the line gives fmax=2.643 whereas 
r-=1.646 thus indicating a maximum of the 
potential curve Uo(r).2* Of course this was to be 
expected since the maximum has already been 
proven for AID. The 7 value for the maximum 
comes out to be the same for AIH and AID 
within less than 0.01A. 

As shown above, the intersection of the straight 
line in Fig. 7 or 8 with the ordinate axis gives the 
energy of the potential maximum of the rotation- 
less curve (apart from the influence of the tunnel 
effect considered below). For AID this energy 
comes out to be 25,173 cm for BH 28,350 cm= 
above the v=0, J=0 level of the ground state. 
These values represent upper bounds for the 
dissociation energies of the respective molecules, 
but since there is a potential hill the actual values 
will be smaller by an appreciable amount. 

If the energies are referred to the minimum of 
the potential curve of the ground state, as is done 
in Fig. 7, the points of intersection of the straight 
lines for AIH and AID with the ordinate axis, that 
is the maxima of the rotationless potential curves 
of the two isotopic molecules, do not coincide 
exactly as one might perhaps expect. The values 
are 25,610 cm and 25,775 cm™, respectively. 
The difference, as has previously been pointed 
out by Olsson,?’ is due to the quantum-mechanical 
penetration through potential barriers (tunnel 
effect) near the top of the hill. In consequence of 
this effect the breaking off points in Fig. 7 and 8 
do not actually give the maxima of the effective 
potential curves but energies somewhat below. 
The difference will be larger for AIH than for 


26 As early as 1933 E. Hulthén and R. Rydberg [Nature 
131, 470 (1933) ] suggested that there is a potential hill in 
the rotationless curve of the 'II state of A1H, but did not 


ive a satisfactory proof for it. In a recent paper [Zeits. f. 

hysik 113, 126 (1939) ] Hulthén, in referring to this 
—— hill, says that it is so narrow and low that the 

reaking off at low pressure is not influenced by it and 
that therefore the breaking off point in the level v’=2 of 
AID gives directly the dissociation limit of the 4II state. 
In our opinion this is by no means the case. Hulthén’s 
discussion refers only to the very top of the potential hill 
that has been established here. 

27 E, Olsson, Zeits. f. Physik 108, 322 (1938). 
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AID because H penetrates the potential hill 
more easily, that is, farther below the top, than D 
and this explains the small difference in the 
energies of the breaking off points for zero 
rotation of AIH compared to AID. While of 
course the difference in energy between the top of 
the potential hill and the breaking off point for 
AIH is larger than the difference between the 
breaking off points for AIH and AID it will 
certainly be of the same order of magnitude, that 
is the influence of the tunnel effect is slight. We are 
therefore justified in applying the above semi- 
classical method to the detection of the potential 
hills. The influence of the tunnel effect will be 
smaller yet for atoms heavier than H or D, and in 
general may be neglected for them except for very 
narrow potential hills, such as the one for S» 
discussed above where the influence is certainly 
present but yet slight. 

The dissociation limit has of course the same 
energy above the minimum of the potential curve 
of the ground state in AIH and AID and conse- 
quently the smaller of the two upper bounds (that 
for AIH) applies to both isotopic molecules. 
Subtracting the zero point energies we obtain for 
the dissociation energies referred to v=0, J=0 
of the ground state Do(ALH) <24,775 cm— (=3.072 
ev) and D,(AID) <25,010 cm (=3.101 ev).*8 

The above considerations concerning potential 
hills and the method of their detection, as has 
been mentioned before, are equally applicable for 
case Ic and III of predissociation® i.e., they are 
independent of whether the potential hill is 
formed by the intersection of two potential 
curves or whether it is present in a single zero 
rotation potential curve. According to von 
Neumann and Wigner if the potential curves of 
two electronic states of the same species cross in 
a zero approximation they will not do so in 
higher approximation. Thus, strictly speaking, 
case Ic is not possible for two states of the same 
species. But a maximum of a single potential 
curve may always be thought of as brought about 
by the intersection of two zero approximation 
curves of the same species (an attractive and a 
repulsive one), and it would depend on the 


*8If Olsson’s data, reference 23, for the breaking off 
= were used very slightly smaller limits for Do would 
obtained. 
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strength of the interaction of the two states 
whether one would consider the predissociation 
as case Ic or III. 

In the cases of AIH (AID) and BH we have 
apparently a case of two states of the same 
species interacting so strongly that these predis- 
sociations are most adequately described by case 
IIT. If the two states were not of the same species, 
or if for like species the interaction were small, 
and case Jc would apply one would expect to 
find, at least in absorption, transitions to higher 
vibrational levels beyond the predissociation 
limit, which at least for AlH is definitely not the 


case. The 'II state that predissociates may result 
from normal atoms Al (or B) ?P)+H (2S). But 
such a state also arises from the first excited ?P 
state of Al (or B)+H (2S). Apparently the former 
II state is “originally,” that is in zero approxi- 
mation, repulsive whereas the upper 'II state is 
attractive. The latter for smaller r presses the 
former downward giving rise to a minimum in its 
potential curve and a maximum at larger r 
values. Information about the height of the 
maximum can be obtained only after the heat of 
dissociation of these molecules has been obtained 
from other data. 
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The Partition Functions and Energy Levels of Molecules with Internal 
Torsional Motions 


Bryce L. CRAWForD, JR.* 
Mallinckrodt Chemistry Laboratory, Harvard University, Cambridge, Massachusetts 
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The complete energy matrix is set up for external and (hindered) internal rotations of 
“pseudo-rigid” molecules; i.e., molecules consisting of a rigid framework to which symmetric 
hindered rotators are attached. The reduced secular equation for the energy levels of such a 
molecule with one internal rotation is given explicitly. By using the perturbation method of 
Wilson, the partition function for this model is found, correct to the second order. Effects due to 
angular momentum couplings are included. A simple method for the calculation of thermo- 
dynamic contributions of hindered internal rotations is given. As an example, the application of 
the method to propylene is described; in this case, the coupling terms here included may amount 
to as much as 19 percent of the effect of the hindering potential on the thermodynamic 


contributions. 


[' is now generally accepted that the internal 
rotations of methyl groups in organic com- 
pounds are restricted by potential barriers of 
appreciable magnitude. In order to determine 
these barrier heights from thermal, spectroscopic, 
and molecular structure data, it is necessary to 
study the dynamics of the molecules and to 
evaluate the partition functions for both the 
external and internal rotational degrees of 
freedom. 

The quantum mechanics has been worked out 
by Nielsen! for the molecule consisting of two 

* Present address: Sterling Chemistry Laboratory, Yale 
University, New Haven, Connecticut. This research was 


carried out while the writer held a National Research 
Fellowship in Chemistry. 


H. H. Nielsen, Phys. Rev. 40, 445 (1932). See also J. B. 
Howard, J. Chem. Phys. 5, 451 (1937). 


tops on a common axis. If these tops are identical, 
the wave-functions are periodic Mathieu func- 
tions; the thermodynamic contributions of such 
systems have been evaluated by Pitzer.? Ethane, 
C2D¢,* and dimethyl acetylene‘ correspond to 
this model. 

The quantum mechanics has not been worked 
out, however, for molecules in which the axis of 
the top does not pass through the center of 
gravity (e.g., methyl alcohol), or for molecules 
with several tops (e.g., propane). In this paper 
we present a study of the ‘‘pseudo-rigid’”’ mole- 

2K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 

3J. D. Kemp and K. S. Pitzer, J. Am. Chem. Soc. 59, 
276 (1937). G. B. Kistiakowsky, J. R. Lacher, and F. Stitt, 
J. Chem. Phys. 7 289 (1939). 


4B. L. Crawford, Jr., and W. W. Rice, J. Chem. Phys. 7, 
437 (1939). 
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cule, the over-all moments of inertia of which 
do not change with internal rotations. Such 
molecules consist of a rigid framework to which 
symmetric tops are attached. Examples of this 
class are propane, isobutene, acetone, and the 
methyl amines; simple molecules falling outside 
this class include ethyl alcohol and n-butane. 


THERMODYNAMIC RESULTS 


In this section we shall describe the nature of 
the partition function obtained, and shall give 
the final expressions for the thermodynamic 
contributions, postponing the rather involved 
derivation. 

The partition function for the external and 
internal rotations may be written 


Q=Qn I Qu(m), (1) 


where Qz is the partition function for over-all 
rotation, and Qy(m) is the partition function 
for the internal hindered rotation of the mth 
top. The effects of angular momentum couplings, 
etc., are all included in the Qu(m). Hence the 
thermodynamic contributions of the over-all 
rotations are the usual ones for a rigid molecule; 
the moments of inertia, J;, Is, I3, are to be 
evaluated for the molecule as a whole. The 
symmetry number for these contributions is 
obtained by considering the internal rotations to 
be frozen; e.g., o=6 for ethane, o=1 for 
propylene and for methyl alcohol. 

The contributions of the internal hindered 
rotations being additive, we may describe Qu(m) 
for a single such motion. Let the mth top have a 
moment of inertia A» about its axis, and let 
Xmiy Am2, Ams be the direction cosines of the top 
axis with respect to the principal axes of the 
molecule as a whole. For the hindering potential, 
we use (Vm/2)(1—COS ¥mam), Where am is the 
angle of rotation of the top about its axis, 
Ym is the number of minima encountered in one 
complete rotation, and V,, is the height of the 
maxima. 

_ The energy levels for the internal motion 
depend on the eigenvalues a,, of the Mathieu 
equation 


(d*p/dx*) + (dmr+24m cos 2x)y=0. 


These dm, depend upon the quantum number 


BRYCE L. CRAWFORD, JR. 


r (r=0, +1, +2, etc.), and upon the quantity 6, 
which is essentially the hindering potential : 


Om = Vin(827A mQm/h?¥m?). (2) 


Here h is Planck’s constant and Q, is an im- 
portant geometrical parameter: 


(3) 


The nature of ©, is best understood from the 
fact that (A,,Q,) may be regarded as a “reduced 
moment” for the top. Tables of the necessary 
Amr eXist. 

The internal energy levels are vibrational in 
character for low values of the quantum number 
r, and rotational for high 7; the angular mo- 
mentum interactions for these two types of 
level are different. The partition function Qy(m) 
is found by summing over the energy levels in 
the usual way; it turns out that the interactions 
can be included if we break the sum into two 
parts, one for the “vibrational” and one for the 
“rotational” levels, and use appropriate factors 
for the two. The point of separation, r=7ro, is 
easily found on inspection of the a», for the 
particular case. 

Writing 


Vine = m&mRT (4) 

where R is the gas constant, we find 
Qu (m) =exp (—gm) (5) 
Gr= exp (—QnY mr), 


exp (— Yn,). 


(6) 


The index R indicates that the sum is over all 
r2Yro, and the index V indicates summation over 
all r <r. The factor exp (— qm) merely adjusts the 


zero of energy: 
Qm= Vm/2RT. (7) 


The thermodynamic contributions due to the 
mth top may be written: 


(E/T)m=R[g+(H/G)], 
(F/T)u=R[g—In G], 
Su G+(H/G)], 
Cu=RL(J/G) 


(8) 


= 
| 
| 
T 
be 
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where 
G=(QniGret+Gy); H=(QniHr+Hy); 
Ymr EXP (— Yr), 
Tr= Qn? exp (—QnYmr), 


R 


exp (— Yur). 
Vv 


(9) 


It will be seen that for cases in which 2, may 
be taken as unity, the distinction between Gp 
and Gy vanishes, and the partition function (5) 
reduces to the simple sum over states for a 
Mathieu system. The corresponding thermo- 
dynamic contributions have been given by 
Pitzer? in his Eqs. (2) and (3) and tables, so 
that no direct calculation is necessary in such 
cases. The ‘“‘reduced moment,” A»,Q2n, becomes 
simply the moment of the top. 

The ‘quantum effects,” which require the 
introduction of the symmetry number in Qp, are 
particularly easy to take into account in Qu(m). 
It turns out that the result obtained by using 
only levels of a certain symmetry, say of class A, 
in our sums Gr, Gy, etc., gives a very good 
approximation, except at very low temperatures. 
The error may be estimated by evaluating 
Gy(A), using A levels only, and Gy(B), using B 
levels only, and noting the difference, ro being 
taken so that Gy(A) and Gy(B) have the same 
number of levels. No symmetry number factor 
is used. 

These A and B eigenvalues are those which are 
given in the ordinary tables such as those of 
Ince ;5 with his notation, the eigenvalues a, and 
b, belong to the A class if is even, and to the 
B class if n is odd. In evaluating Gr, Hr, Jr, the 
following approximate expression is useful for 
finding higher values of amr; in Ince’s notation, 


]. 


This formula is valid for small @ and large n; 
for n=6, 9=16, the error in a, is less than 1 


*E. L. Ince, Proc. Roy. Soc. Edinburgh 52, 355 (1931), 
especially p. 363. Ince’s @ corresponds to our @m; his dn, 
» correspond to our @mr, for positive and negative r, 
respectively, if m is even, and vice versa if m is odd. 
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percent. A convenient method of finding the 
eigenvalues accurately is given in reference 11. 

For very small values of Q,,, it may be neces- 
sary to use the exact expressions for Gr and Gy 
given in Eq. (42); the error due to use of the 
above expressions (6) may be estimated by noting 
the effect of changing the position of ro to include 
one more or less level in Gy. 

The application of this method to the case of 
propylene is discussed in the last section. We now 
proceed to the derivation of the energy matrix 
and the partition function described above. 


KINETIC ENERGY 


Let there be V symmetric tops in the molecule, 
attached to a rigid framework. Let J;, Iz, I3 be 
the moments of inertia of the molecule as a 
whole (framework and tops) about the principal 
axes. These axes and moments are independent 
of the rotations of the tops. We retain the nota- 
tions described in the previous section. Then it 
is not hard to show® that the classical kinetic 
energy is given by the expression 


N 3 3 ON 

m=1 i=1 i=1 m=1 (10) 
where w; is the instantaneous angular velocity of 
the molecule about the principal axis i. 

This expresses the kinetic energy as a quadratic 
form in the (V+3) velocities, w; and dm. The 
Hamiltonian form of the kinetic energy is a 
quadratic form in the conjugate momenta 


P;=(0T/0w;), Pm=(dT/Odm). 

The matrix of the Hamiltonian quadratic form is 
the inverse of the matrix of (10); carrying out 
this inversion (see below) yields the expression 


1 


™m m nm 


+2> LmiPmP i, (11) 


where the terms have the following meanings. 
The first term gives the energy of over-all rota- 


tion of the molecule; the 8;; are analogous to 


6 L.S. Kassel, J. Chem. Phys. 4, 276 (1936) ; or one ma 
use the method of E. B. Wilson, Jr., and J. B. Howard, if 
Chem. Phys. 4, 260 (1936). 
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reciprocal moments of inertia, and are best 
defined through the matrix 


(gij) = (6i;— Tia), (12) 
where 


and 6;; is the Kronecker symbol. Then 
Bij = (14) 


where (g*’) is the matrix inverse to (g;;). If g is 
the determinant of (g;;), then the determinant 
of B is 


(15) 


The second term in (11) gives the energy of the 
internal motions; the F,, are corrected moments 
of inertia: 


where 


Ann= (17) 


The third term expresses the coupling between 
different tops (note that only terms between 
different tops are included in the summation), 
while the last term gives the angular momentum 
coupling between the top motions and the over- 
all rotations; here 


Since the coefficients in (12) are constants, 
this equation gives the quantum-mechanical 
kinetic energy operator, if for P; we write P;, and 
for Pm, Pm.7? Here the P; are the angular mo- 
mentum operators studied by Klein,® and 
P,, = (h/271)(0/dam). The wave functions must 
be periodic in a, with period 27. 

The inversion of (10) is conveniently carried 
out by matrix technique.® Denoting the matrix of 
the quadratic form (10) by (J), we may break 


7B. Podolsky, Phys. Rev. 32, 812 (1928). This procedure 
has been justified for P,, P,, P. by Wilson and Howard, 
reference 6. 

80. Klein, Zeits. f. Physik 58, 730 (1929). 

*See, for example, R. A. Frazer, W. J. Duncan, and 
A. R. Collar, Elementary Matrices (Cambridge University 
Press, 1938). 


it up into four sub-matrices: 


Here 


(T1)mn=Amébma, (T3)mi=AmAmis 


and 7; is the transpose of 73, these two being 
rectangular matrices. If we break up (7-) in 
the same way, writing 


and use the fact that (7)(7-)=(T-)(T)=E, 
where £ is the identity matrix, we can find 
relationships between the sub-matrices. We 
obtain 


T 
and the problem is reduced to the inversion of the 
3-by-3 matrix (8-'), since 7;-' and T.~! can be 
written down at once. 
Case in which N=1 


For a molecule with a single top, (11) re- 
duces to 


1 
PP Pm +22 is (19) 


The elimination of the sum over m results in 
simple expressions for the parameters: 


pm, 
Fn =AmQn = Pm/Amm; (20) 
Ymi= —Xmi/TiQm. 


For the molecule composed of two tops on a 
common axis, (19) can be transformed into 
Nielsen’s equation,' which involves no _ inter- 
action terms. Unless the tops are identical, 
however, the new boundary conditions demand 
nonperiodic solutions, as Nielsen points out. 
It is more convenient to work with periodic 
solutions, retaining the interaction, therefore, 
except in the case of ethane-like molecules. 


(20) 
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Case in which N>1 


For a molecule containing several tops, it 
proves convenient to regard the 7;; as small, and 
neglect terms of second order in 7;;. In most 
applications of interest, methyl groups will be 
involved ; if there are two or more methyl groups 
attached to a framework, the ratio A,,/J; will 
indeed be small, and it will be reasonable to 
neglect 7;;°. On this basis we find 


Qn =(1— pm), 
Fn=AmQn; 
Ymi= —Xmi/ Ti, 
(AmiAni/ Ti), 


F 
where 


Pm = 


i=1 


ENERGY MATRIX 


Potential energy 


The energy operator consists of the kinetic 
energy (11) plus the potential energy, which we 
write 

ymam), (22) 


neglecting possible interaction terms of the type 
Vinn(Qm,@n). The use of the cosine form is de- 
sirable for mathematical simplicity, and is sup- 
ported (though not conclusively) by the heat 
capacity data on ethane.* The parameters de- 
fining the potential barrier, are vm, the number 
of minima encountered in one complete rotation, 
and V,,, the height of the maxima, as described 
in the first section. 
The complete Hamiltonian may be written: 


H=H,+Hhi, (23) 
where 
Ho=P?+>H.., 
Hi= DP 2" AmnPmP a, 
2H,,= (1/Fn)(Pm)?+ — cos VmO&m) 
P?=(1/2))) 


Pn = LeymiPi. 


Zeroth-order energy 


The term P? gives the rotational energy of a 
top which has the moment of inertia tensor (8~'). 
The rotational energy matrix may be written 
down in terms of the usual quantum numbers 
J, K, 

The term H,, gives the energy of internal rota- 
tion for the mth top. The wave-equation ob- 
tained from the operator H,, may be written in 
the form 


r/ OX? + 26m COS 2x)~mr=O0, (25) 


which is one form of Mathieu’s equation,” with 
2X=VmOm. For the quantum number, 7, we use 
that of the limiting free rotator derived by setting 
6 =0; hence r=0, +1, +2, ---. The degeneracy 
of the “positive” and “‘negative’’ levels is lifted 
for 0,0, if r is a multiple of (v,/2); the 
“‘negative”’ level has the lower energy. In general, 
the energy will be 


= (Amr +20m) (h?¥m?/327? Fm). (26) 


The parameter @,, is given by Eq. (2). Unfortu- 
nately, dm, is not a simple function of r and 8,,. 
As was mentioned above, the lower levels of 
this Mathieu system closely approximate those 
of a harmonic oscillator, while as 7 increases, the 
levels change and become like those of a free ro- 
tator. This transition takes place quite abruptly, 
so that we may speak of “vibrational” and 
“rotational” levels in this connection. 


Energy matrix 


It will be convenient to write the energy matrix 
H in terms of sub-matrices, each sub-matrix 
corresponding tocertain sets of Mathieu quantum 
numbers, 7, and being itself a function of J, K, 
and M. (Throughout we shall use 7 and s as 
Mathieu quantum numbers for the mth and mth 
tops, respectively.) If we write 


10 This form is chosen to conform to Ince, reference 5. 
as his tables of the eigenvalues @m, are the most extensive 
we have found. The discrepancy in the sign of the term in 
cos 2x arises because we choose to regard Vm and @, as 
positive; this difference does not affect our results. 
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then the sub-matrices of H are 


+ { jE, 


(28) 


where E is the identity matrix. This expression 
is valid for any form of potential energy which 
does not include interactions Vinn(a@m,@n). 

The large H matrix can be factored by using 
the symmetry properties of the Mathieu func- 
tions. The levels can be classified into the sym- 
metry classes of the group D>," (e.g., D3 for 
the methyl group in C2H¢). Using the fact that 
P,, has symmetry Ae, we can show that Ym,rr must 
vanish unless the r and 7’ levels belong to one 
of the pairs (A1,A2), (B1,Be), (Ex,Ex). By expand- 
ing the E; functions in exponentials instead of 
sine and cosine terms, we may eliminate inter- 
actions between E;,+ and E;-. Thus we have H 
factored in IInvm factors, of symmetries (A:+Az2), 
(B,+B:), E,*+, for each top. 


Limiting values 


For the ‘“‘rotational’’ levels in each factor, the 
limiting value of dm, is 47?/vm?. If we use expo- 
nential expansions for the ‘‘rotational’’ wave 
functions of the (Ai+Az2) factor, as we may 
without difficulty, the limiting value of Ymrr- is 
the same for all factors; viz., 


Omer? =(h/2m)r (29) 


For the “vibrational” levels, there will be v», 
levels, one of each symmetry factor, with the 
same energy. Writing ¢ for the “vibrational” 
quantum number which we may use in place of 
r, we find 


Amt = 2(2t+1) (Om)! — 24m 


1 E. B. Wilson, Jr., J. Chem. Phys. 6, 740 (1938), gives 
the symmetries for v»,=6, corresponding to C2Hg in his 
coordinates. The general case is in exact analogy. It may 
be pointed out that these symmetry properties do not 
depend on our use of the cosine-form potential energy ; 
they demand only that V,, be an even function of am with 
period (27/vm). 

2 As Om—> ©, or r—>0, a group of vm levels will coalesce 
to give the vm-fold degenerate vibrational level designated 
by the vibrational quantum number t¢. (¢=0, 1, 2, --+.) In 


(30) 
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and the nonvanishing elements v4 are 


Um, t, t+1= —Um, t41, t= (h/4ni) Vm(E+1)!Om?. (31) 


Energy levels 


To find the separate energy levels from the 
energy matrix, we must write out P? and p,, in 
terms of J, K, M. First, however, we shall find 
it convenient to rotate axes so that the matrix 
B becomes diagonal. This rotation also affects 
the pn. If C is the orthogonal transformation. 
(rotation of axes) which diagonalizes 6, so that 


then we can easily show that we must take 
Pn= i's 
mi = 
i 


The insertion of Klein’s matrices* for the P,’ 
yields the energy sub-matrices H”’ and H””’ as 
explicit functions of J, K, and M. The complete 
energy matrix H is diagonal in J and M. By 
regarding the matrices H"”’ as perturbations, 
and applying perturbation methods," the energy 
matrix may be rendered diagonal in r, correct 
to the third order in the elements of H”’. The 
expression for the energy of a molecule with a 
single top, for example, may be written (the 
subscript m has been dropped) 


Ws, x, M, r= (h?/8x") 


(32) 


where 7 is the root of the secular equation whose 


this group there will be one level belonging to each of the 
symmetry factors mentioned above. The relation between 
7 and t may be written 


Vmt = |2r| —k 


with 0<¢ k<vm for r positive, and O<k< vm for r negative. 
The plot of wm, against @, (reference 1) makes this clear. 

13 J. H. Van Vleck, Phys. Rev. 33, 484 (1929); O. M. 
Jordahl, Phys. Rev. 45, 87 (1934); Wilson and Howard, 
reference 6. In applying this method, we have regarded 
the spacing between interacting Mathieu levels 7, 7’ as 
large compared to the over-all rotational spacing; 1.¢. 
W,: — Ww, > B(2K +1). This will be true for the ‘‘vibrational 
levels; for the “rotational” levels, the value of 
decreases to zero and the interaction becomes negligible 
as (w,-—w,) becomes small. 
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nonvanishing elements are 
(K,K) = —9+K?+ 
(K,K+1) = 
(33) 


The matrix is of course Hermitian; the symbols 
have the following meanings: 


= 
X(J+K)(J+K+1)]}, 
$=1, 2, 3, 
By=B;—(B,+Bs)/2. 


The secular equation for the molecule with 
several tops may be found by the same technique. 


PARTITION FUNCTION: SINGLE Top 


We shall use the method given by Wilson" to 
find the partition function from the energy 
matrix (28). It is not necessary to diagonalize . 
We shall first develop the case N=1, and then 
extend our results to the general case. 


Zeroth-order expression 


If we include only the terms of Ho, as given 
in (24), we see at once that the partition function, 
Q, for internal and external rotational motion, is 
given by 


Qu (Gms Om), (34) 


where Qz is the usual rotational partition func- 
tion for an asymmetric top of moment of inertia 
product |6-'| (Eq. (15)). 


Qr = 3g)! (35) 
Qu is the sum-of-states for a Mathieu system, 


* E. B. Wilson, Jr., J. Chem. Phys. 7, 948 (1939). This 
aes is based on the perturbation technique of reference 


as evaluated by Pitzer.’ 


= exp (—tm,/RT) 


=exp (—ge) exp (36) 
gm = Vn/2RT. 


First-order correction 


The interaction term in the term 
is linear in the P;. We may take account of this 
term by a change of variable which eliminates 
linear terms; the transformation 


P;=0;— (37) 
reduces H”’, after simplification, to 
H"’=0?+E(w,— (38) 


Here 


Since Qz is classical, and since the transformation 
(37) has a Jacobian equal to unity and leaves the 
matrix of P? unaltered, we may leave Qe un- 
changed and write 


Q=Qr(Iilelsg) Qu (Gm; Om). (39) 


The evaluation of Qy™ is more complex. If we 
define two parameters km, and €m,, such that 


Kmr = 1 — €mr = (162? /h? 
then we may write, after algebraic manipulation, 
Ou? = exp — mre?)/RT 


=exp exp 


pat) 


OmrQmQm 
=exp (—gm) exp 
r=0 4 20m 


x( (40) 
On 


the two forms being equivalent. Now for ‘‘rota- 
tional’ Mathieu levels, for “‘vibrational”’ 
levels, km-20. We may fix a transition point 
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r=ro, and write 


Qu =exp (—gm)[Gr(Om; 


+Gv (4m, dm) J; (41) 


where 


QmQm QmPm 
——— | CXp ’ 


(42) 


Qm 


This expression is exact ; the sharpness of transi- 
tion between the two types of Mathieu level is 
such that we may neglect the second exponentials 
and write, as in Eq. (6), 


r> 10 43) 
Gy= exp [- 


r<ro 


Second-order correction 


In taking account of the H”” matrices, we 
assume at once that the levels with r<7o are 
strictly vibrational, and those with r>7 strictly 
rotational. Using the limiting values (29) to (31), 
we find the second-order correction (H2)"’ 
vanishes for 7>7o (rotational) ; and for r<7o, we 
may write, using the quantum-number ?,” 


(He)**= pm? {L| Ome, Wm, ] 
+L] Ome, (mt — Wm, ]} 
= — Fypm?/2, 
where 


PP. 


The transformed diagonal sub-matrix, X”’, 
correct to the third order, is now given by 


"=H ro (H2)” r 
= 0? — (44) 


for ‘“‘vibrational” levels, and by (38) for ‘‘rota- 
tional” levels. Since the transformation (37), 
involving Um,rr, justifies us in taking @?=P? for 
vibrational levels, we may combine the matrices 
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of the first two terms of (44), writing 


Hr’ FinPm?/2 
= X(1/21;)P,?. 


Hence for the vibrational levels we must use 
Qr(IiIeI3) for the over-all rotational partition 
function, and not Qpr(IiJ2I3g). Thus we have, 
recalling (41) and (35) 


Q=exp Gv Gr] 
= Qu, 


where 


Qu =exp (— [Qn'Ge+Gy ]. (45) 


PARTITION FUNCTION: GENERAL CASE 


In extending the treatment to the case in 
which N>1, we assume that terms in 7;;? are 
negligible. The neglected terms encountered in 
the algebra do in fact tend to cancel, reducing 
the error. 


One of the most unfavorable cases is dimethyl ether, for 
which 71: =0.23, 722 =0.089, other 7;;=0. Here, fortunately, 
the error is reduced to a negligible amount by cancellation 
due to symmetry; such beneficial effects are found for any 
molecule in which the tops are symmetrically placed 
(symmetry C2, or higher). In cases of such large 7ij, 
however, it is better to calculate the needed parameters 
from Eq. (12) to (18) and the definition Qn = Fn/An, 
rather than from (21). 


The zeroth-order approximation is quickly 
written down: 


TT Qu‘ (am On). 


The first-order correction is found as it was for 
the single-top case; here the transformation (37) 
becomes 


The transformation results in the cancellation of 
terms in Ama, giving 


This is exactly similar to (38); proceeding as 
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TABLE I. 


(Cm 
Eq. (45) 
148.1 0.85 0.69 
225.35 1.16 1.04 
333.9 0.99 1.02 
367.1 0.90 0.97 


PERCENT 
fo. Deveateen 


before, we find 
IL (Gm, Om, Qm) 


with Qy™ given by (41). 

In finding the second-order correction, we 
neglect terms in Am, (the corrections from these 
terms would be of order 7;;*); proceeding as for 
the former case, we get one term for each 
“vibrating” top: 


Hy" "= pn’, 
the sum being over all tops with r<ro. We 


neglect the error caused by adding this matrix 
to that of &?, and write 


Evaluating the moment of inertia product corre- 
sponding to Hr’, we find the corresponding 
over-all rotational partition function to be 


2,), 


the product being over all ‘‘vibrating” tops. It 
follows that the complete partition function is 


Qu being given by (45). 


SYMMETRY EFFECTS 


The “quantum effects” due to symmetry have 
not yet been mentioned. Except at very low 
temperatures (below 50°K for ordinary methyl 
group rotations), these may be included just as 
in the case of over-all rotation, by summing 
over levels of all symmetries in (42) or (43) and 
then dividing the partition function by the 
Symmetry number ym. However, if we take 
Q,=1, and write G(A) for the sums taken over 


levels of symmetry A; and Az only, etc., then 
except at very low temperatures we find that 
G(A) =G(B)=G(E;*), etc. While it does not 
follow that Gr(A)=Gp(B), etc., it turns out that 
this is a good approximation. Therefore we may 
adopt the procedure discussed in the first section, 
using only the A levels, and omitting the sym- 
metry number factor. 


AN EXAMPLE: PROPYLENE 


The results of the application of this method 
to the internal hindered rotation in propylene 
have been published ;!° a few details of the calcu- 
lation will be given here by way of illustration. 
The interatomic distances were taken to be: 
C-—H, 1.09; C—C, 1.54; C=C, 1.33A; the 
angles of the CH; group were assumed to be 
tetrahedral, and all other angles were taken to 
be 120°. The value of 2, was found to be 
0.7622, appreciably different from unity; the 
“corrected moment”’ of the methyl groups was 
thus F,,=4.014(10-“ g cm?). The value 6,,=12 
(Vm=2119 cal./mole) gave the best agreement 
with thermal data. Two levels each were taken 
for Gy(A) and Gy(B), values of a», being taken 
from Ince;> these sums were found to be 7.0918 
and 7.0923, respectively, at the boiling point, 
225.35°K. It is apparent that we may take 
Gy(A)=Gy(B) without appreciable error. When 
three levels were included in Gy(A), the value 
of Sy increased by 0.01 e.u. 

A comparison of the results obtained from the 
zeroth-order expression (36) with those obtained 
from the more complete form (45) is of interest. 
At 225.35°K, the entropy difference due to the 
barrier, (S;—S), where S; is the entropy con- 
tribution for free rotation, was found to be 
1.08 e.u., from (45) ; the zeroth-order approxima- 
tion gave (S;—Sy) =1.21 e.u., the difference 
amounting to 12 percent. A similar comparison 
of the calculated heat capacity contributions is 
given in Table I. 

In concluding, I wish to express my thanks to 
Professor E. B. Wilson, Jr., for many helpful 
suggestions during the course of this study, and 
for permission to quote the calculations on 
propylene, carried out in his laboratory. 

1% B. L. Crawford, Jr., G. B. Kistiakowsky, W. W. Rice, 


A. J. Wells, and E. B. Wilson, Jr., J. Am. Chem. Soc. 61, 
2980 (1939). 
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An Experimental Determination of the Electron Affinity of Chlorine 


J. J. anp JosepH E. Mayer** 
Chemical Laboratory of The Johns Hopkins University, Baltimore, Maryland 


(Received December 8, 1939) 


The electron affinity of chlorine has been determined to be about 92.7 kcal. by measuring the 
ratio of ions to electrons and atoms leaving a hot tungsten surface. Owing to the reactivity of the 
chlorine a certain erratic behavior leads to some possible scepticism as to the accuracy of the 


result. 


INTRODUCTION 


Y using a hot (~2000°K) tungsten filament 
exposed to an atmosphere of dilute (P~10-? 
mm Hg) chlorine gas in the center of a magnetron 
tube the ratio of negative ions to electrons 
leaving the filament can be determined, whereas 
the number of neutral atoms may be calculated 
from the pressure of surrounding chlorine gas. 
These data permit the calculation’? of the 
equilibrium constant for the reaction 


Cl,+e,-—Cl,- 


2m.T 


(1) 


2) 


in which: J; and J,- are the negative ion and 
electron currents, Mci and m,.- the masses of 
chlorine atoms and electrons, Pci, the pressure 
of Cl. gas in the apparatus measured at the 
temperature Tz, and 7; is the temperature of the 
filament. 

Since, for the reaction (1) 


ci”; (3) 
in which the free energies F° are at unit pressure 
measured from free energy zero at absolute zero, 
the electron affinity Ey may be calculated from 
the data used in (2). 

The free energy of the chlorine monatomic gas 
is 


RT In K= —AF°=Eo+F 


1 4 


This a author wishes to acknowledge a grant-in-aid 
received from The Hynson, Westcott and Dunning Re- 
search Fund. 

** Now located at Columbia 

1P. P. Sutton and J. E. Mayer, J. Chem. Phys. 3, 20 


(1935). 
2H. B. Weisblatt, Johns Hopkins, Dissertation (1938). 
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The internal partition function, 


Qicr = = 4-4 
i 


(5) 


depends on the temperature. Call 


x= (6) 


and then 


In Qici=In 4+1n x. (7) 


Now if In x is expanded in a McLaurin’s series 
around 2140°K, one finds that ~ second-order 
terms are negligible. 

Finally one obtains 


logio x= T(2.5964X 10-5) +0.04077. (8) 


The resulting complete equation from (4), (3) 
and (2) for Eo is 


I; 
Eo =—4,5 738 logio Pei 
1000 I 


T 
+2 logue +1(2.5964x 10-4) +9.832] (9) 


Here Ep is the electron affinity in kcal. and Pci: is 
the pressure of chlorine in mm of Hg times a 
thousand. 


EXPERIMENTAL 


The chlorine was fractionally distilled eight 
times in a baked out all glass vacuum system. It 
was stored in 5-liter bulbs which were connected 
to the apparatus as shown in Fig. 1. The chlorine 
was released by breaking the glass tip with a 
magnetic hammer after the pressure had been 
reduced to the order of a millimeter or less by 
cooling the lower bulb in a slush of dry ice and 
petroleum ether which was under reduced 
pressure. 


ELECTRON AFFINITY OF CHLORINE 


The pressure of the chlorine in the reaction 
chamber was controlled by obtaining a steady- 
state flow through the reaction tube between the 
reservoir and a liquid-air-cooled trap. The flow 
was regulated by means of a capillary device 
shown in Fig. 2 which was between the chlorine 
reservoir and the reaction vessel. When the rod 
was inserted all of the way into the capillary, the 
pressure in the reaction tube was reduced to 1/10 
of the value obtained when the rod was free. 

The chlorine pressure was measured with a 
quartz fiber manometer designed and built by 
H. B. Weisblatt.? 

The magnetron tube in which the measure- 
ments were made is shown in Fig. 3. 

The tube and the solenoid were packed in dry 
ice in order to freeze out any reaction products 
formed by the Cl, and the metal parts. Before use 
these metal parts were baked out with an 
induction furnace. The plate, 4.7 cm diameter 
and 3.0 cm long, the guard plates, 4.7 cm 
diameter and 4.0 cm long, and the grid, 2.5 cm 
diameter and 7.5 cm long, were made of tantalum. 
The leads were made of tungsten and nickel. 

The first filament used had been previously 
operated in an atmosphere of bromine gas in 
measuring the electron affinity of bromine. After 
several days of operation the filament was 


% Porenrwmerer 


Fic. 3. 


removed and it was found that it had grown to 
about twice its original diameter, presumably 
due to the decomposition of some volatile metallic 
chloride on its surface. During these experiments 
the reaction tube had not been packed in dry ice. 

Since the temperatures are determined from 
the current voltage characteristics of the filament 
one may legitimately feel considerable scepticism 
concerning the value of these results. 

Joseph Hickman, of this laboratory, kindly 
made an x-ray powder photograph of the fila- 
ment. The resulting interplanar distances are 
given in Table I. They corresponded to none of 
the expected materials.* Before the apparatus 
was opened a yellow powder was noticed inside; 
on exposure to air it turned white. This is 
characteristic of tantalum pentachloride which 
hydrolyzes in moist air.‘ 

A second filament led to erratic results, the 


3 Hanawalt, Rinn and Frevel, Ind. Eng. Chem. Analyt. 
Ed. 10, 457 (1938). 
*H. Rose, Pogg. Ann. 99, 65 (1856). 
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calculated electron affinity being much lower 
than any probable value. 

Through the kindness of Dr. Irving Langmuir 
a pure tungsten filament of diameter 0.0199 cm 
was obtained, filament III. 

The first 31 measurements on this filament 
were taken in the course of 7.5 hours, and led to 
reasonable results, which are arranged in chrono- 
logical order in Table II. A slight tendency for 
the calculated Ey values to decrease with time is 
observable. 

Some 50 measurements taken later gave erratic 
results with calculated Eo values as low as 70 
kcal., and few results as high as 85 kcal. 


ERRORS 


This method of determination of electron 
affinity is based on only one assumption other 
than the obvious ones of accuracy in the measure- 
ments themselves: namely, that the reflection 
coefficients of ions, atoms, and electrons on the 
filament surface are unity. This means that the 
probability of a chlorine nucleus leaving the 
surface as a negative ion, or as a neutral atom, is 
independent of its past history, that is whether it 
struck the filament as molecule, ion or atom. The 
correctness of this assumption can only be tested 


by the independence of the resulting Ey on the 


temperature and Cl, pressure, which determine 
the nature of the filament surface. 

Since the results were apparently independent 
of temperature and pressure, but strongly de- 
pendent on the time for which the filament had 
been exposed to chlorine these authors believe 
that the error is to be ascribed to some difficulties 
in the data themselves, rather than on the 
reflection coefficients. 

Of the data, only the filament temperature, T;, 
appears to be in reasonable doubt. These temper- 
atures were calculated by the method suggested 
by Langmuir and Jones,’ in which the voltage 


TABLE I. d=0.7093/2 sin 6. 


5 Langmuir and Jones, Gen. Elec. Rev. 419, 310 (1927). 
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aol 1s uncertain when its value exceeds unity. Average Eo =92.7 


TABLE III. 


Filament I 
35 observations covering the range 
Ii/le Pc}, T Eo 


1.1 X10 0.5 X10%mm  1972°K 79.9 to 99.1 kcal. 
5.3 ~ 22. X10-mm 2203°K 85% within 90.1 to 95.1 


Filament observed to be coated thickly with foreign metallic coating 
after observations. 


Filament II 
74 observations covering the range 
li/Te Po, Ts Eo 
6 X10 0.24 x 10-3 1717° 75 to 106 kcal. 
to 10.4 14. X10-3 2309°K 
Calculated Eo values decreased uniformly with time at first, then rose 
abruptly at end of experiments. 


Filament III 
first 31 observations in Table II 
next 49 observations covering the range 
Ti/Te Pa, Ty Eo 
9  £0.34X10-mm 1869°K 66 to 87.5 kcal. 
to 3 X107 11.11 X10"? mm 2249°K 
Calculated Eo values decreased uniformly with time at first, then rose 
at end of the experiments. 


Ry | Ti/Ie | Pc}, X10? | T; | 
Filament III 
: 1.39 1845 94.5* 
1.39 1845 95.4* 
1.57 2034 92.5* 
1.57 |_| 2033 109.3* 
1.56 2032 91.4 
1.58 2034 90.8 
1.63 2094 91.7 
1.62 2092 91.5 
i. 2141 94.2 
1. 2126 94.0 
a, 2132 95, 
1. 2134 95 
a 2.3107 2180 94 
1. 1.8107 2179 92 
1. 3.4X107 2186 94) 
1 2187 
1 2258 
1 2263 
: 1 2320 
1 2315 
1.96 1.69 2325 
1.93 1.73 2299 
1.86 2.46 2231 
1.85 2.72 2229 
1.78 3.48 2164 
1.79 3.48 2167 
1.79 3.68 2166 
1.65 4.03 2054 
1.63 4.17 2049 4 
1.49 4.62 1980 89.2* 
3.20 1.61 0.876 
2.40 1.47 0.774 
| 2.12 1.28 0.732 
1.83 1.04 0.644 
0.600 
it 
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times cube root of current over the length of 


filament is used. This function gives the tempera- 


ture provided the surface is that of clean tungsten. 
The cooling due to the chlorine gas itself at the 
pressures used can be calculated to be negligible. 
It was assumed that a monomolecular Cl layer on 
the surface would not affect its optical properties. 

Presumably this is true, and only after con- 
siderable exposure when the filament surface has 
become covered, as noted particularly with fila- 


ment I, does the temperature calculated differ 
appreciably from the true filament temperature. 

The average of the first 31 measurements, on 
filament III, recorded in Table II are taken as 
giving a probable value for Eo. 

In Table III a brief summary of the other 
results are indicated for completeness. 

The authors wish to acknowledge the help of 
Dr. Harold B. Weisblatt in much of the experi- 
mental work. 
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The Structure of Phosphine and Related Hydrides 


D. P. STEVENSON* 
California Institute of Technology, Pasadena, California 


(Received November 20, 1939) 


It is shown that Badger’s rule when applied to the frequencies of some hydrides gives values 
for the interatomic distances which are in agreement with values deduced from the moments of 
inertia of the hydrides and with values predicted from the covalent radius suggested for 
hydrogen, 0.315A. Using the value 1.415A. for the phosphorus-hydrogen distance in phosphine 
and a reliable spectroscopic value for one of the moments of inertia of the molecule one finds the 
bond angle to be 93°. Values are suggested for the bond angles and bond distances in a number 
of hydrides whose structures have not been experimentally determined. 


N the course of calculating the thermodynamic 
functions of some compounds from spectro- 
scopic data it became apparent that the structure 
of phosphine as well as the structures of the 
hydrides of the heavier elements of the fourth, 
fifth, and sixth groups of the periodic system are 
either unknown or only partially known. Due to 
the complicated nature of the vibration-rotation 
interactions it is very difficult if not impossible 
to determine the structure of these hydrides 
from the fine structure of the rotation-vibration 
absorption bands in the infra-red. Although it 
has been proved possible to determine the inter- 
atomic distances in methane! and ammonia" by 
the electron diffraction method, it is, at present, 
not feasible to use this method for the heavier 
hydrides because of the excessive atomic scatter- 
ing of the heavy atom and the very severe 
temperature factors that arise when the vibra- 
tional frequencies become relatively small. 


* Contribution from the Gates and Crellin Laboratories 
of Chemistry, California Institute of Technology, No. 749. 

%¢V. Schomaker and D. P. Stevenson, J. Am. Chem. 
Soc. 62 (1940). 

» P. P. Debye, Physik. Zeits. 40, 404 (1939). 


Inasmuch as the structures of the hydrides are 
of interest from several points of view it seems 
worth while to report the results of an indirect 
determination of the structure of phosphine and 
of some of the heavy hydrides. 

Wright and Randall? have found 6.22 x10~-* 
g cm? for the perpendicular moments of inertia 
of phosphine. These moments of inertia are 
related to the structural parameters of the 
molecule by the equation 


a 3M 4 a 
sin? “+——_(1 —- sin? (1) 
2 M+3m 3 2 


where m is the mass of the hydrogen atom, M is 
the mass of the phosphorus atom, fro is the 
phosphorus-hydrogen distance, and a is the bond 
angle. In Table I are given the values of 7p and 
a which are compatible with the value J,=6.22 
X10-* g cm?, as well as the corresponding 
values of J;, the moment of inertia about the 


? Wright and Randall, Phys. Rev. 44, 391 (1933). 
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symmetry axis. 


(2) 


I,;=4mr,? sin? — 


The angle a is plotted against ro in Fig. 1. 

The interatomic distances found in methane,* 
ammonia,” water,‘ and hydrogen sulfide® taken 
with the covalent radii of carbon, nitrogen, 
oxygen, and sulfur as given by Pauling® give 
0.315+0.005A for the covalent radius of hydro- 
gen. The covalent radius of phosphorus is 
1.100+0.010A, hence 7o>=1.415+0.015A. Inter- 
polating by means of the graph of 7 ws. a 
(Fig. 1), one obtains a=93+2°. 

One may obtain a check on the reliability of ro 
as estimated by the use of the covalent radius 
table by the application of Badger’s rule’ to the 
known frequencies of phosphine.* Furthermore a 
measure of the accuracy of Badger’s rule may be 
obtained by comparing the interatomic distances 
given by the frequencies of methane, ammonia, 
water, and hydrogen sulfide with those found 
from the moments of inertia. In Fig. 2 a plot of 
ro vs. (1/k.)' is given for a number of diatomic 
molecules of the type 1-H, 2-H, and 3-H, the 
Arabic numerals indicating rows of the periodic 
system. The values of k, were calculated from 
the formula 


k.=0.0586M'»?, dynes per cm, 


where v is the observed frequency in wave 
numbers and M’ is the reduced mass in atomic 
weight units. The values of 79 and v were taken 
from Sponer. It was assumed that for the 


TABLE I. 


1; X104° G cm? 


6.51 
7.35 
8.26 
9.22 
10.24 


roINA 


1.399 
1.426 
1.455 
1,484 
1.515 


‘as3 33). Ginsburg and E. F. Barker, J. Chen. Phys. 3, 668 
— and R. Mecke, Zeits. f. Physik. 81, 465 
5 P. C. Cross, Phys. Rev. 46, 536 (1934) ; 47, 7 (1935). 

*L. Pauling, The Nature of the Chemical Bond (Cornell 
University —— 1939), p. 154. 

7R. M. Badger, J. Chem. Phys. 2, 128 (1934). 

8 Hibben, The Raman Effect (Rheinhold Press, 1939). 

*H. Sponer, Molekiilspektren (Springer, 1935). 


D. P. STEVENSON 


110 


| 


1.45 
To 


Fie. 1. 


1.40 


polyatomic hydrides k, has a form similar to 
that for the diatomic hydrides, 


v1) 
+p. 
=0.0586.M’>? dynes per cm, 


k.=0.0586.M’ ( 


where M’=Mm/(M+m), M is the mass of the 
central atom, m is the mass of the hydrogen 
atom, and p, is the degeneracy of the per- 
pendicular vibration. This assumption is reason- 
able for the hydrides since the interaction 
between the various modes of vibration is quite 
small. There is practically no difference between 
the force constants calculated according to Eq. 
(4) and those calculated for more complicated 
potential functions in which allowance is made 
for interaction. The values of ro obtained for 
the hydrides whose vibrational frequencies are 
known are given in column 4 of Table II with 
the values from the moments of inertia in column 
5 and those predicted from the covalent radius 
table in column 6. 

The quality of the agreement indicates that in 
general the interatomic distances predicted for 
hydrides by the covalent radius table or by 
Badger’s rule are accurate to within 0.02A. 

The values obtained for the Si-H distance in 
silane are considerably higher than the value 
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TABLE II. 


ro A (ops.)| ro A (COVALENT) 
1.09 


ke~4X102| ro A (catc.) 


a 


BERESERS 


1.42A calculated from the moment of inertia 
found by Steward and Nielsen.’® Their value is 
probably less reliable than that given by the 
covalent radius table, inasmuch as an inspection 
of their data reveals that the fine structures of 
the bands in question are less accurately known 
than those of the corresponding bands of 
methane. 

It has been suggested by Pauling® that the 
elements of the fifth and sixth groups should tend 
to form p bonds rather than hybrid s—p bonds 
because of the much greater stability of the s 
electrons in the atoms. The fact that the bond 
angles in ammonia and water are so much 


ass 3 B. Steward and H. H. Nielsen, Phys. Rev. 47, 828 
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greater than the 90° predicted for pure p bonds 
was explained on the basis of partial ionic 
character due to the large electro-negativities of 
nitrogen and oxygen. If this is so one would 
expect the bond angles to decrease when the 
nitrogen or oxygen atom is replaced by heavier 
atoms. The value 93° obtained for phosphine 
adds to the confirmation of this hypothesis 
already provided by Cross’s* value for the bond 
angle in hydrogen sulfide, 92°. One would thus 
expect that the bond angles in arsine, stibine, 
hydrogen selenide, and hydrogen telluride to be 
within 1° of 90°. 

On the basis of the conclusion reached in the 
preceding paragraphs we may predict the inter- 
atomic distances and bond angles for the hydrides 
shown in Table ITI. 

Cameron, Sears and Nielsen™ have recently 
obtained approximate values for the moments of 


TABLE III. 


MOLECULE 


inertia of hydrogen selenide by assuming the 
molecule to be almost a symmetrical top, 
I:2I3;=I/2. From the spacing of lines of the 
4.3u band of H2Se they find 3.6 10-* g cm?. 
The structure suggested in Table III leads to 
3.69 X 10~*° g cm? for this moment of inertia. 

In conclusion, the author expresses his thanks 
to Professor L. Pauling for his interest and 
helpful suggestions. 


1D. M. Cameron, W. C. Sears and H. H. Nielsen, J. 
Chem. Phys. 7, 994 (1939), 
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The Relation Between the Energy of a Hydrogen Bond and the 
Frequencies of the O—H Bands* 


RicHarp M. BADGER 


VOLUME 8 


Gates and Crellin Laboratories of Chemistry, California Institute of Technology, Pasadena, California 


T has previously been pointed out! that there 
appears to be a relation between the energy 
of a hydrogen bond and the shift in frequency of 
the O—H bands which accompanies the forma- 
tion of the bond. Recently additional data have 
been obtained in this laboratory”* which confirm 
the existence of this useful relation and it seems 
worth while to discuss them in this connection. 
The data here presented all relate to linkages 
in which a hydroxyl hydrogen is concerned, 
though the atom to which it is weakly bound 
may be oxygen, chlorine, or carbon in an aro- 
matic ring. It is surprising that the same connec- 
tion between energy and frequency shift should 
be found in these various cases but such appears 
to be the fact. For other types of bond, for ex- 
ample those involving the hydrogen of an amino 
group, a slightly different relation may well be 
expected. 

The data plotted in Fig. 1, with one exception, 
relate to cases in which a hydroxyl group which 
is initially in a relatively ideal, unperturbed 
condition enters into a hydrogen bond. The 
observed shift in frequency of the third harmonic 
O—H band which accompanies this change is 
plotted against the energy of the bond.‘ There 
is a fairly strong indication that the relation 
between frequency and energy is not linear, 
though the curve drawn is indeed somewhat 
arbitrary. The inaccuracy in the energy values 
may be at least as great as 500 cal. and the 
precision of the frequency measurements is un- 
certain. The frequencies employed are those of 
the band maxima. This is satisfactory in the case 
of narrow bands but the association bands are 
broad and unsymmetrical and it is not certain 
to what point one should measure. It seems quite 


*Contribution No. 744 from the Gates and Crellin 
Laboratories of Chemistry. 

1 Badger and Bauer, J. Chem. Phys. 5, 839 (1937). 
ow and Dickinson, J. Am. Chem. Soc. 61, 54 
3 Zumwalt and Badger, J. Am. Chem. Soc. (in press). 

‘The data for water and methyl alcohol are these of 
reference 1. 


(Received November 6, 1939) 
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certain that resolvable features which these 
bands have sometimes been supposed to possess 
are merely superposed C—H bands or other 
harmonic or combination bands. At present the 
maximum seems to be the least arbitrary point 
of reference for measurement. 

The abscissae of the points for o-chloro phenol 
and ethylene chlorohydrin are differences in 
energy of two molecular configurations in one 
of which an intramolecular hydrogen bond is 
present. The energy difference of the two states 
may be partly due to other factors beside the 
hydrogen bond, but at present it is not possible 
to make allowance for them. Similar difficulties, 
however, are encountered in the estimation of 
the energy of any hydrogen bond. For example, 
in obtaining the energy of the two bonds in the 
dimer organic acids from the dissociation energy 
one should probably make some allowance for 
van der Waals interactions. 

It is not certain that it is permissible to include 
the point for phenol in benzol solution in the 
diagram. The spectroscopic data® indicate a very 
low O—H frequency for the phenol monomer in 


600 
MeOH 


2000 4000 


6000 CALORIES 


Fic. 1. A plot of the frequency shift in the third harmonic 
O-—H bands against energy of the hydrogen bond. 


5 A. Naherniac, Ann. de physique 7, 528 (1937). 
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FREQUENCIES OF THE O-H BANDS 


benzol solution, which may be taken as evidence 
for solvation through the agency of a hydrogen 
bond. The energy quantity evaluated by Las- 
settre and Dickinson? may thus be regarded as 
the difference in energy of two hydrogen bonds. 
If the energy frequency relation is nonlinear one 
should expect the point for phenol to lie some- 
what above the curve fitting the other points. A 
similar situation may be expected for the car- 
boxylic acids. The low O—H frequency and the 
rather rigid planar COOH group of the acid 
monomers suggest some tendency for the forma- 
tion of a weak hydrogen bond in a four-membered 
ring. 

The frequency shifts plotted in the figure are 
all for the third harmonic O—H bands in the 
region 9500-11,000. The photograph region, 
beside being convenient for investigation, has 
the advantage that shifted and unshifted bands 
are better separated and more accurately meas- 
ured than in the fundamental region. It would 
no doubt be preferable to deal with frequencies 
extrapolated to zero amplitude but the data are 
at present too incomplete to permit this. In any 
case the frequency shifts in the third harmonic 
bands will usually be nearly proportional to 


6 Bauer and Badger, J. Chem. Phys. 5, 852 (1937). 
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those in the w,’s, since the convergence of most 
O—H band progressions is about the same. 

A very important exception to this is the case | 
of the carboxylic acids. By a somewhat unjusti- | 
fiable method of plotting Badger and Bauer 
brought a point for acetic acid in line with other 
data, but actually the bands of the associated 
organic acids are quite anomalous. Badger and 
Bauer have discussed the character of the O—H 
vibrations in these substances and the nature of 
the potential function. As they pointed out, an 
unsymmetrical O—H vibration is responsible for 
the infra-red association bands. Owing to the 
somewhat special character of the potential func- 
tion the succession of energy levels corresponding 
to this vibration must diverge rather than con- 
verge, as is customary. One is consequently not 
surprised to find that the progression of O—H 
bands for the double molecules diverges at about 
the same rate as other O—H progressions con- 
verge. Such anomalous cases are probably rare. 

Many applications of the relations here dis- 
cussed must be obvious and need not be dis- 
cussed. Recently spectroscopic data relating to 
hydrogen bonds have been accumulating rapidly 


and in many cases they may greatly assist in 
interpreting observations of other kinds. 
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LETTERS TO THE EDITOR 


This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section should not exceed 600 words in length and must 


The Second Virial Coefficient of Ethane 


This is to make known an error in the value of a constant 
in the second virial temperature coefficient equation for 
ethane of Eucken and Parts.! The equation is given by 
these workers as 


9.5102 

B=+59.0— (1) 
and tabulated data appended for B observed and B calcu- 
lated, over the temperature range investigated. However, 
comparing values of B calculated from Eq. (1) it will be 
found that they do not agree with those given in the table 
of Eucken and Parts. Upon inspection of the situation, it 
becomes apparent that the equation should have been 
given as 


3.9X10 9.5108 


that is, the exponent for the constant of the last term should 
be 13 instead of 12. The values in Table I, calculated for 
three of the temperatures of Eucken and Parts demonstrate 
the situation. This comparison can leave but little doubt 
that Eq. (2) above was the one employed by Eucken and 
Parts for calculating the values given in the third column of 
the above table. The use of the incorrect equation will lead 
to an error of ca. 8 percent at the lowest temperature and 
of 2.6 percent at the highest, in B. 


TABLE I. 


VALUES OF —B IN CC PER MOLE 
From E. & P. TABLE CALC. FROM ABOVE EQUATIONS 
Oss. CALc. (1) (2) 


4aB/dT, CC/MOLE/DEG. 
Ea. (1) (2) 


d?B/dT?, ccC/MOLE/DEG.? 
(1) ( 


—0.1058 
—0.0207 
—0.00629 


reach the office of the Managing Editor not later that the 15th 
of the month preceding that of the issue in which the letter 
is to appear. No proof will be sent to the authors. The usual 
publication charge ($3.00 per page) will not be made and no 
reprints will be furnished free. 


It is of importance to ascertain the resulting errors in 
dB/dTand in @B/dT? that will occur from employing the 
incorrect equation. The results of some calculations are 
given in Table II. It will be observed that appreciable 
error will be made in the derivatives, especially at the 
lower temperatures. 

Mr. W. W. Rice of these laboratories has kindly checked 
the calculations in the above two tables. 

A search of volumes of Zeitschrift fiir physikalische 
Chemie later than B19 failed to reveal any recorded cor- 
rection of this error. 


Epwin E. 
Mallinckrodt Chemical Laboratories, 
Harvard University, 
Cambridge, Massachusetts, 
February 6, 1940. 


1 Eucken and Parts, Zeits. f. physik. Chemie B20, 184 (1933). 


The Crystal Structure of Sodium Sulfate III 


Of the five modifications of anhydrous sodium sulfate,'? 
thenardite is the only form that has been examined 
structurally. To help understand the nature of the temper- 
ature transition from thenardite to NazSO, III, the struc- 
ture of the latter modification has been determined. 
Na2SO, III belongs to the orthorhombic space group 
Pb nn. The unit cell contains 4Na2SO, and has the fol- 
lowing dimensions at room temperature: - 


a=5.59+0.02A, b=8.93+0.02A, c=6.98+0.02A. 


The tetrahedral SO," groups are located at +(4, v, }) and 
+(4, }-—v, 3) with v=0.097+0.001. The two sets of 
sodium ions have the following parameters: +(}, 1, 2); 
+(4, 4-0, 3) with v:=0.43540.010; and +(u, }, 9), 
+(3—u, 4, 4) with u=—}. The sixteen oxygens occupy 
two sets of general positions +(x, y, z), +(x, 
4-2), 2+4) with m1, % 
=0.030+0.002, 0.197+0.003, and x2, ys, 9, 
0.072 +0.002. 


Lupo K. FREVEL 


The Dow Chemical Company, 
Midland, Michigan, 
january 23, 1940. 


1F. C. Kracek and R. E. Gibson, J. Phys. Chem. 33, 1304 (1929). 
2F. C. Kracek and C. J. Ksanda, J. Phys. Chem. 34, 1741 (1930). 
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4 
191.86 498 499 458.4 498.5 
224.50 354 353 335.5 352.7 
273.20 227.5 228 222.9 228.8 
TABLE II. 
272.07 1.783 1.90 —0.0179 
| 
| 
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Zeeman Effect From Liquids, Solutions of 
Salts of Europium 


Among liquids, the solutions of the salts of europium 
appear to be unique in the sharpness of their absorption 
spectra. The absorption lines cluster in groups which owe 
their origin to the splitting of the energy levels of the 
europium ions by the electric fields of their environment. 
Should the symmetry of the fields about the ions be 
reduced, as by superimposing an external field, a further 
splitting is to be expected in general and also an increase 
in the number of lines within the groups. It is in the study 
of the structures of the groups that we have turned to the 
effect of a magnetic field on the spectra. 

Zeeman effects from crystals at low temperatures have 
indicated that the lines must be very sharp before the 
effect of the magnetic field can be resolved. It has actually 
been found, notably from aqueous solutions of europium 
nitrate, that some of the Zeeman components at the fields 
available could be separated with an instrument of rather 
modest resolving power and dispersion, a three-prism 
Steinheil spectrograph having a dispersion of about 13A 
per mm at 4650A. 

The group of absorption lines at 4650A which, in the 
absence of the field, had consisted of three lines at 4645.0A, 
4649.;A, 4654.sA increased to five, when the radiation 
passed through the solution transverse to the field. How- 
ever, the existence of one of these lines was somewhat 
uncertain. As the magnetic field was increased, a compo- 
nent was observed to split off the line at 4649.s toward 
long wave-lengths until at 33,000 gauss, the separation 
from the undisplaced component amounted to 7 cm 
+1.5 cm™ or about six to eight times the normal Larmor 
precession. Another component also appeared to separate 
off from the same line toward shorter wave-lengths, but it 
was much less sensitive to the magnetic field and indeed 
the displacement was very small relative to the dispersion 
of the spectrograph. However, its presence seemed to be 
confirmed on examining the polarization of the spectra. 
The sensitive component represented absorption of radia- 
tion having its electric vector at right angles to the mag- 
netic field; the insensitive component was polarized in the 
same sense. However, the undisplaced component remained 
unpolarized. The line at 4645.oA displaced by the field to 
4644.;A also remained unpolarized. 

The effect of a magnetic field on the spectrum from a 
solution of europium chloride was altogether different from 
what occurred with the nitrate. 

At some later occasion, a more. detailed description of 
the spectra arising in the magnetic field will be given. 

We gratefully acknowledge the support of the Penrose 
Fund of the American Philosophical Society. 
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Raman Effect of Symmetric Difluorotetrachloroethane 


We have observed twenty-two Raman shifts in sym- 
metric difluorotetrachloroethane. The results are shown in 
Table I. Three exposures were made on Eastman Spectro- 
scopic Plates Type 1-J using our usual set-up of a three- 
prism Steinheil spectroscope and eight concentric mercury 
discharge tubes as a source of radiation. Exposure times 
employed were 22 and 42 hours with a slit width of 0.05 
mm, and 46 hours with a slit at 0.1 mm. The spectroscope 
was thermostated at 33.4+2°C. The average temperature 
of the sample was about 40°C. 

A separate exposure was necessary in order to obtain the 
line at 81.4 cm™ since the region of the plate close to 4358 
and 4046 is blocked off by dihedral light stops to prevent 
overexposure and consequent halation of the strong lines. 


TABLE I. Raman lines of symmetric tetrachlorodifluoroethane. 


ESTIMATED 
Av RELATIVE % MEAN EXCITING No. or 

cm"! INTENSITY DEVIATION LINES READINGS 
81.4 5 0.25 c 2 
165.5 3 36 c 4 
181.2 5 .38 c 6 
231.0 7 30 a c 12 
264.3 7 45 a c 12 
297.2 5 A7 a c 12 
329.1 7 21 a c 12 
402.1 3 32 a c 5 
421.2 10 .26 a é 12 
444.7 8 .22 a c 10 
455.5 2 .50 a c 7 
523.8 8 34 a c 12 
578.2 6 .27 a c 10 
625.9 8 .23 a c 12 
778.9 0 .25 c 2 
854.8 4 19 a c 8 
906.8 2 40 a c 6 
1027.2 4 .20 a c & 
1046.0 0 .26 a c 7 
1086.0 0 .23 c 2 
1143.0 2 17 ¢ 3 
1167.5 3 08 3 


a =4046, c =4358A. 


An exposure of 13.5 hours with a slit width of 0.1 mm was 
made with the stops out. The line was found between the 
halated excitation line and the circular halo which is 
formed about 100 cm~ from it. 

The material used was furnished by the E. I. duPont 
de Nemours & Company and was specified as ‘‘carefully 
fractionated a number of times, practically 100 percent 
pure.” We wish to thank Dr. A. F. Benning of the Jackson 
Laboratory for the loan of this substance. Further con- 
siderations of these results will be published later. 

Gero. GLOCKLER 
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University of Minnesota, 
Minneapolis, Minnesota, 
February 1, 1940 


